Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,001 C→A 6.6% intergenic (‑100/+167) ygfF ← / ← gcvP putative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,046,0010CA 6.6% 356.4 / 12.2 122intergenic (‑100/+167)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (44/70);  new base A (4/4);  total (48/74)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-02

Only 100 of 122 total aligned reads displayed.

TTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAG  >  NC_000913/3045853‑3046147
                                                                                                                                                    |                                                                                                                                                  
ttGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACACTTTTCCCCCTCGCCCTTTCa                                                                                                                                                   <  2:1432976/150‑1 (MQ=255)
 tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAg                                                                                                                                                  >  1:1185254/1‑150 (MQ=255)
 tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTAAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAg                                                                                                                                                  >  1:127871/1‑150 (MQ=255)
  gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAgg                                                                                                                                                 >  1:704093/1‑150 (MQ=255)
   cccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAggg                                                                                                                                                >  2:82907/1‑150 (MQ=255)
   cccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACATATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAggg                                                                                                                                                >  1:372817/1‑150 (MQ=255)
    ccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGga                                                                                                                                               <  2:1074141/150‑1 (MQ=255)
      gcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGgaga                                                                                                                                             <  2:586327/150‑1 (MQ=255)
      gcccaATGCCTCGACTGCCACCAGGCACTAGTGCTATAGCCATTTGTTTTTACGAGAAAAATTAATTCGTTGAGGTATCTTTGAGATTGCTTACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGgaga                                                                                                                                             <  1:1036818/150‑1 (MQ=255)
       cccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGagag                                                                                                                                            >  1:103711/1‑150 (MQ=255)
       cccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGagag                                                                                                                                            <  1:1422297/150‑1 (MQ=255)
         cgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTTGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCTCCCTGAAAGGGCGAggg                                                                                                                                          <  2:149438/150‑1 (MQ=255)
         cgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCTCCCTGAAAGGGCGAggg                                                                                                                                          <  1:262217/150‑1 (MQ=255)
          gATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCTCCCTGAAAGGGCGAgggg                                                                                                                                         <  2:592575/150‑2 (MQ=255)
          gATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGc                                                                                                                                         <  1:207386/150‑1 (MQ=255)
           aTGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCTCCCTGAAAGGGCGAggggg                                                                                                                                        <  2:1571910/150‑3 (MQ=255)
            tGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCg                                                                                                                                       <  2:1566662/150‑1 (MQ=255)
              ccGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCggg                                                                                                                                     >  1:1438211/1‑150 (MQ=255)
               cgcgACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCgggg                                                                                                                                    <  2:1617860/150‑1 (MQ=255)
               cgcgACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCgggg                                                                                                                                    <  1:579004/150‑1 (MQ=255)
                 cgACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCTCGCTTTTCCACCTCGCCCTTTCCGGGAGAGGGCCGGGGTg                                                                                                                                  >  1:269777/1‑150 (MQ=255)
                   aCTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAg                                                                                                                                <  2:1163078/150‑1 (MQ=255)
                   aCTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAg                                                                                                                                >  1:1399403/1‑150 (MQ=255)
                     tGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTACCGAGAAAAATTAATTCTTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTACCCCTCGCCCTTTCAGTGAGAACGCCGGGGTGAggg                                                                                                                              >  1:1248121/1‑150 (MQ=255)
                      gCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGt                                                                                                                             >  1:827220/1‑150 (MQ=255)
                        caccaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTaa                                                                                                                           <  1:1537190/150‑1 (MQ=255)
                         accaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTaaa                                                                                                                          >  2:626842/1‑150 (MQ=255)
                         aacaGTCACTAGTGCTATATCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGTCTTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTaaa                                                                                                                          <  2:819678/148‑1 (MQ=255)
                          ccaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAat                                                                                                                         <  1:368094/150‑1 (MQ=255)
                            aGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAatat                                                                                                                       >  2:36482/1‑150 (MQ=255)
                               aacaAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcg                                                                                                                    <  2:1607280/149‑1 (MQ=255)
                                 caAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcgcc                                                                                                                  >  1:156331/1‑149 (MQ=255)
                                   aGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGcc                                                                                                                <  2:1683650/150‑1 (MQ=255)
                                    gTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCa                                                                                                               >  2:1425792/1‑150 (MQ=255)
                                     tGCTTTAGCCATTTGTTTTTCCTAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTGGGCACGCTGTTCCCCCTCGCCCTTTCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAg                                                                                                              <  1:1305974/150‑1 (MQ=255)
                                     tGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAg                                                                                                              <  1:1322178/150‑1 (MQ=255)
                                          taGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGTGTAAATATTCGCGCCAGTGCCg                                                                                                         <  1:1429308/150‑1 (MQ=255)
                                             ccATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGcc                                                                                                      <  1:567399/150‑1 (MQ=255)
                                             ccAATTGTTTTTCCGAGAAAAATTAATTCGGTGACGTATCTTTGAGATTGCTAACAACCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGcc                                                                                                      <  2:306947/150‑1 (MQ=255)
                                               aTTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTg                                                                                                    >  2:870394/1‑150 (MQ=255)
                                                  tGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTTACTAACTGGCGCATGCTTCATTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTc                                                                                                 <  1:1258636/150‑1 (MQ=255)
                                                    tttttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTccc                                                                                               <  1:1630737/150‑1 (MQ=255)
                                                      tttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAGCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                                             <  1:1169911/150‑1 (MQ=255)
                                                          cGAGAAAAATTATTTCGTTGACGTATCTTTGAGTTTGCTAACAATCTGGCTCAATCTTCTATGTTGGCACGCTGTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                         <  2:720157/150‑1 (MQ=255)
                                                          cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                         >  1:318772/1‑150 (MQ=255)
                                                          cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCATCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                         <  2:606925/150‑1 (MQ=255)
                                                          cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGATCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCTCGCAAGTGCCGGCCTGTTCCCCACAcc                                                                                         >  2:877460/1‑150 (MQ=255)
                                                          cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTCCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                         <  2:1399403/150‑1 (MQ=255)
                                                           gagaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAccc                                                                                        <  2:804130/150‑1 (MQ=255)
                                                             gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTa                                                                                      >  1:1658611/1‑150 (MQ=255)
                                                                   ttaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctc                                                                                <  1:601485/150‑1 (MQ=255)
                                                                      attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                             <  1:381814/150‑1 (MQ=255)
                                                                      attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                             <  2:514693/150‑1 (MQ=255)
                                                                         cGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGTCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                                                          >  2:232148/1‑150 (MQ=255)
                                                                         cGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                                                          <  2:23355/150‑1 (MQ=255)
                                                                           ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                        <  1:1564530/150‑1 (MQ=255)
                                                                                gTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACACTAACCCTCTACCCAAAGGGGc                                                                   >  1:842958/1‑150 (MQ=255)
                                                                                 tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                                  <  2:44453/150‑1 (MQ=255)
                                                                                  aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCATTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                 <  1:1327991/150‑1 (MQ=255)
                                                                                     tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAggg                                                              <  1:626842/150‑1 (MQ=37)
                                                                                      ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgggg                                                             <  2:827220/150‑1 (MQ=255)
                                                                                       tGAGATTGCTAACAAACTCGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                            <  1:1581155/150‑1 (MQ=255)
                                                                                         agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgggttag                                                          <  1:980665/150‑5 (MQ=37)
                                                                                            ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTc                                                       <  2:1571291/150‑1 (MQ=255)
                                                                                            ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTc                                                       <  2:826801/150‑1 (MQ=255)
                                                                                             tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCTGTCt                                                      >  2:1448137/1‑150 (MQ=255)
                                                                                              gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGGGAGGGTAAATATTCGCGCCAGTGCCGGCATGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTg                                                     <  2:858942/150‑1 (MQ=255)
                                                                                               cTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGa                                                    <  1:1456596/150‑1 (MQ=255)
                                                                                                tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                   <  1:82907/150‑1 (MQ=255)
                                                                                                 aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                  >  1:633865/1‑150 (MQ=255)
                                                                                                    aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGGCTGTCTGAGCACt                                               >  1:218198/1‑150 (MQ=255)
                                                                                                         ggttCAATCTTCAAGGTGGGCAGGCTTTTCCCCCTGGCCCTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                         <  2:198395/147‑1 (MQ=25)
                                                                                                          gttCAATCTTCAATGTTGGCTAGCTTTTCCCCGTCTCCCTTTCAGGGAGAGGGCCGGGGTGAGTGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCGCGAAAGGGGCGAGGGGAGTGTCTGAGCACTTTTGTa                                         <  2:1248121/148‑1 (MQ=25)
                                                                                                            tCAATCTTCAATGTTGGCACGGTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGGCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACt                                       <  2:704093/150‑1 (MQ=255)
                                                                                                               aTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTg                                    <  2:775930/150‑1 (MQ=255)
                                                                                                                 cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGt                                          <  1:252929/142‑1 (MQ=255)
                                                                                                                 cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGt                                          >  2:252929/1‑142 (MQ=255)
                                                                                                                  ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCa                                 >  2:1195379/1‑150 (MQ=255)
                                                                                                                      aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                             <  1:683557/150‑1 (MQ=255)
                                                                                                                      aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                             <  2:719102/150‑1 (MQ=255)
                                                                                                                      aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                             <  2:588361/150‑1 (MQ=255)
                                                                                                                      aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                             <  2:103711/150‑1 (MQ=255)
                                                                                                                        gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCTTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                           <  1:1669821/150‑1 (MQ=255)
                                                                                                                        gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                           <  1:36482/150‑1 (MQ=255)
                                                                                                                            gCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaa                       >  2:926482/1‑150 (MQ=255)
                                                                                                                               cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg                    >  2:473256/1‑150 (MQ=255)
                                                                                                                                    ccggcccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGcc               >  2:1167512/5‑150 (MQ=25)
                                                                                                                                        cccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg           >  1:526649/1‑150 (MQ=25)
                                                                                                                                        cccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg           >  2:1547790/1‑150 (MQ=25)
                                                                                                                                        cccTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAAAATTCGCGCCAGAGCCGGCCTGTTCCCCTCACCCTAACCCTCACCCCACAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg           <  2:1055312/150‑1 (MQ=37)
                                                                                                                                         ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc          >  1:1103176/1‑150 (MQ=255)
                                                                                                                                          cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGTGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg         <  2:44417/150‑1 (MQ=255)
                                                                                                                                          cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg         >  2:1161476/1‑150 (MQ=255)
                                                                                                                                              cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt     >  1:923127/1‑150 (MQ=255)
                                                                                                                                              cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt     >  1:1004519/1‑150 (MQ=255)
                                                                                                                                              cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt     >  1:1249374/1‑150 (MQ=255)
                                                                                                                                               ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCttt    <  1:698250/150‑1 (MQ=255)
                                                                                                                                                cTTTCCGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCGCTCCCCAAAGGGGCGAGGGGACTGGCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCCCCGAAGCGCCTTTa   <  1:999524/150‑1 (MQ=255)
                                                                                                                                                 tttCAGGGAGCGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCGGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAg  >  1:275561/1‑150 (MQ=255)
                                                                                                                                                 tttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCttttt  <  1:1195379/150‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAG  >  NC_000913/3045853‑3046147

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