Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,627,813 G→T 13.8% E167* (GAA→TAA)  nadR → trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,8130GT13.8% 282.7 / 2.9 115E167* (GAA→TAA) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (49/46);  new base T (2/17);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 8.62e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

AGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTC  >  NC_000913/4627664‑4627962
                                                                                                                                                     |                                                                                                                                                     
agcCGACCGTGCCGGATCGTCTGCGTTGTTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTg                                                                                                                                                       <  1:1363605/150‑1 (MQ=255)
agcCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACTGCATCAAAAAGTTTATGGCTg                                                                                                                                                       >  2:463516/1‑150 (MQ=255)
          gCCGGATCGTCTGCGTTGGATATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGAt                                                                                                                                             <  2:114652/150‑1 (MQ=255)
           ccGTATCGTCTGCGTTGTTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATc                                                                                                                                            <  2:1745822/150‑1 (MQ=255)
            cGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATcc                                                                                                                                           <  1:356654/150‑1 (MQ=255)
             ggATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCa                                                                                                                                          >  1:1347519/1‑150 (MQ=255)
               aTCGTCTGCGTTGGTTATTTCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACTAAGAGGGCATGGAGCCGTTTCCGCACGGCTGGGATGTGTGGAGCAACGGCAGCAAAAAGTTTATGGCTGAAAAAGGGATCCAGc                                                                                                                                        <  1:155383/150‑1 (MQ=255)
                 cGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCg                                                                                                                                      <  1:1418879/150‑1 (MQ=255)
                  gTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCTCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCgg                                                                                                                                     <  2:1640279/150‑1 (MQ=255)
                  gTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCgg                                                                                                                                     >  2:1309378/1‑150 (MQ=255)
                     tGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCTGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatc                                                                                                                                  <  2:1625739/150‑1 (MQ=255)
                     tGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAGAGGGATCCAGCCGgatc                                                                                                                                  >  2:189467/1‑150 (MQ=255)
                     tGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatc                                                                                                                                  <  1:975683/150‑1 (MQ=255)
                      gCGTTGGTTATTGCAAACTTTTATATATCAGAAAAATATTCGCATTCATGCTTTCATCGAAGAGGGCCTGGAGCCGTATCCGCACGGCTGGGAGGGGGGGAGCAACGGCATCAAAAAGTTAATGGCTGAAAAAGTGATCCAGCCGgatct                                                                                                                                 <  2:73816/150‑1 (MQ=255)
                      gCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATACGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGATGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatct                                                                                                                                 >  1:914735/1‑150 (MQ=255)
                          tGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCTCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                             <  2:1367886/150‑1 (MQ=255)
                          tGGTTATTGCAAACATATAAATATCAGAAAAAAATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                                             <  2:1373186/150‑1 (MQ=255)
                           ggTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                                                                                                            <  1:824491/150‑1 (MQ=255)
                             ttATTGCAAACTTTTAAATATCAGAAAAATATTCCCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTTAAAAAGGGATCCAGCCGGATCTGATCTac                                                                                                                          <  2:316492/150‑1 (MQ=255)
                               aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacac                                                                                                                        <  1:1049445/150‑1 (MQ=255)
                               aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacac                                                                                                                        >  1:233881/1‑150 (MQ=255)
                                    aaaCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCgg                                                                                                                   >  2:1479973/1‑150 (MQ=255)
                                      aCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaa                                                                                                                 >  1:847149/1‑150 (MQ=255)
                                         tttAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaagaa                                                                                                              >  2:867548/1‑150 (MQ=255)
                                         tttAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaagaa                                                                                                              >  2:199475/1‑150 (MQ=255)
                                            aaaTATCAGAAAACTATTCGCATTCATGCTTTCAACGACGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTTTAAAGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGcc                                                                                                           <  2:1762778/150‑1 (MQ=255)
                                             aaTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCg                                                                                                          >  1:1565817/1‑150 (MQ=255)
                                              atatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGa                                                                                                         <  1:249669/150‑1 (MQ=255)
                                               tatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGAt                                                                                                        >  2:1037437/1‑150 (MQ=255)
                                               tatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAAGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAATTTTATTTCTAAAAAAGTCATCCATCCGGATCTGATCTACACCTCGGAAGAAGCCGAt                                                                                                        >  2:972501/1‑150 (MQ=255)
                                               aatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGACGGCATGGACCAGAATCCGCACCGCTGCGAAGTGTGGAGCAACGCCATCAAAAAGTTTATCGCTTAAAAAGCGAAACAGCCGGATCTGAACAACACCTCCGAACAAGCCGaa                                                                                                        >  2:1092574/2‑149 (MQ=255)
                                                atCAGAAAAATATTCGCAGTCCTGCTTTCACCGCAGAGGCCATGGACCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTGTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATg                                                                                                       <  2:1619560/150‑1 (MQ=255)
                                                    gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCGTAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                                   >  2:705449/1‑150 (MQ=255)
                                                    gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAGGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGACGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                                   >  1:1189306/1‑150 (MQ=255)
                                                            ttCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGtata                                                                                           >  1:1327301/1‑150 (MQ=255)
                                                              cGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                                                         <  1:675368/150‑1 (MQ=255)
                                                               gCATTCCTGCTTTCAAGGAACAGGGCATGGAGCCGGAGCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAATTTTATTTTTTTAAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATgg                                                                                        <  2:494302/150‑1 (MQ=255)
                                                               gCATTCATGCTTTCAACGAAGAGGGCCTGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTTTAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATgg                                                                                        <  2:1851444/150‑1 (MQ=255)
                                                               gCATTCATGCTTTAAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGGGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATgg                                                                                        <  1:99755/150‑1 (MQ=255)
                                                               gCATTAATGCTTGCAACGAAGAGGTCATGGAGACGTATCCGCACGGCTGGGATGTGTGGATCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTAGGAAGAAGCCGATGCGCCACAGTATATgg                                                                                        <  1:298402/150‑1 (MQ=255)
                                                                 aTTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGGGTGCAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGaa                                                                                      <  1:203683/150‑1 (MQ=255)
                                                                    cATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCCCAGTATATGGAACAt                                                                                   >  2:207054/1‑150 (MQ=255)
                                                                    cATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACAt                                                                                   >  1:1835323/1‑150 (MQ=255)
                                                                      tGCTTTCAACCAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTTAAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCt                                                                                 <  2:1307382/150‑1 (MQ=255)
                                                                        cTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTgg                                                                               >  2:1534618/1‑150 (MQ=255)
                                                                         attCAACGCAGAGCGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTggg                                                                              <  2:1489702/149‑1 (MQ=255)
                                                                            cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAt                                                                           <  1:1449843/150‑1 (MQ=255)
                                                                            cAACGAAGAGGGCATGGAGCCCTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCAACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAt                                                                           >  1:1005281/1‑150 (MQ=255)
                                                                             aaCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                                                          >  2:875363/1‑150 (MQ=255)
                                                                              cccaaGAGGGCAGGGAGCCGTATCCGCACGGCTGGTATGTGTGCAGCAACGGCATCAAAAAGTTTATTGTTTATAATGATATGGAGACGTAGGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                                                         <  2:288876/147‑1 (MQ=255)
                                                                              ccGAAGCGGGCATGGAGCCGTATCCGCACGGCGGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAACAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCg                                                                         <  1:355759/149‑1 (MQ=255)
                                                                                gAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCgag                                                                       <  1:825174/150‑1 (MQ=255)
                                                                                   gagGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACATCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACg                                                                    >  1:536871/1‑150 (MQ=255)
                                                                                      ggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTg                                                                 >  2:1348485/1‑150 (MQ=255)
                                                                                       gCATGGAGCCGTATCCGGACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTAATGGCTGATAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGc                                                                <  1:936563/150‑1 (MQ=255)
                                                                                       gCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGc                                                                >  2:1826515/1‑150 (MQ=255)
                                                                                        cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCACCGGCAGCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGCAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCt                                                               <  1:260483/150‑1 (MQ=255)
                                                                                        cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGTGGATCGAGACGGTGCt                                                               >  2:1412702/1‑150 (MQ=255)
                                                                                           ggAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGt                                                            <  1:867548/150‑1 (MQ=255)
                                                                                           ggAGCCGTATCCGCACGGCGGGGATGTGTGGACCCACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCAACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGt                                                            <  1:1643890/150‑1 (MQ=255)
                                                                                            gAGCCTTATCCGCACGGCTGGGATGTGTGGAGCAACGTCATCAAAAAGTTTATTGTTTTAAATGTGATGCAGACGTATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTc                                                           <  2:89505/150‑1 (MQ=255)
                                                                                              tccGTATCCGCACGGCTGGGATGTGTGTAGCAAGGTCATAAAAAAGTTTTTTTTTTATAATGTTATGGAGACGTATGTGATCTACAGCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGa                                                         <  2:1574789/149‑1 (MQ=255)
                                                                                                cGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGACCATCTGGGGATCGAGACGGTGCTGGTCGATc                                                       >  2:709449/1‑150 (MQ=255)
                                                                                                 gTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTATTAAAAAGTTATCCAGCCGTATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATcc                                                      <  2:1565306/150‑1 (MQ=255)
                                                                                                    tCCGCACGGCTGGGATGTGTGGAGCACCGGCATCAAAAAGTTTATGGCTGAAAAAGGGAACCAGCCGGAGCTGATCTACTCCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGaa                                                   <  1:68827/150‑1 (MQ=255)
                                                                                                      cGCACGGCTGGGATGTGTGGAGCATCGGCATCAAAAAGTTTTTTTATTTAAATGTGATGCAGCCGTATCTGATCTACTCCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAAc                                                 <  2:1356069/150‑1 (MQ=255)
                                                                                                      cGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAAc                                                 <  2:1356072/150‑1 (MQ=255)
                                                                                                         ccGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTa                                              <  1:830527/149‑1 (MQ=255)
                                                                                                         aCGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTGTTGAAAAAGTGATCCAGCCGGATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTa                                              <  2:64976/150‑1 (MQ=255)
                                                                                                            gCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCt                                           >  1:505081/1‑150 (MQ=255)
                                                                                                              tGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCttt                                         >  2:1050518/1‑150 (MQ=255)
                                                                                                                   tgtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATTCGCCACAGTATAAGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCATTATGAg                                    >  2:357804/1‑150 (MQ=255)
                                                                                                                   tgtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAg                                    >  2:1190527/1‑150 (MQ=255)
                                                                                                                      gtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTAt                                 <  2:1496002/150‑1 (MQ=255)
                                                                                                                      gtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTAt                                 >  1:1869989/1‑150 (MQ=255)
                                                                                                                       tgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATc                                >  1:623295/1‑150 (MQ=255)
                                                                                                                        ggAGCACCGGCATCAAAAAGATTATGGCTGAAAAAGCGATCCAGCCGGATCTGAACTACACCTCGAAAGAAGCCGATGCGCCACAGAAAATGGAACATCTGGGGAACGAGACGGTCCTGGTCGATCCGAAACGTACCTTTATGAGTATCa                               >  2:577342/1‑150 (MQ=255)
                                                                                                                        ggAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCa                               >  2:1871912/1‑150 (MQ=255)
                                                                                                                          aGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGc                             >  2:557851/1‑150 (MQ=255)
                                                                                                                           gCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCg                            <  1:1701871/150‑1 (MQ=255)
                                                                                                                           gCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCg                            <  1:1348485/150‑1 (MQ=255)
                                                                                                                             aaCGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTAAATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGt                          >  1:1080220/1‑150 (MQ=255)
                                                                                                                              actgCATAAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCAAAACGTACCTTTATGAGTATCAGCGGTg                         <  1:463516/147‑1 (MQ=255)
                                                                                                                                 tcATCAAAAAGTTTTTGTTTTATAATGTTATGCAGCCGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcgc                      <  2:1478741/149‑1 (MQ=255)
                                                                                                                                  cATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCa                     <  1:1033154/150‑1 (MQ=255)
                                                                                                                                     cAAAAAGTTTATGTGTGAAAAAGGGATCCAGCCGGATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAt                  <  2:1332152/150‑1 (MQ=255)
                                                                                                                                     cAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAt                  >  2:1264935/1‑150 (MQ=255)
                                                                                                                                      aaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                 >  2:369363/1‑150 (MQ=255)
                                                                                                                                       aaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTCCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATcc                >  2:1098622/1‑150 (MQ=255)
                                                                                                                                          aGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCcgcg             >  2:1295399/1‑150 (MQ=255)
                                                                                                                                          aGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCcgcg             <  1:1037437/150‑1 (MQ=255)
                                                                                                                                              tATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGaaaa         <  2:1559780/150‑1 (MQ=255)
                                                                                                                                               aTGGCTGCAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAAc        <  1:194598/150‑1 (MQ=255)
                                                                                                                                                tGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGCTCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAAcc       >  1:1842661/1‑150 (MQ=255)
                                                                                                                                                 ggCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAAccc      >  1:1285598/1‑150 (MQ=255)
                                                                                                                                                 ggCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAAccc      >  2:920748/1‑150 (MQ=255)
                                                                                                                                                  gCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCg     >  2:166393/1‑150 (MQ=255)
                                                                                                                                                   cTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAACCCCGt    >  2:1218959/1‑150 (MQ=255)
                                                                                                                                                    tGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGtt   >  2:509196/1‑150 (MQ=255)
                                                                                                                                                     gAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACAGCTGGGGATCGAGACGTTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGCTCCGCGAAAACCCGTTc  >  2:1331137/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTC  >  NC_000913/4627664‑4627962

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: