Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,627,809 G→T 11.0% M165I (ATG→ATT 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,8090GT11.0% 293.6 / 2.1 117M165I (ATG→ATTnadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (49/50);  new base T (2/15);  total (52/65)
Fisher's exact test for biased strand distribution p-value = 3.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

AGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCC  >  NC_000913/4627661‑4627958
                                                                                                                                                    |                                                                                                                                                     
agcagcCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATgg                                                                                                                                                      >  1:1332152/1‑150 (MQ=255)
 tcagcCGACCGTGCCGGATCGTCTGCGATGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTGc                                                                                                                                                     >  1:1356069/2‑150 (MQ=255)
 gcagcTGACCGTGCCGGATCGTCTGCCTTGGTTATTCCAAACTTATAAATAACAGAAAAATATTCCCATTCATGCTTTCAACCAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGTCATCAAAAAGATTATgga                                                                                                                                                     >  2:1644351/1‑149 (MQ=255)
 gcagcCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGc                                                                                                                                                     <  1:608263/150‑1 (MQ=255)
 gcagcCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGc                                                                                                                                                     >  1:1356072/1‑150 (MQ=255)
  cagcCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCt                                                                                                                                                    >  1:1559780/1‑150 (MQ=255)
  cagcCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCt                                                                                                                                                    >  1:1060070/1‑150 (MQ=255)
   agcCGACCGTGCCGGATCGTCTGCGTTGTTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTg                                                                                                                                                   <  1:1363605/150‑1 (MQ=255)
   agcCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACTGCATCAAAAAGTTTATGGCTg                                                                                                                                                   >  2:463516/1‑150 (MQ=255)
             gctggATCCTCTGCGTTGGTTATAGCAAACTTTAAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGAt                                                                                                                                         <  2:1537023/147‑1 (MQ=255)
             gCCGGATCGTCTGCGTTGGATATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGAt                                                                                                                                         <  2:114652/150‑1 (MQ=255)
              ccGTATCGTCTGCGTTGTTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATc                                                                                                                                        <  2:1745822/150‑1 (MQ=255)
              ccGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATc                                                                                                                                        <  2:1641245/150‑1 (MQ=255)
               cGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATcc                                                                                                                                       <  1:420448/150‑1 (MQ=255)
               cGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATcc                                                                                                                                       <  1:356654/150‑1 (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)
                        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)
                                ttATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAGGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTac                                                                                                                      <  1:1754100/150‑1 (MQ=255)
                                ttATTGCAAACTTTTAAATATCAGAAAAATATTCCCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTTAAAAAGGGATCCAGCCGGATCTGATCTac                                                                                                                      <  2:316492/150‑1 (MQ=255)
                                  aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacac                                                                                                                    >  1:233881/1‑150 (MQ=255)
                                  aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTacac                                                                                                                    <  1:1049445/150‑1 (MQ=255)
                                       aaaCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCgg                                                                                                               >  2:1479973/1‑150 (MQ=255)
                                         aCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaa                                                                                                             >  1:847149/1‑150 (MQ=255)
                                            tttAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaagaa                                                                                                          >  2:199475/1‑150 (MQ=255)
                                             ttAAATATGAGAAAACAATTCGCATTCAAGCTATCACCGAAGCCGCCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGTTATCCAGCCGTATGTGATCTACACCTCGGaagaag                                                                                                         <  2:185630/150‑1 (MQ=255)
                                                aaTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCg                                                                                                      >  1:1565817/1‑150 (MQ=255)
                                                 atatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGa                                                                                                     <  1:249669/150‑1 (MQ=255)
                                                  tatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGAt                                                                                                    >  1:1268104/1‑150 (MQ=255)
                                                  tatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGAt                                                                                                    >  2:1037437/1‑150 (MQ=255)
                                                  tatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAAGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAATTTTATTTCTAAAAAAGTCATCCATCCGGATCTGATCTACACCTCGGAAGAAGCCGAt                                                                                                    >  2:972501/1‑150 (MQ=255)
                                                  aatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGACGGCATGGACCAGAATCCGCACCGCTGCGAAGTGTGGAGCAACGCCATCAAAAAGTTTATCGCTTAAAAAGCGAAACAGCCGGATCTGAACAACACCTCCGAACAAGCCGaa                                                                                                    >  2:1092574/2‑149 (MQ=255)
                                                       gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCGTAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                               >  2:705449/1‑150 (MQ=255)
                                                       gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAGGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGACGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                               >  1:1189306/1‑150 (MQ=255)
                                                               ttCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGtata                                                                                       >  1:1327301/1‑150 (MQ=255)
                                                                tCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGACGCCGATGCCACACAGTatat                                                                                      >  2:338414/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)
                                                                    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)
                                                                             ttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATCTgggg                                                                         >  1:238072/1‑150 (MQ=255)
                                                                               cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAt                                                                       <  1:1449843/150‑1 (MQ=255)
                                                                               cAACGAAGAGGGCATGGAGCCCTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCAACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAt                                                                       >  1:1005281/1‑150 (MQ=255)
                                                                                aaCGAAGAGGGCCTGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTTAAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                                                      <  2:542306/150‑1 (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)
                                                                                         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)
                                                                                                                cTGGGATGTGTGTAGCAACGGCATCAAAAAGTTTATGTTTTTAAATGTTATGGAGACGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCtt                                      <  2:1060070/150‑1 (MQ=255)
                                                                                                                      tgtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAg                                >  2:1190527/1‑150 (MQ=255)
                                                                                                                         gtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTAt                             >  1:1869989/1‑150 (MQ=255)
                                                                                                                          tgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATc                            >  1:623295/1‑150 (MQ=255)
                                                                                                                           ggAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCa                           >  2:1871912/1‑150 (MQ=255)
                                                                                                                             aGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGc                         >  2:557851/1‑150 (MQ=255)
                                                                                                                              gCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCg                        <  1:1348485/150‑1 (MQ=255)
                                                                                                                                 actgCATAAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCAAAACGTACCTTTATGAGTATCAGCGGTg                     <  1:463516/147‑1 (MQ=255)
                                                                                                                                    tcATCAAAAAGTTTTTGTTTTATAATGTTATGCAGCCGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcgc                  <  2:1478741/149‑1 (MQ=255)
                                                                                                                                     cATCAAAAAGTTTTTTGTTTAAAAAGTTATCCAGCCGTATGTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCa                 <  2:1464474/150‑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)
                                                                                                                                                ttATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCAACAGTATATGGAACATCTGGGGATCGAGACTGTGCTGGTCGATACGAAACGAACCTTTAAGAGTATCAGCGGTGCGCAGATCCCCGaaa      >  2:767123/1‑150 (MQ=255)
                                                                                                                                                 tATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGaaca     >  1:774252/1‑148 (MQ=255)
                                                                                                                                                 tATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGaaaa     <  2:1559780/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)
                                                                                                                                                    |                                                                                                                                                     
AGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCC  >  NC_000913/4627661‑4627958

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: