Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,068 T→A 16.9% intergenic (‑167/+100) 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,0680TA16.9% 744.4 / 40.0 319intergenic (‑167/+100)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (66/185);  new base A (56/3);  total (128/191)
Fisher's exact test for biased strand distribution p-value = 1.47e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

ATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045919‑3046217
                                                                                                                                                     |                                                                                                                                                     
     ttCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAACCTCTccc                                                                                                                                                  >  2:2534533/1‑150 (MQ=255)
      tCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCACTcccc                                                                                                                                                 >  2:1222083/1‑150 (MQ=255)
         ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                              <  1:3151233/150‑1 (MQ=255)
          tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                             <  1:5926233/150‑1 (MQ=255)
              gTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGc                                                                                                                                         >  1:4426402/1‑150 (MQ=255)
                aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                                       <  2:2907239/150‑1 (MQ=255)
                 tCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                      <  1:1083311/150‑1 (MQ=255)
                   tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCCGGCGTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGGCTGGGCCCCTGCCCCTAAGCCTCTCCCCCGAGGGGGGCggg                                                                                                                                    >  2:5251883/1‑150 (MQ=255)
                   tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTATACCTCTCCCCAAAGGGGCGAggg                                                                                                                                    >  2:2407774/1‑150 (MQ=255)
                   cttcggcgcctttttagtcagatgacaaagtacaaaagtgctcagacagtCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                    >  1:1761457‑M1/51‑150 (MQ=255)
                    ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAACGGGCGCgggg                                                                                                                                   >  1:5276333/1‑150 (MQ=255)
                     tGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                                                  <  1:2143350/150‑1 (MQ=255)
                      gagaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgggtta                                                                                                                                 <  1:4506689/150‑4 (MQ=37)
                       agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAAACCTCACCCCAAAGGGGCGAGGGGACt                                                                                                                                >  1:3260845/1‑150 (MQ=255)
                       agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTCACCCTAAACCACTCCCCCAAGGGGCGAGGGGCCt                                                                                                                                >  1:217008/1‑150 (MQ=255)
                       agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACCCCCAAACCCACACCCCAAACGGGCGAggggccg                                                                                                                                >  1:2617081/1‑147 (MQ=255)
                          ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGCGGGGCCTGTc                                                                                                                             >  1:2883847/1‑150 (MQ=255)
                               aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTCACCCAAAACCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                                                                                        >  2:5715637/1‑150 (MQ=255)
                                  aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACACCCCAAAGGGGCGAGGGGACTGTCTGAGCAca                                                                                                                     >  1:1838410/1‑149 (MQ=255)
                                    aCTGGCTCAATCTTCAATGTTGGCACGCTGTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACatt                                                                                                                   <  2:3354121/150‑1 (MQ=255)
                                    aCTGGCTCAATCTTCAATGTTGGCACGCTATTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACttt                                                                                                                   <  2:5345451/150‑1 (MQ=255)
                                      tGGCTCACTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTg                                                                                                                 <  2:4657608/150‑1 (MQ=255)
                                      tGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTg                                                                                                                 <  1:5513821/150‑1 (MQ=255)
                                      tGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTg                                                                                                                 >  1:3385262/1‑150 (MQ=255)
                                          tCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACt                                                                                                             >  2:5902407/1‑150 (MQ=255)
                                            aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                                           <  2:5139559/150‑1 (MQ=255)
                                            aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                                           >  1:754911/1‑150 (MQ=255)
                                             aTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTg                                                                                                          <  2:862016/150‑1 (MQ=255)
                                              tCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAACTCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGt                                                                                                         <  2:2220409/150‑1 (MQ=255)
                                               cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACAACCCAAAGGGGCGGGGGGGCTGTCTGAGCACGTTTTGACATTGTc                                                                                                        >  2:2269962/1‑150 (MQ=255)
                                                 tCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATc                                                                                                     <  1:2336891/150‑1 (MQ=255)
                                                  cAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCGAACCCGCAACCCCAAGGGGCGAGGGGACTGTCGGAGCCCGTTTGTACTTTGTCATc                                                                                                     >  2:2146006/1‑150 (MQ=255)
                                                     tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                  <  1:3639293/150‑1 (MQ=255)
                                                        tGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTa                                                                                               <  2:1856281/150‑1 (MQ=255)
                                                        tGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTa                                                                                               <  2:4585193/150‑1 (MQ=255)
                                                           cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaa                                                                                            <  2:3073538/150‑1 (MQ=255)
                                                             cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg                                                                                          <  1:5402606/150‑1 (MQ=255)
                                                                  ccggCCCTCTCCCTGAAAGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGcc                                                                                     <  2:5940852‑M1/146‑1 (MQ=255)
                                                                      CCCTCGCCCTTTCAGGGAGAGGGCCGGGTTGAGGGTAAATATTCGCGCCAGTGCCTGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg                                                                                 >  1:408401‑M1/1‑150 (MQ=255)
                                                                        CTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg                                                                               >  2:2829167‑M1/1‑150 (MQ=255)
                                                                         TCGCCCTGTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                                              >  2:1564893‑M1/1‑150 (MQ=255)
                                                                           GCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTACCCAAAGGGGCGGGGGGACTGTCTGAGCACTTTTTTACTTTTTCATCTTAATAAAAAGGCGCCGAAGCGCCt                                                                            <  2:5106413‑M1/150‑1 (MQ=255)
                                                                             ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCCAACCCCCAACCCCAACGGGCGCGGGGGCGGGCTGAGCACGTTTTTACATTGTGATATGACTAAAAATTCGCCGACACGccggt                                                                          >  1:4066613/1‑147 (MQ=255)
                                                                                ttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                                       <  2:1421148/150‑1 (MQ=255)
                                                                                gtCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                                       <  2:5757589/149‑1 (MQ=255)
                                                                                   aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                                    <  1:5710971/150‑1 (MQ=255)
                                                                                   aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                                    <  1:4924527/150‑1 (MQ=255)
                                                                                     ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCGCACCCAAACCCACAACCCCAAGGGGCGGGGGGGCTGTGGGAGCACGTTTTTAATTTTTGATAAGACTAAAAATGCGCCGAAGCGCCGTTAGAAAATa                                                                  >  1:1425668/1‑150 (MQ=255)
                                                                                     ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                  <  1:2524795/150‑1 (MQ=255)
                                                                                         agGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCg                                                              <  1:4610201/150‑1 (MQ=255)
                                                                                          gggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGa                                                             <  2:3051192/150‑1 (MQ=255)
                                                                                          gggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGa                                                             <  2:1462291/150‑1 (MQ=255)
                                                                                           ggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                                            >  2:3504761/1‑150 (MQ=255)
                                                                                            gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAAt                                                           <  1:5650503/150‑1 (MQ=255)
                                                                                              cGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGAGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCa                                                         <  1:5039945/150‑1 (MQ=255)
                                                                                                gggTGAGGGTAAAGATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                                       <  1:1029770/150‑1 (MQ=255)
                                                                                                    gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAAt                                                   <  2:797811/150‑1 (MQ=255)
                                                                                                     aGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAAtt                                                  <  2:1338179/150‑1 (MQ=255)
                                                                                                      gtgTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTa                                                 <  1:565802/148‑1 (MQ=255)
                                                                                                      gggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTa                                                 <  1:4357630/150‑1 (MQ=255)
                                                                                                       ggTTAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAc                                                <  1:2049023/150‑1 (MQ=255)
                                                                                                       ggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCGCTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAc                                                <  1:1806331/150‑1 (MQ=255)
                                                                                                       ggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAc                                                <  1:3904781/150‑1 (MQ=255)
                                                                                                        gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                               <  2:5416042/150‑1 (MQ=255)
                                                                                                         tAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTg                                              <  2:1146369/150‑1 (MQ=255)
                                                                                                         tAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTg                                              <  1:2357883/150‑1 (MQ=255)
                                                                                                          aaaTATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTgg                                             >  2:1959524/1‑150 (MQ=255)
                                                                                                              atTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTAtt                                         <  2:4817715/150‑1 (MQ=255)
                                                                                                               ttCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTc                                        <  1:4452832/150‑1 (MQ=255)
                                                                                                               ttCGCGCCAGTGCCGGCCTGGTGCCCACACCCAAAACCACTCCCCAAAGGGGCGGGGGGAAAGTCGGGGCACGTTTGTACGTTGTCATATTACTAAAAATGCGCCGAAGCGCCGTTAGAAAAAATTAGAATAAGTTAAATAATGGTATTc                                        >  2:4034754/1‑150 (MQ=255)
                                                                                                               ttCGCGCCAGTGCCGGCCTGGTACCCACACCCAAAACCACAACCCAAAGGGGGGGGGGGACTGTGTTCTCACGTTTTTACATTGTAATAATACAAAAAAAGCGCCCAATCGCCGTTAGAAAAAATTGGAATCAGTGTATTACTTGTATTc                                        >  2:2575784/1‑150 (MQ=255)
                                                                                                                 cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTTTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                                      <  2:2095497/150‑1 (MQ=255)
                                                                                                                 cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                                      <  1:1508867/150‑1 (MQ=255)
                                                                                                                 cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                                      <  1:1297738/150‑1 (MQ=255)
                                                                                                                 cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                                      <  2:4426402/150‑1 (MQ=255)
                                                                                                                       aGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCg                                <  2:2602030/150‑1 (MQ=255)
                                                                                                                         tGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGt                              <  1:1326157/150‑1 (MQ=255)
                                                                                                                          gCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTa                             <  1:5362322/150‑1 (MQ=255)
                                                                                                                          gCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTa                             <  2:647527/150‑1 (MQ=255)
                                                                                                                          gCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTa                             <  1:5887134/150‑1 (MQ=255)
                                                                                                                            cGGCCTGTTACCCACACCCAAAACCACTCCCCAAAGGGGCGAGGGGACTTTCTGATCACGGTTGTACATTGTCATAAGACAAAAAAAGCGCCCAAGCGCCATTAGGAAATATTCGAAACAGTGAATTACAGGTTTTCTCTAAACGGTaag                           >  2:5686595/1‑148 (MQ=255)
                                                                                                                              gCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAATGGTATTAGCTAATCGGTACGCa                         >  2:421571/1‑150 (MQ=255)
                                                                                                                               ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg                        <  1:3071358/150‑1 (MQ=255)
                                                                                                                               ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg                        <  1:3036138/150‑1 (MQ=255)
                                                                                                                                 tggtCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa                      <  1:4582241/147‑1 (MQ=255)
                                                                                                                                 tGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa                      <  1:5748341/150‑1 (MQ=255)
                                                                                                                                 tGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa                      >  2:1926963/1‑150 (MQ=255)
                                                                                                                                    tCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCa                   <  2:4475830/150‑1 (MQ=255)
                                                                                                                                      cccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGa                 <  2:3771608/150‑1 (MQ=255)
                                                                                                                                        cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGACGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAAGCGGTACGCAGGAGCAGAAc               >  2:1504177/1‑150 (MQ=255)
                                                                                                                                        cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc               <  2:3669348/150‑1 (MQ=255)
                                                                                                                                         tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACa              <  1:892860/150‑1 (MQ=255)
                                                                                                                                            cccTAACCCGCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGt           <  2:4847921/150‑1 (MQ=255)
                                                                                                                                              cTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGttc         >  2:3528478/1‑149 (MQ=255)
                                                                                                                                              cTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTa         >  1:3763757/1‑150 (MQ=255)
                                                                                                                                               tAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTAc        >  2:1066554/1‑150 (MQ=255)
                                                                                                                                               tAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTAc        >  2:1575571/1‑150 (MQ=255)
                                                                                                                                                    ctctCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCg   <  2:2133843/150‑1 (MQ=255)
                                                                                                                                                     tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  >  1:1345551/1‑150 (MQ=255)
                                                                                                                                                     tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  <  1:1198271/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045919‑3046217

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: