Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,335,050 T→G intergenic (+221/‑74) cysB → / → ymiA N‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,0500TG97.6% 314.7 / ‑1.6 84intergenic (+221/‑74)cysB/ymiAN‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (42/40);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

TATACCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGA  >  NC_000913/1334902‑1335197
                                                                                                                                                    |                                                                                                                                                   
tataCCGTGCACAATATTTTTATCCCTCATATTTATAGGGCAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCg                                                                                                                                                   >  1:768716/1‑151 (MQ=255)
 ataCCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATgcgc                                                                                                                                                  >  1:673970/1‑151 (MQ=255)
 ataCCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATgcgc                                                                                                                                                  >  2:876631/1‑151 (MQ=255)
  taCCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCt                                                                                                                                                 >  1:709157/1‑151 (MQ=255)
     cGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTgg                                                                                                                                               >  1:445081/1‑150 (MQ=255)
        gCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGAttt                                                                                                                                           <  2:935155/151‑1 (MQ=255)
         cacaATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGAtttt                                                                                                                                          >  2:49845/1‑151 (MQ=255)
         cacaATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGAtttt                                                                                                                                          <  2:1063961/151‑1 (MQ=255)
                 tttaaacccTCATATTAATAGGGTAAATTCCCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTa                                                                                                                                    <  1:1239273/143‑1 (MQ=255)
                  tttATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTaa                                                                                                                                   >  1:740922/1‑149 (MQ=255)
                    tATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGtt                                                                                                                               >  2:923290/1‑151 (MQ=255)
                       cccTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCt                                                                                                                            >  1:564254/1‑151 (MQ=255)
                        ccTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCtt                                                                                                                           <  2:1205556/151‑1 (MQ=255)
                        ccTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCtt                                                                                                                           >  1:1113049/1‑151 (MQ=255)
                          tCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCttt                                                                                                                          >  1:682209/1‑150 (MQ=255)
                           cATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAg                                                                                                                        <  2:891391/151‑1 (MQ=255)
                            atatTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAgg                                                                                                                       <  2:1202295/151‑1 (MQ=255)
                              atTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGAt                                                                                                                     <  1:929382/151‑1 (MQ=255)
                              atTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGAt                                                                                                                     <  1:222759/151‑1 (MQ=255)
                               tttATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGAtt                                                                                                                    >  2:622896/1‑151 (MQ=255)
                               tttATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGAt                                                                                                                     <  2:662842/150‑1 (MQ=255)
                               tttATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGAt                                                                                                                     <  2:529278/150‑1 (MQ=255)
                                ttATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGAttt                                                                                                                   >  2:1251558/1‑151 (MQ=255)
                                 tataGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGAttt                                                                                                                   <  1:971351/150‑1 (MQ=255)
                                 tataGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTg                                                                                                                  >  2:167264/1‑151 (MQ=255)
                                  ataGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTAGTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTg                                                                                                                  <  2:63546/150‑1 (MQ=255)
                                        tAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATc                                                                                                           >  1:1118439/1‑151 (MQ=255)
                                         aaaTTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCg                                                                                                          <  1:644649/151‑1 (MQ=255)
                                          aaTTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGt                                                                                                         >  1:1047349/1‑151 (MQ=255)
                                           aTTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGt                                                                                                         >  1:149109/1‑150 (MQ=255)
                                                   tAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGattaa                                                                                                 <  2:1243122/150‑1 (MQ=255)
                                                    aaaGTAGTGTCTAATTCCCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGattaat                                                                                                <  2:761721/150‑1 (MQ=255)
                                                    aaaGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGattaa                                                                                                 <  1:902269/149‑1 (MQ=255)
                                                    aaaGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGattaat                                                                                                <  1:412899/150‑1 (MQ=255)
                                                     aaGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGattaatta                                                                                              <  2:76811/151‑1 (MQ=255)
                                                      agtagtGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTac                                                                                             >  1:745256/1‑151 (MQ=255)
                                                       gtagtGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTaca                                                                                            <  2:155782/151‑1 (MQ=255)
                                                             tCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATc                                                                                      >  2:550974/1‑151 (MQ=255)
                                                               tAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGt                                                                                    >  1:532614/1‑151 (MQ=255)
                                                                     cccGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAAc                                                                               <  2:1113118/150‑1 (MQ=255)
                                                                      ccGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCCTTCAACg                                                                              <  1:560091/150‑1 (MQ=255)
                                                                           tCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGc                                                                        >  2:159464/1‑151 (MQ=255)
                                                                               aaaCTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGAt                                                                    >  2:979371/1‑151 (MQ=255)
                                                                                aaCTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATa                                                                   <  2:1184144/151‑1 (MQ=255)
                                                                                 aCTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATaa                                                                  <  2:1085313/151‑1 (MQ=255)
                                                                                 aCTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATaa                                                                  <  2:849219/151‑1 (MQ=255)
                                                                                   tCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGt                                                                >  2:101189/1‑151 (MQ=255)
                                                                                    ccAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAg                                                                 <  1:209385/149‑1 (MQ=255)
                                                                                     cAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTg                                                               >  2:1255789/1‑150 (MQ=255)
                                                                                      aGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGAt                                                             >  1:617971/1‑151 (MQ=255)
                                                                                        ttttCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATgg                                                           >  2:1240942/1‑151 (MQ=255)
                                                                                        ttttCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATgg                                                           >  1:312014/1‑151 (MQ=255)
                                                                                        ttttCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATg                                                            <  2:745525/150‑1 (MQ=255)
                                                                                          ttCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCg                                                         >  1:509785/1‑151 (MQ=255)
                                                                                             aaaTATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGAtt                                                      >  2:248329/1‑151 (MQ=255)
                                                                                              aaTATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATta                                                     <  2:940587/151‑1 (MQ=255)
                                                                                              aaTATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATta                                                     <  2:321472/151‑1 (MQ=255)
                                                                                               atatCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATtat                                                    <  2:661378/151‑1 (MQ=255)
                                                                                                tatCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATtat                                                    >  2:365660/1‑150 (MQ=255)
                                                                                                 atCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATtata                                                   <  1:1075973/150‑1 (MQ=255)
                                                                                                  tCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTAGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATg                                                 >  2:86191/1‑151 (MQ=255)
                                                                                                   cTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATg                                                 >  1:808829/1‑150 (MQ=255)
                                                                                                   cTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATg                                                 >  2:1070093/1‑150 (MQ=255)
                                                                                                   cTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATg                                                 >  2:448258/1‑150 (MQ=255)
                                                                                                    tATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATg                                                 >  2:1213990/1‑149 (MQ=255)
                                                                                                     aTTTCAGGTCAACAGATTGAATTATTCACTGGAGACGATTCGTCAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCCTCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAgg                                              >  1:1036076/1‑151 (MQ=255)
                                                                                                     aTTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGCCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATg                                                 <  2:60360/148‑1 (MQ=255)
                                                                                                         cAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTa                                           <  1:90104/150‑1 (MQ=255)
                                                                                                           ggTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCa                                        <  1:479794/151‑1 (MQ=255)
                                                                                                              cAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAAt                                      <  2:1239768/150‑1 (MQ=255)
                                                                                                              cAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATg                                     >  1:755279/1‑151 (MQ=255)
                                                                                                                 aaGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCt                                  >  2:1091763/1‑151 (MQ=255)
                                                                                                                 aaGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCt                                  >  1:176279/1‑151 (MQ=255)
                                                                                                                  aGATTGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCtt                                 <  2:243186/151‑1 (MQ=255)
                                                                                                                      tGAATTATTCACTGGAGTCGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCGTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTg                              <  2:1133007/150‑1 (MQ=255)
                                                                                                                      tGAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTgg                             >  1:193009/1‑151 (MQ=255)
                                                                                                                       gAATTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTgg                             >  1:412063/1‑150 (MQ=255)
                                                                                                                        aaTTATTCACTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGa                            >  2:540032/1‑150 (MQ=255)
                                                                                                                                aCTGGAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGaa                   <  2:20219/151‑1 (MQ=255)
                                                                                                                                   ggAGACGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAAc                  <  2:445700/149‑1 (MQ=255)
                                                                                                                                    gagaCGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAAcccc               <  2:157985/151‑1 (MQ=255)
                                                                                                                                       aCGATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCccgcc            <  2:1211025/151‑1 (MQ=255)
                                                                                                                                         gATTCGTAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCcgccgc          >  2:401391/1‑151 (MQ=255)
                                                                                                                                               tAAATGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCt    >  1:1168368/1‑151 (MQ=255)
                                                                                                                                                aaaTGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCTTTCTGGAAATCAGGAACCCCGCCGCGATCCt    <  1:49844/150‑1 (MQ=255)
                                                                                                                                                 aaTGCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGACCCCCGCCGCGATCCTGa  <  1:101187/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATACCGTGCACAATATTTTTATCCCTCATATTTATAGGGTAAATTACCTGTAAAGTAGTGTCTAATTACCCGCCTCGCAAACTCCAGTTTTCAAATATCTATTTCAGGTCAAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGA  >  NC_000913/1334902‑1335197

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