Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,006 A→C 7.3% intergenic (‑105/+162) 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,0060AC 7.3% 301.8 / 13.7 110intergenic (‑105/+162)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base A (39/62);  new base C (4/4);  total (43/67)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

Only 100 of 117 total aligned reads displayed.

CCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAG  >  NC_000913/3045857‑3046154
                                                                                                                                                     |                                                                                                                                                    
ccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGga                                                                                                                                                      <  2:1074141/150‑1 (MQ=255)
ccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTAAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGga                                                                                                                                                      <  2:1079732/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)
     cgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAggg                                                                                                                                                 >  2:679534/1‑150 (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                                                                                                                                           <  1:579004/150‑1 (MQ=255)
           cgcgACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCgggg                                                                                                                                           <  2:1617860/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)
                           aacaAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcg                                                                                                                           <  2:1607280/149‑1 (MQ=255)
                             caAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcgcc                                                                                                                         >  1:156331/1‑149 (MQ=255)
                             caAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACATATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcgcg                                                                                                                         <  2:372817/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)
                                     ataGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGcc                                                                                                                 >  2:1456596/1‑150 (MQ=255)
                                      taGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGTTAGGGTAAATATTCGCGCCAGTGCCg                                                                                                                <  1:387496/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)
                                               gTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTcc                                                                                                       >  1:287711/1‑150 (MQ=255)
                                                  tttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAGCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                                                    <  1:1169911/150‑1 (MQ=255)
                                                      cGAGAAAAATTATTTCGTTGACGTATCTTTGAGTTTGCTAACAATCTGGCTCAATCTTCTATGTTGGCACGCTGTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                                <  2:720157/150‑1 (MQ=255)
                                                      cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                                >  1:318772/1‑150 (MQ=255)
                                                      cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                                >  2:561321/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)
                                                       gagaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAccc                                                                                               <  1:1513701/150‑1 (MQ=255)
                                                         gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCCCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAATTATTCGCGCCCGTGCCGGCCTGTTCCCCTCACCCTa                                                                                             >  2:158256/1‑150 (MQ=255)
                                                         gaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTa                                                                                             >  1:23355/1‑150 (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)
                                                                      gTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaa                                                                                >  1:1596217/1‑150 (MQ=255)
                                                                       ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                               <  1:1564530/150‑1 (MQ=255)
                                                                        tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                              <  2:210625/150‑1 (MQ=255)
                                                                            gTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACACTAACCCTCTACCCAAAGGGGc                                                                          >  1:842958/1‑150 (MQ=255)
                                                                             tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                                         <  2:791479/150‑1 (MQ=255)
                                                                             tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                                         <  2:44453/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)
                                                                                         tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCTGTCt                                                             >  2:1448137/1‑150 (MQ=255)
                                                                                          gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTg                                                            >  1:1183645/1‑150 (MQ=255)
                                                                                          gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGGGAGGGTAAATATTCGCGCCAGTGCCGGCATGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTg                                                            <  2:858942/150‑1 (MQ=255)
                                                                                            tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTTTCTGAg                                                          >  2:851496/1‑150 (MQ=255)
                                                                                             aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                         >  1:633865/1‑150 (MQ=255)
                                                                                                aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGGCTGTCTGAGCACt                                                      >  1:218198/1‑150 (MQ=255)
                                                                                                 accTGGCGCAATCTTCAATGTTGGCACGCTGTTCCCCCTCGCCCTTTCCGGGTGAGGGCCGGGGTCTGGGTAAATATTCGCGCCATTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACtt                                                     <  1:1272724/148‑1 (MQ=25)
                                                                                                     ggttCAATCTTCAAGGTGGGCAGGCTTTTCCCCCTGGCCCTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                                <  2:198395/147‑1 (MQ=25)
                                                                                                           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)
                                                                                                                 aaTGTTGGCACGCTTTTCCCCCTCACCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCCGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                                     <  1:1375787/150‑1 (MQ=255)
                                                                                                                  aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGGGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                                    <  1:764522/150‑1 (MQ=255)
                                                                                                                  aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                                    <  2:588361/150‑1 (MQ=255)
                                                                                                                  aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                                    <  1:683557/150‑1 (MQ=255)
                                                                                                                  aTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTg                                    <  2:719102/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)
                                                                                                                           cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg                           >  2:473256/1‑150 (MQ=255)
                                                                                                                               cccggcccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGCGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                       >  2:1389466/6‑150 (MQ=25)
                                                                                                                                ccggcccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGcc                      >  2:1167512/5‑150 (MQ=25)
                                                                                                                                    cccTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg                  >  1:526649/1‑150 (MQ=25)
                                                                                                                                     ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc                 >  1:1103176/1‑150 (MQ=255)
                                                                                                                                      cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGTGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg                <  2:44417/150‑1 (MQ=255)
                                                                                                                                      cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg                >  2:1161476/1‑150 (MQ=255)
                                                                                                                                          cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt            >  1:1004519/1‑150 (MQ=255)
                                                                                                                                          cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt            >  1:923127/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)
                                                                                                                                                  gggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCTGTCTGCGCACTATTGTAGTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAAt    >  2:701002/1‑150 (MQ=255)
                                                                                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAtata   >  1:263454/1‑150 (MQ=255)
                                                                                                                                                    gagaGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAg  <  1:1425792/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAG  >  NC_000913/3045857‑3046154

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: