Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 300,698 A→C 7.8% A79A (GCT→GCG)  paoC ← aldehyde dehydrogenase: molybdenum cofactor‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,6980AC 7.8% 470.7 / 29.2 167A79A (GCT→GCG) paoCaldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (87/66);  new base C (8/5);  total (96/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 170 total aligned reads displayed.

AGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTT  >  NC_000913/300549‑300847
                                                                                                                                                     |                                                                                                                                                      
aGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACa                                                                                                                                                        >  2:1490896/1‑150 (MQ=255)
 gCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAg                                                                                                                                                       <  1:775111/150‑1 (MQ=255)
  cGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTTGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGc                                                                                                                                                      >  2:1370341/1‑150 (MQ=255)
     ggCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAg                                                                                                                                                   >  1:2108555/1‑150 (MQ=255)
     ggCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAg                                                                                                                                                   >  1:92073/1‑150 (MQ=255)
     ggCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAg                                                                                                                                                   >  2:1160467/1‑150 (MQ=255)
        cgccgcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑Ac                                                                                                                                               >  1:612572/1‑150 (MQ=255)
         gccgcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCTGTCATTACCGCCAGT‑aac                                                                                                                                              <  1:15645/150‑3 (MQ=255)
         gccgcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCG‑TACTGGCTGTCATTACCGCCAGT‑aacg                                                                                                                                             <  1:1017451/150‑4 (MQ=255)
          ccgcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGc                                                                                                                                             >  1:2186738/1‑150 (MQ=255)
            gcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCTGTCATTACCGCCAGT‑aacgcc                                                                                                                                           <  2:28269/150‑6 (MQ=255)
                  gcCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCt                                                                                                                                     >  2:1233993/1‑150 (MQ=255)
                  gcCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCt                                                                                                                                     <  2:667457/150‑1 (MQ=255)
                   ccTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCtt                                                                                                                                    >  1:2192541/1‑150 (MQ=255)
                      gTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCttttt                                                                                                                                 >  2:1263148/1‑150 (MQ=255)
                      gTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCttttt                                                                                                                                 <  2:2375365/150‑1 (MQ=255)
                        tCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTgc                                                                                                                               >  1:1440778/1‑150 (MQ=255)
                        tCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTgc                                                                                                                               <  2:2147940/150‑1 (MQ=255)
                           aaGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGcgcg                                                                                                                            <  1:2472865/150‑1 (MQ=255)
                           aaGGGCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTGTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGcgcg                                                                                                                            <  1:1363277/150‑1 (MQ=255)
                              gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCgcggcg                                                                                                                         >  1:1629689/1‑150 (MQ=255)
                                ctcGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTc                                                                                                                       >  2:653828/1‑150 (MQ=255)
                                 tcGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTcc                                                                                                                      >  2:838883/1‑150 (MQ=255)
                                   ggCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGt                                                                                                                    <  2:510966/150‑1 (MQ=255)
                                    gCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTa                                                                                                                   <  2:173511/150‑1 (MQ=255)
                                        cTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCaa                                                                                                               <  2:271486/150‑1 (MQ=255)
                                          aCCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAgg                                                                                                             <  2:1791665/150‑1 (MQ=255)
                                          aCCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAgg                                                                                                             >  1:387427/1‑150 (MQ=255)
                                           ccAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAggg                                                                                                            >  2:1813921/1‑150 (MQ=255)
                                              gcgcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCgg                                                                                                         <  1:2012322/150‑1 (MQ=255)
                                                gcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTg                                                                                                       <  2:2108555/150‑1 (MQ=255)
                                                 cAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGa                                                                                                      <  2:1258183/150‑1 (MQ=255)
                                                   aTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAgg                                                                                                    <  1:910759/150‑1 (MQ=255)
                                                      gCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGt                                                                                                 >  1:224842/1‑150 (MQ=255)
                                                            tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTg                                                                                           >  2:1277333/1‑150 (MQ=255)
                                                            tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTg                                                                                           <  1:1682139/150‑1 (MQ=255)
                                                             gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTgg                                                                                          >  1:549155/1‑150 (MQ=255)
                                                              atAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGc                                                                                         <  2:92073/150‑1 (MQ=255)
                                                              atAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGc                                                                                         >  2:844759/1‑150 (MQ=255)
                                                              atAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGc                                                                                         <  1:1518989/150‑1 (MQ=255)
                                                               tAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCa                                                                                        >  1:729514/1‑150 (MQ=255)
                                                                 gTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAAt                                                                                      >  1:2123340/1‑150 (MQ=255)
                                                                   gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgg                                                                                    >  2:1697836/1‑150 (MQ=255)
                                                                   gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgg                                                                                    >  1:1435081/1‑150 (MQ=255)
                                                                   gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgg                                                                                    >  1:1093365/1‑150 (MQ=255)
                                                                      cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCgg                                                                                 >  1:2428744/1‑150 (MQ=255)
                                                                        aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGaa                                                                               >  1:1475785/1‑150 (MQ=255)
                                                                         tAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAc                                                                              >  1:1381066/1‑150 (MQ=255)
                                                                           gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcg                                                                            <  1:2379934/150‑1 (MQ=255)
                                                                              gggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgacg                                                                         <  1:1492328/150‑1 (MQ=255)
                                                                                gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGat                                                                       <  2:625722/150‑1 (MQ=255)
                                                                                  cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat                                                                     >  2:1622479/1‑150 (MQ=255)
                                                                                   gccTAACAGCCTGGCGGTGTTTTTGTCTCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGAtata                                                                    <  1:869311/150‑1 (MQ=255)
                                                                                   gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGAtata                                                                    <  1:1233993/150‑1 (MQ=255)
                                                                                    ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAg                                                                   >  1:441081/1‑150 (MQ=255)
                                                                                     cTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGc                                                                  <  1:538424/150‑1 (MQ=255)
                                                                                      tAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGCACCGACGATATAGcc                                                                 >  2:1587897/1‑150 (MQ=255)
                                                                                        aCAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCAt                                                               >  1:2247878/1‑150 (MQ=255)
                                                                                        aCAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCAt                                                               <  2:724080/150‑1 (MQ=255)
                                                                                         cAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATa                                                              >  2:1041805/1‑150 (MQ=255)
                                                                                          aGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAg                                                             >  2:1499223/1‑150 (MQ=255)
                                                                                           gCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAgg                                                            >  1:1674050/1‑150 (MQ=255)
                                                                                            ccTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggc                                                           <  1:1213492/150‑1 (MQ=255)
                                                                                              tGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg                                                         <  1:1160467/150‑1 (MQ=255)
                                                                                                    tgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgg                                                   <  1:2427113/150‑1 (MQ=255)
                                                                                                    tgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgg                                                   <  1:601883/150‑1 (MQ=255)
                                                                                                          tGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGctt                                             >  2:2056381/1‑150 (MQ=255)
                                                                                                           gTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttc                                            <  1:322451/150‑1 (MQ=255)
                                                                                                             cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGGAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttctt                                          >  1:1773367/1‑150 (MQ=255)
                                                                                                             cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttctt                                          <  1:1274117/150‑1 (MQ=255)
                                                                                                                 tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATg                                      >  2:25601/1‑150 (MQ=255)
                                                                                                                   tGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGcc                                    <  2:1762538/150‑1 (MQ=255)
                                                                                                                      cGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCAtt                                 >  2:55938/1‑150 (MQ=255)
                                                                                                                       gAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTc                                >  2:1887993/1‑150 (MQ=255)
                                                                                                                         gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGt                              <  2:1435081/150‑1 (MQ=255)
                                                                                                                         gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGt                              <  2:1955895/150‑1 (MQ=255)
                                                                                                                          tGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTa                             <  1:1803905/150‑1 (MQ=255)
                                                                                                                            ccccGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                           >  1:1189666/1‑150 (MQ=255)
                                                                                                                            ccccGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                           >  2:1363882/1‑150 (MQ=255)
                                                                                                                            ccccGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                           <  2:1233096/150‑1 (MQ=255)
                                                                                                                              ccgggcgtACTGGCTGTCATTACCGCCAGTAACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGTTATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGc                          >  1:1077393/8‑150 (MQ=38)
                                                                                                                               cGGCGTTACTGGCAGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGt                        >  2:1902619/1‑150 (MQ=255)
                                                                                                                                   gTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcg                    >  2:1709106/1‑150 (MQ=255)
                                                                                                                                   cgtACTGGCTGTCATTACCGCCAGTAACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGc                     >  2:1177959/3‑150 (MQ=255)
                                                                                                                                     tACTGGCTGTCATTACCGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtg                  >  2:398077/1‑150 (MQ=255)
                                                                                                                                       cTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgcc                >  2:519881/1‑150 (MQ=255)
                                                                                                                                         ggCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATg                                      >  2:153001/1‑126 (MQ=255)
                                                                                                                                          gCGGTAATGACGGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccgtg             >  1:877026/1‑150 (MQ=255)
                                                                                                                                           cGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccgtgc            >  1:305143/1‑150 (MQ=255)
                                                                                                                                            ggTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCgtgccgtgcc           >  1:128982/1‑150 (MQ=255)
                                                                                                                                             gTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGtgccgtgccg          <  1:1813921/150‑1 (MQ=255)
                                                                                                                                             gTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGtgccgtgccg          >  1:378490/1‑150 (MQ=255)
                                                                                                                                              tAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGTCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCgg         <  2:2247878/150‑1 (MQ=255)
                                                                                                                                              tAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTGGCAATCGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCTTAGCGTGCCGTGCCgg         <  2:940791/150‑1 (MQ=255)
                                                                                                                                               aaTGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGt        >  2:1594157/1‑150 (MQ=255)
                                                                                                                                                aTGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGtt       >  1:1245191/1‑149 (MQ=255)
                                                                                                                                                aTGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGATTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTa       >  2:1794787/1‑150 (MQ=255)
                                                                                                                                                 tGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg      >  2:2380505/1‑150 (MQ=255)
                                                                                                                                                 tGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg      <  1:1370341/150‑1 (MQ=255)
                                                                                                                                                     aGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGtttt  <  1:439120/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGT‑ACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTT  >  NC_000913/300549‑300847

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑