Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,283,465 G→A P480S (CCG→TCG)  actP ← acetate/glycolate:cation symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,283,4650GA95.8% 138.7 / ‑4.1 48P480S (CCG→TCG) actPacetate/glycolate:cation symporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (29/17);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

TCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATACCCGCAGCTCTTCACGTTC  >  NC_000913/4283221‑4283711
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
ccccagggcAGACGACCGGGGAAGGGGGGAGATTAAGGCGCTCGGCCGTTCGAAACGCCAAAGCCGGTCGGGGGGGGGATAAACTGCGCGCGGAACAGTTCACGCTCAGGCGCGCCTTCCGCTGATTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGGACCCAAATCGTCGAGCCGAg                                                                                                                                                                                                                                                  <  2:416375/243‑1 (MQ=255)
   ccGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGat                                                                                                                                                                                                                                               >  1:182355/1‑251 (MQ=255)
                 ggggAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTATGGGAGCGGATAAACTGCGCGAGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCCCGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCAAAATCGTCGAGCCGAGGATAATCAGTACTACTg                                                                                                                                                                                                                                 >  2:436129/1‑251 (MQ=255)
                          gcgagaTTCATGCGCGAGGCCTTGCTCAACGAAAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTACGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCACCTGGTATTCATACGGGAAGATTGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACAACTGCGGTAATc                                                                                                                                                                                                                         >  2:91917/3‑250 (MQ=255)
                                                   tCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCCCGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCAACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATcgc                                                                                                                                                                                               >  1:427806/1‑251 (MQ=255)
                                                     aaCGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGAACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCgcgc                                                                                                                                                                                             >  1:317728/1‑251 (MQ=255)
                                                              gCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTgg                                                                                                                                                                                    >  1:68387/1‑251 (MQ=255)
                                                                         ggCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTAATCGGTGGCCGAGAAGAACCATATGCCGATGAAGGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGGAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGGCAGGTTCGAc                                                                                                                                                                         <  2:427806/250‑1 (MQ=255)
                                                                            ggATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGATGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGAcc                                                                                                                                                                        >  1:223169/1‑249 (MQ=255)
                                                                                     gcgcgcGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAg                                                                                                                                                               <  1:344598/249‑1 (MQ=255)
                                                                                           ggAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGATGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGa                                                                                                                                                       <  2:223169/251‑1 (MQ=255)
                                                                                              aCAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATg                                                                                                                                                    >  1:270492/1‑251 (MQ=255)
                                                                                                       gctcacgcGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGTCCAGTACATAGAAAGCAGAATGATCGGGa                                                                                                                                             >  2:110569/1‑249 (MQ=255)
                                                                                                        ctcacgcGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGGAAGCAGAATGATCGGGAAGt                                                                                                                                          >  1:483150/1‑251 (MQ=255)
                                                                                                               gcgcCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCCTCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGGACATAGAAAGCAGAATGATCGGGAAGTTACAGCt                                                                                                                                   >  2:275793/1‑251 (MQ=255)
                                                                                                                cgcCTTCCGCTGAGTTATCGGTGGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTc                                                                                                                                  <  2:182355/251‑1 (MQ=255)
                                                                                                                 gcCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTAc                                                                                                                                       >  1:303915/1‑245 (MQ=255)
                                                                                                                                  cGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCCAACGcc                                                                                                                >  2:77044/1‑251 (MQ=255)
                                                                                                                                  cGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGcc                                                                                                                <  1:119518/251‑1 (MQ=255)
                                                                                                                                       gCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGGATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGGGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCAAGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAACCGCCAGccc                                                                                                           >  2:375676/1‑251 (MQ=255)
                                                                                                                                             aagaaCCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATcg                                                                                                                                                                                                <  2:90050/160‑1 (MQ=255)
                                                                                                                                             aagaaCCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATcg                                                                                                                                                                                                >  1:90050/1‑160 (MQ=255)
                                                                                                                                                          cGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATACAGGCGAGGTTc                                                                                        >  1:235336/1‑251 (MQ=255)
                                                                                                                                                           gAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCt                                                                                       >  1:374098/1‑251 (MQ=255)
                                                                                                                                                            aGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGaaa                                                                                                                                                            <  2:99813/181‑1 (MQ=255)
                                                                                                                                                            aGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGaaa                                                                                                                                                            >  1:99813/1‑181 (MQ=255)
                                                                                                                                                                 tGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGc                                                                                               <  1:169924/237‑1 (MQ=255)
                                                                                                                                                                 tGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGc                                                                                               >  2:169924/1‑237 (MQ=255)
                                                                                                                                                                    cACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGGTCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCa                                                                               >  2:400930/1‑250 (MQ=255)
                                                                                                                                                                          gATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATcgc                                                                                                                                                                                               <  1:102628/132‑1 (MQ=255)
                                                                                                                                                                          gATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTCCAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGc                                                                         >  2:423137/1‑250 (MQ=255)
                                                                                                                                                                          gATAGAGAACAGCGCCGGGGATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATcgc                                                                                                                                                                                               >  2:102628/1‑132 (MQ=255)
                                                                                                                                                                                        ccGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATc                                                          >  2:286215/1‑251 (MQ=255)
                                                                                                                                                                                               ttCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAACCGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATc                                                    >  1:434068/1‑250 (MQ=255)
                                                                                                                                                                                                     cGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCg                                              >  1:429318/1‑250 (MQ=255)
                                                                                                                                                                                                              ggCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTa                                    <  1:148484/251‑1 (MQ=255)
                                                                                                                                                                                                              ggCTTTTTCGTGACCAAGGATCTGTACCCAAATCGGCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGt                                     <  1:286215/250‑1 (MQ=255)
                                                                                                                                                                                                                 tttttCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACgg                                 >  1:157536/1‑251 (MQ=255)
                                                                                                                                                                                                                   tttCGTGACCAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGt                                <  1:275039/250‑1 (MQ=255)
                                                                                                                                                                                                                          cccatggATCTGTACCCAAGCCGTCGGGACGAGGATCATCAGTACTACTGCGGTAAGCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGAACATAGAAAGCAGAATGATCGGGAAGTTACAGCGCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTg                        <  2:91023/246‑1 (MQ=255)
                                                                                                                                                                                                                            cAAGGATCTGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTCCTGATTCTCAACCAGCACGCCGCGGAGAATCGCAATCACGCCGAGGATAAGTACGGTGATTTTTggt                      >  2:71702/1‑249 (MQ=255)
                                                                                                                                                                                                                               ggATCTGTACCCAAATCGTCGAGCCGGGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGTATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATAcc                   <  1:375676/251‑1 (MQ=255)
                                                                                                                                                                                                                               agatatGTACACAAATCGTCGAGCCGCGGATCATCCGGACTCCTGCAGTAAACAGCCACAGCCAGCAACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGCTACAGCTCGCCGCGGTGGCAAACGCCAGCCCCACCATACCGGCGATGTTCTGATTCTCAAACAGCACGCCCAGGGCAATCGCCAACACGCCGGGGATCACTACGGGGAGTTTTGATAcc                   >  2:207709/6‑251 (MQ=255)
                                                                                                                                                                                                                                gATCTGTACCCAAATCGTCGAGCCGAGGATCATCCGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATAccc                  >  1:354531/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tGTACCCAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGCTGTTCTGATTCCCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCCGGACGGTGATTTTTGAAACCCGccg              >  2:159298/1‑249 (MQ=255)
                                                                                                                                                                                                                                          cAAATCTTCGAGCCGAGGATCATCAGTACTACTGCTGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGACAGTTTCGCCCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTTCGGTGATTTTTGATTCCCGCAGCTCTTc        <  1:297703/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aatggtcgAGCCTAGGATCATCAGTACTACTGCGTGAATCATCCCATGCCAGCCCCCCATCATCGCGACACGCTGGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTGACAGCACGCCGGGATGGAAAACGCCAGCCCCACCATAAATGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCTCAATCACGCCGAGGATCAGTACGGTGATTTGTGATACCCGCAGctctcctt      <  2:268141/247‑5 (MQ=255)
                                                                                                                                                                                                                                                gtcgAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATACCCGCAGCTCTTCACGTTc  <  1:77044/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
TCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACCCAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCAAACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATACCCGCAGCTCTTCACGTTC  >  NC_000913/4283221‑4283711

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: