Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,355,090 Δ1 bp 8.0% 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,355,0900G. 8.0% 309.5 / 21.7 100coding (246/681 nt)modBmolybdate ABC transporter permease ModB
Reads supporting (aligned to +/- strand):  ref base G (41/50);  new base . (3/5);  total (44/56)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 102 total aligned reads displayed.

CCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAGCCA  >  NC_002947/4354942‑4355238
                                                                                                                                                    |                                                                                                                                                     
ccTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGgcgc                                                                                                                                                      >  2:1575731/1‑150 (MQ=255)
 gggcgggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCt                                                                                                                                                     >  1:1892355/6‑150 (MQ=255)
  tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                                                                                                                                                    <  1:218923/150‑1 (MQ=255)
  tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCCGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa                                                                                                                                                    <  1:220791/150‑1 (MQ=255)
   ggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGTTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAg                                                                                                                                                   <  1:2308021/150‑1 (MQ=255)
   ggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAg                                                                                                                                                   <  1:880483/150‑1 (MQ=255)
        cTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTCGATCACACCGCGCACGCACGGCTCGCTCAGTAAGGCCACCAAGGCCCTAgcgctg                                                                                                                                              <  2:341912/150‑6 (MQ=255)
         acgcgctCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg                                                                                                                                            >  2:2212979/2‑150 (MQ=255)
          cgcgcgCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTggg                                                                                                                                           <  2:53738/150‑1 (MQ=255)
               gcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAgtg                                                                                                                                      >  2:2224335/1‑150 (MQ=255)
                cACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtggcct                                                                                                                                     <  2:2038318/150‑7 (MQ=255)
                       tctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTc                                                                                                                              <  1:2212979/150‑1 (MQ=255)
                       tctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTc                                                                                                                              >  1:1936766/1‑150 (MQ=255)
                        ctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCa                                                                                                                             <  1:224542/150‑1 (MQ=255)
                         ttGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                    <  2:1055650/142‑1 (MQ=255)
                         ttGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                    >  1:1055650/1‑142 (MQ=255)
                          tGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGc                                                                                                                           <  2:1852065/150‑1 (MQ=255)
                            gCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGGtt                                                                                                                         <  1:2400839/150‑1 (MQ=255)
                                 gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAg                                                                                                                    >  1:2338377/1‑150 (MQ=255)
                                   gggCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGc                                                                                                                  <  2:175045/150‑1 (MQ=255)
                                       cGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTgg                                                                                                              >  2:1753974/1‑150 (MQ=255)
                                         gtgggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCCGGtt                                                                                                            <  2:1261418/150‑2 (MQ=255)
                                            ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGGTCGGCTTTTACCTGTTGATCGCCCTTGGCCCCCACGGCTGGCTAGGTCAGGCCACCCAGCCGCTAGGCCTGGGCCGTGTGGTGTTCAGCTTCACAGGCCGGGTGAtt                                                                                                         >  2:596378/1‑150 (MQ=255)
                                            ggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAtt                                                                                                         <  1:1575731/150‑1 (MQ=255)
                                              tgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtgg                                                                                                                                        <  2:10796/117‑4 (MQ=255)
                                              tgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtgg                                                                                                                                        >  1:10796/1‑114 (MQ=255)
                                               gCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGc                                                                                                      <  1:1359015/150‑1 (MQ=255)
                                                cGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                                     <  1:1952026/150‑1 (MQ=255)
                                                 ggtggtGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGTCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTc                                                                                                    <  2:805797/150‑1 (MQ=255)
                                                  gtgttgGCACTGCTGCTGGTAGTGCCACCGACGGTGATCGGGTTTTACCTGTTGATCGCCCTTGGCCCGCAGGGCTGGCTTGGTTAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCg                                                                                                   >  2:941282/1‑150 (MQ=255)
                                                   tggtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGAGTGGCTCGa                                                                                                  <  2:563544/150‑1 (MQ=255)
                                                    ggtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGAc                                                                                                 <  2:238622/150‑1 (MQ=255)
                                                     gtggCACTGCCGCTGGTACTGCCACCGACGGTGCAAGCTTTTTACCTGGTGTTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACg                                                                                                <  1:541146/150‑1 (MQ=255)
                                                     gtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACg                                                                                                <  2:1555072/150‑1 (MQ=255)
                                                       ggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                              <  2:2338377/150‑1 (MQ=255)
                                                        gCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg                                                                                             >  1:519365/1‑150 (MQ=255)
                                                           cTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTAc                                                                                          <  2:284697/150‑1 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGc                                                                                   <  1:821934/144‑1 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGc                                                                                   >  2:821934/1‑144 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCg                                                                             >  2:581790/1‑150 (MQ=255)
                                                                         gCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgt                                                                            <  2:1936766/150‑1 (MQ=255)
                                                                          ccACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                           <  1:1550586/150‑1 (MQ=255)
                                                                          ccACCGACAGCGAACGAATGTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTGCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                           <  2:519365/150‑1 (MQ=255)
                                                                           cACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtgg                                                                          <  1:2345690/150‑1 (MQ=255)
                                                                           cACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTGCAGCTTCATAGGCCGGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtgg                                                                          <  2:1892355/150‑1 (MQ=255)
                                                                            aCCGCCGGTGATTGGTTGGGTCCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACTCGGGCGCTAGGCCTGGGCAGTGTGGTGTCCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                         <  2:61859/150‑1 (MQ=255)
                                                                            aCCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                         <  1:1721190/150‑1 (MQ=255)
                                                                                  ggTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAcc                                                                   >  2:511763/1‑150 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggtg                                                                                                                                 <  2:607414/84‑1 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggtg                                                                                                                                 <  2:607625/84‑1 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggtg                                                                                                                                 >  1:607414/1‑84 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggtg                                                                                                                                 >  1:607625/1‑84 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggt                                                                                                                                  <  2:2030968/83‑1 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTgtggt                                                                                                                                  >  1:2030968/1‑83 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTg                                                               >  1:859003/1‑150 (MQ=255)
                                                                                        tggTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGTCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCa                                                             <  1:1831911/149‑1 (MQ=255)
                                                                                          gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGa                                                                                                           >  2:51476/1‑102 (MQ=255)
                                                                                          gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGa                                                                                                           <  1:51476/102‑1 (MQ=255)
                                                                                          gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGAGTACTCCAGGCCGTGCGTGGTGCAACCCCTGCAGa                                                           >  2:1831911/1‑150 (MQ=255)
                                                                                            ttttACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAAt                                                         <  1:511763/150‑1 (MQ=255)
                                                                                                cccTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTGAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCt                                                     <  2:1346836/149‑1 (MQ=255)
                                                                                                aCCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCt                                                     >  1:274780/1‑150 (MQ=255)
                                                                                                  cTGTTGATCGCCCTTGGCCCGCCCGGATGGCTTGGTCTGCCGTGCGGTTCGCTAGGCCTGGGCAGTGTGGTGTGCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                   <  1:540020/150‑1 (MQ=255)
                                                                                                  cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                                   >  2:2271384/1‑150 (MQ=255)
                                                                                                       gATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgcgc                                              <  1:2049820/150‑1 (MQ=255)
                                                                                                        aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGcgcg                                             >  2:1754595/1‑150 (MQ=255)
                                                                                                        aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGcgcg                                             >  1:1432712/1‑150 (MQ=255)
                                                                                                         tCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGc                                                              >  1:730450/1‑132 (MQ=255)
                                                                                                         tCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGc                                                              <  2:730450/132‑1 (MQ=255)
                                                                                                          cGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTGGGCGCGAt                                           >  1:103039/1‑150 (MQ=255)
                                                                                                             ccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                           >  2:1429217/1‑115 (MQ=255)
                                                                                                             ccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg                                                                           <  1:1429217/115‑1 (MQ=255)
                                                                                                             ccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCgg                                        >  2:1056971/1‑150 (MQ=255)
                                                                                                             ccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCgg                                        >  2:2073880/1‑150 (MQ=255)
                                                                                                              cTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGc                                       <  2:103039/150‑1 (MQ=255)
                                                                                                                tGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTAATACAGCCCTGGTATGGTGGAATCCCTGCAGAATGCCTTTGGCGCGATCGGCCa                                     <  1:2224335/150‑1 (MQ=255)
                                                                                                                 ggCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAg                                    <  2:859003/150‑1 (MQ=255)
                                                                                                                    ccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgcgc                                 >  2:345618/1‑150 (MQ=255)
                                                                                                                    ccGCACGGCTGGCTTGGTCAGGCCAACGAGTCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgcgc                                 >  1:676895/1‑150 (MQ=255)
                                                                                                                            ctggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                              >  1:2318609/1‑81 (MQ=255)
                                                                                                                            ctggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgt                                                                                              <  2:2318609/81‑1 (MQ=255)
                                                                                                                            ctggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCATGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGa                         >  1:883051/1‑150 (MQ=255)
                                                                                                                               gctTGGTCGGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGCTTGGCTCGACGGTGTACTCCATGCCGTGCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGCTCGGCCAGCGCCCCCTGGAAGt                      <  2:949570/150‑1 (MQ=255)
                                                                                                                                 ttGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgcgc                                              >  1:628481/1‑124 (MQ=255)
                                                                                                                                 ttGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGgcgc                                              <  2:628481/124‑1 (MQ=255)
                                                                                                                                 ttGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcg                    >  1:10976/1‑150 (MQ=255)
                                                                                                                                    gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGAc                                                                                                 >  2:2172326/1‑70 (MQ=255)
                                                                                                                                    gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGAc                                                                                                 <  1:2172326/70‑1 (MQ=255)
                                                                                                                                    gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg                 >  2:1202420/1‑150 (MQ=255)
                                                                                                                                         gCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTa              <  2:549209/148‑1 (MQ=255)
                                                                                                                                         gCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTa              >  1:549209/1‑148 (MQ=255)
                                                                                                                                           tgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg                                                  <  2:889144/108‑1 (MQ=255)
                                                                                                                                           tgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg                                                  >  1:889144/3‑110 (MQ=255)
                                                                                                                                            aCCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActctt         >  2:1973122/1‑149 (MQ=255)
                                                                                                                                             cccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCa                                                                                      <  1:1833034/72‑1 (MQ=255)
                                                                                                                                             cccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCa                                                                                      >  2:1833034/1‑72 (MQ=255)
                                                                                                                                              ccAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCg       >  2:954175/1‑150 (MQ=255)
                                                                                                                                              ccAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCg       >  2:876327/1‑150 (MQ=255)
                                                                                                                                               cAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGa      <  1:1753974/150‑1 (MQ=255)
                                                                                                                                                   cgcTAGCCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACGCCATGCCGGTCGTGGTGCAACCCCTGCACAATGCCGTTGGCGCGATTGCCCAGCGCCCCCTGGAAGTCGCCGCTACTCCGCGagcca  <  1:1384757/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGAGCCA  >  NC_002947/4354942‑4355238

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