Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,357,012 Δ1 bp 6.1% coding (246/681 nt) modB → molybdate ABC transporter permease ModB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,357,0120G. 6.1% 599.6 / 32.5 180coding (246/681 nt)modBmolybdate ABC transporter permease ModB
Reads supporting (aligned to +/- strand):  ref base G (65/104);  new base . (6/5);  total (71/109)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 180 total aligned reads displayed.

CTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAG  >  NC_002947/4356865‑4357157
                                                                                                                                                   |                                                                                                                                                  
 tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                                                                                                                                                 <  1:4853818/150‑1 (MQ=255)
 tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                                                                                                                                                 >  1:1767674/1‑150 (MQ=255)
 tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTAGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                                                                                                                                                 <  2:2365946/150‑1 (MQ=255)
   gtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGTCTTGGTCAGGCCACCCAGGCGCTAgg                                                                                                                                              >  1:1151543/1‑150 (MQ=255)
    tggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGc                                                                                                                                             <  2:4109041/150‑1 (MQ=255)
      gCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCt                                                                                                                                           <  1:149990/150‑1 (MQ=255)
           cgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtg                                                                                                                                      <  2:4157348/150‑3 (MQ=255)
                aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                 <  1:2699664/150‑1 (MQ=255)
                aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                 <  1:4930356/150‑1 (MQ=255)
                aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                 >  2:1544639/1‑150 (MQ=255)
                aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtggc                                                                                                                                    <  1:333488/147‑5 (MQ=255)
                aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtggc                                                                                                                                    >  2:333488/1‑143 (MQ=255)
                aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGGGGTGGCACTGCCGCGGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                 <  2:4194632/150‑1 (MQ=255)
                  gcgcTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggt                                                                                                                               >  1:2000335/1‑150 (MQ=255)
                     ctctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGtt                                                                                                                            <  2:4152837/150‑1 (MQ=255)
                      tctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTc                                                                                                                           <  2:3930037/150‑1 (MQ=255)
                      tctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtg                                                                                                                                      <  2:1689977/139‑3 (MQ=255)
                       ctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCa                                                                                                                          <  2:4558605/150‑1 (MQ=255)
                             tGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTca                                                                                                                    >  2:4665724/1‑150 (MQ=255)
                               cgcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTcaca                                                                                                                  <  1:3746935/150‑1 (MQ=255)
                                      cGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTgg                                                                                                           <  1:2374557/150‑1 (MQ=255)
                                        gtgggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGTTCGGTTTTTAGCTGTTGATCGCGCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCGTGGGCAGTGTGGTGTTCAGGTTCTGAGGCCTGGTg                                                                                                         >  2:860241/1‑150 (MQ=255)
                                           ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                      <  2:3041964/150‑1 (MQ=255)
                                           ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                      <  1:1534283/150‑1 (MQ=255)
                                            gtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTg                                                                                                     <  1:4744991/150‑1 (MQ=255)
                                               cGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                                  <  1:1862049/150‑1 (MQ=255)
                                               cGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                                  <  1:4112051/150‑1 (MQ=255)
                                                     tggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACgg                                                                                            >  1:999229/1‑150 (MQ=255)
                                                      ggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                           <  2:105855/150‑1 (MQ=255)
                                                      ggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                           <  1:4903789/150‑1 (MQ=255)
                                                      ggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                           <  1:3466285/150‑1 (MQ=255)
                                                      ggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCGGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                           <  2:257030/150‑1 (MQ=255)
                                                       gCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg                                                                                          >  2:356697/1‑150 (MQ=255)
                                                       gCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg                                                                                          <  1:1426579/150‑1 (MQ=255)
                                                              cGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCa                                                                                   <  1:2163469/150‑1 (MQ=255)
                                                                gtGGTACTGCCACCGACGGTGATCGGATTTTACCTGTTGATCGCCCTTGGCCCGCACGGATGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGTCAGTGTGGTGTTCTGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATg                                                                                 <  1:2335220/149‑1 (MQ=255)
                                                                cTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATg                                                                                 <  1:2825752/150‑1 (MQ=255)
                                                                       tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCg                                                                          >  2:1684269/1‑150 (MQ=255)
                                                                       tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCg                                                                          >  2:3608129/1‑150 (MQ=255)
                                                                       tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCg                                                                          >  2:4000564/1‑150 (MQ=255)
                                                                         ccACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                        <  1:25819/150‑1 (MQ=255)
                                                                           aCCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                      <  1:2220970/150‑1 (MQ=255)
                                                                           aCCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                      <  1:581601/150‑1 (MQ=255)
                                                                           aCCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                      <  1:4897719/150‑1 (MQ=255)
                                                                            ccGACGGTGATCGGTTTTTACCTGTTGATCGCCCTGGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggtg                                                                     <  2:3378504/150‑1 (MQ=255)
                                                                              gACGGTGATCGGTTTTTAGGTGTTGATCGCGCTTGGCCCGCACGGCTGGCTTGGTCAGGCCAGCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGGCTGGTGATTGGCTCGACGGTGTACTACATGCCGTTCGTGGTGCa                                                                   >  2:2915737/1‑150 (MQ=255)
                                                                                cGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAc                                                                 <  1:3417108/150‑1 (MQ=255)
                                                                                   tGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAcccc                                                              <  1:4517533/150‑1 (MQ=255)
                                                                                      tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGc                                                           <  1:1823588/150‑1 (MQ=255)
                                                                                      tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGc                                                           >  1:3043886/1‑150 (MQ=255)
                                                                                       cGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCa                                                          <  2:3376192/150‑1 (MQ=255)
                                                                                           ttttACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAAt                                                      <  2:789671/150‑1 (MQ=255)
                                                                                           ttttACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAAt                                                      <  2:3449971/150‑1 (MQ=255)
                                                                                                ccTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCtt                                                 >  1:1163684/1‑150 (MQ=255)
                                                                                                 cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                <  1:2315470/150‑1 (MQ=255)
                                                                                                 cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                <  1:2315405/150‑1 (MQ=255)
                                                                                                 cTGTTGATCGCCCTTGGCAAGCACGGCTGGATTTGTCATGCCACCCAGGCGCTAGGCCTTGGCAGTGTGGTGTTCAGCTTCACAGGCCTTGTGATTGGCTCGACGGTGTACTCCATGCAGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                <  1:1621374/150‑1 (MQ=255)
                                                                                                  tggtGATCGCCCTTGGCCCGCACGGCGGGCGTGGTCAGGCCACCCAGGCGCGAGGCCTGGGCAGTGTGGTGTTCAGCGTCACAGGCCTGCTGATTGGCTCGCCGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg                                               <  2:4976794/147‑1 (MQ=255)
                                                                                                  tGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGTCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg                                               <  2:696906/150‑1 (MQ=255)
                                                                                                  tGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg                                               >  1:2772829/1‑150 (MQ=255)
                                                                                                  tGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg                                               <  1:3788650/150‑1 (MQ=255)
                                                                                                  tGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTACAGAATGCCTTTg                                               >  2:2160533/1‑150 (MQ=255)
                                                                                                   gTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTgg                                              >  2:593747/1‑150 (MQ=255)
                                                                                                   gTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTgg                                              >  1:3318213/1‑150 (MQ=255)
                                                                                                    ttGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgc                                             <  2:922346/150‑1 (MQ=255)
                                                                                                      gATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGAGCTAGTCCTGGGCAGGGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgcgc                                           >  2:4617203/1‑150 (MQ=255)
                                                                                                        tCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGa                                         >  2:2113042/1‑150 (MQ=255)
                                                                                                        tCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCCGCAGAATGCCTTTGGCGCGa                                         >  2:2113043/1‑150 (MQ=255)
                                                                                                         cGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGAt                                        >  1:4668057/1‑150 (MQ=255)
                                                                                                           cccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCg                                      >  1:3220531/1‑150 (MQ=255)
                                                                                                           cccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCg                                      <  2:2508099/150‑1 (MQ=255)
                                                                                                           cccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCg                                      <  2:4745751/150‑1 (MQ=255)
                                                                                                           cccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCg                                      >  1:4413082/1‑150 (MQ=255)
                                                                                                            ccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCgg                                     >  2:2828254/1‑150 (MQ=255)
                                                                                                             cTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGc                                    <  2:2338205/150‑1 (MQ=255)
                                                                                                             cTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGc                                    <  2:539864/150‑1 (MQ=255)
                                                                                                               tGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCACGCGCTGCAGAATGCCTTTGGCGCGATCGGCCa                                  <  2:1261025/150‑1 (MQ=255)
                                                                                                                ggCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAg                                 <  1:89049/150‑1 (MQ=255)
                                                                                                                    cGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGcc                             >  2:2291429/1‑150 (MQ=255)
                                                                                                                    cGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGcc                             >  1:245052/1‑150 (MQ=255)
                                                                                                                     gCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGccc                            <  1:548633/150‑1 (MQ=255)
                                                                                                                        cGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCt                         <  1:2575155/150‑1 (MQ=255)
                                                                                                                        cGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCt                         <  1:4665724/150‑1 (MQ=255)
                                                                                                                           ctggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGa                      <  2:66719/150‑1 (MQ=255)
                                                                                                                            tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaa                     <  2:4358924/150‑1 (MQ=255)
                                                                                                                            gggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCCCGATCGGCCAGCGCCCCCTGGaa                     <  1:128181/149‑1 (MQ=255)
                                                                                                                              gctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                           <  1:3261374/78‑1 (MQ=255)
                                                                                                                              gctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                           >  2:3261374/1‑78 (MQ=255)
                                                                                                                                ttGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcg                 >  2:1108835/1‑150 (MQ=255)
                                                                                                                                  ggTCAGGCCACCCAGGCGCAAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGCTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTGTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgcc               <  2:2073872/150‑1 (MQ=255)
                                                                                                                                   gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg              <  1:3086995/150‑1 (MQ=255)
                                                                                                                                      acactgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTa           >  1:3192702/7‑150 (MQ=255)
                                                                                                                                        gCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAct         >  2:1726824/1‑150 (MQ=255)
                                                                                                                                         ccACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCATCGCCCCCTGGAAGTCGCCGCTActc        <  1:1508713/150‑1 (MQ=255)
                                                                                                                                         ccACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActc        <  1:103062/150‑1 (MQ=255)
                                                                                                                                          cACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActct       >  2:1800989/1‑150 (MQ=255)
                                                                                                                                           aCCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTg      >  2:2811913/1‑150 (MQ=255)
                                                                                                                                              cAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGa   <  2:3907805/150‑1 (MQ=255)
                                                                                                                                              cAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGa   >  1:4090544/1‑150 (MQ=255)
                                                                                                                                               aGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGag  >  2:2317344/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAG  >  NC_002947/4356865‑4357157

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