Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,355,090 Δ1 bp 7.5% 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. 7.5% 247.6 / 15.0 80coding (246/681 nt)modBmolybdate ABC transporter permease ModB
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base . (4/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGA  >  NC_002947/4354942‑4355234
                                                                                                                                                    |                                                                                                                                                  
ccTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGgcgc                                                                                                                                                   >  2:151614/1‑150 (MQ=255)
   ggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAg                                                                                                                                                <  2:1068424/150‑1 (MQ=255)
         tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGTACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg                                                                                                                                         <  1:396366/150‑1 (MQ=255)
         tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg                                                                                                                                         <  2:994059/150‑1 (MQ=255)
           gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggt                                                                                                                                       >  1:345671/1‑149 (MQ=255)
           gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTCGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGc                                                                                                                                       >  1:172628/1‑150 (MQ=255)
              cgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAgt                                                                                                                                    <  2:320655/150‑1 (MQ=255)
               ttccgatctctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAgtg                                                                                                                                   <  2:1826664/143‑1 (MQ=255)
                cACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCCTAGCGCCTgggtggcct                                                                                                                                  <  2:345671/150‑7 (MQ=255)
                 aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                 <  2:1673468/150‑1 (MQ=255)
                      ctctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCAACCCAGGCGCTAGGCCTGGGCAGTGTGGTGtt                                                                                                                            >  1:1735702/1‑150 (MQ=255)
                         ttGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAg                                                                                                                         <  1:260031/150‑1 (MQ=255)
                           gtcggcgcGGGCTGGTGTGTTCGCAGGTGTCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTTATCGCCCTTGGCCCGCACGGCTCGCTTGGTCAGGCCACCCATGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCt                                                                                                                       <  2:806803/146‑1 (MQ=255)
                            gCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCtt                                                                                                                      >  2:1428147/1‑150 (MQ=255)
                              tGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTca                                                                                                                    >  1:236849/1‑150 (MQ=255)
                                 gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAg                                                                                                                 <  2:236849/150‑1 (MQ=255)
                                    ggCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGcc                                                                                                              <  1:1078759/150‑1 (MQ=255)
                                         gtgggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTg                                                                                                         >  1:458921/1‑150 (MQ=255)
                                                cGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                                  <  2:192390/150‑1 (MQ=255)
                                                cGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCt                                                                                                  <  2:1215037/150‑1 (MQ=255)
                                                  gtggtgGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGACAAGTGTGGTGTTCAGCTTCACAGGCCTGGAGAtcggaaga                                                                                                >  1:18505/1‑143 (MQ=255)
                                                         cACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCGCGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtgt                                                                                         >  2:422328/1‑150 (MQ=255)
                                                          aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGATTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTa                                                                                        >  1:18013/1‑150 (MQ=255)
                                                          aCTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCACTTCGGAAGCACGGCTGGCTGAGCTAGGGTACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTa                                                                                        >  1:1206734/1‑150 (MQ=255)
                                                             gCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTc                                                                                     <  2:428960/150‑1 (MQ=255)
                                                             gCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTCGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTc                                                                                     <  2:172628/150‑1 (MQ=255)
                                                              ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                 <  2:1372246/105‑1 (MQ=255)
                                                              ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGtgtg                                                                                                                                 >  1:1372246/1‑105 (MQ=255)
                                                              ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTgg                                                                                                           <  1:1225192/127‑1 (MQ=255)
                                                              ccGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTgg                                                                                                           >  2:1225192/1‑127 (MQ=255)
                                                                  tGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCAACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGc                                                                                <  2:1735702/150‑1 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCt                                                                                                                                                  <  1:871945/79‑1 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCt                                                                                                                                                  >  2:871945/1‑79 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg                                                                                                                                         >  2:265217/1‑87 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTgg                                                                                                                                         <  1:265217/87‑1 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAg                                                                                                                 >  2:1846656/1‑111 (MQ=255)
                                                                        tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAg                                                                                                                 <  1:1846656/111‑1 (MQ=255)
                                                                            cccGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                     <  2:1786293/149‑1 (MQ=255)
                                                                             ccGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggtg                                                                    <  2:1128045/150‑1 (MQ=255)
                                                                                  ggTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAgt                                                                                                                                    >  2:268600/1‑82 (MQ=255)
                                                                                  ggTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAgt                                                                                                                                    <  1:268600/82‑1 (MQ=255)
                                                                                    tGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAAcccc                                                             <  2:458921/150‑1 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTg                                                                                                            <  1:656047/102‑1 (MQ=255)
                                                                                      aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTg                                                                                                            >  2:656047/1‑102 (MQ=255)
                                                                                       tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGc                                                          >  1:833251/1‑150 (MQ=255)
                                                                                        cGGTTTTTACCTGTTGATCGCCCTTGGCCCGCGCGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCa                                                         <  1:422328/150‑1 (MQ=255)
                                                                                            ttttACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAAt                                                     <  1:1857394/150‑1 (MQ=255)
                                                                                             tttACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCACCCCCTGCAGAATg                                                    >  2:1019089/1‑150 (MQ=255)
                                                                                                aCCTGTTGATCGCCATTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCt                                                 >  2:1563030/1‑150 (MQ=255)
                                                                                                  cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGGGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGCTTGGCACGAAGGTGTGCTCCACGCCGGTTAGTGGGCATCTCCGTCACAATGCCttt                                               >  2:211669/1‑150 (MQ=255)
                                                                                                  cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                               <  1:1620869/150‑1 (MQ=255)
                                                                                                    gTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGtt                                                                           >  2:886146/1‑120 (MQ=255)
                                                                                                    gTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGtt                                                                           <  1:886146/120‑1 (MQ=255)
                                                                                                     ttGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGa                                                                                                        >  2:599734/1‑91 (MQ=255)
                                                                                                     ttGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGa                                                                                                        <  1:599734/91‑1 (MQ=255)
                                                                                                        aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGcgcg                                         >  1:160710/1‑150 (MQ=255)
                                                                                                         tCGCCCTTGGCCCGCACTGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGGGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGa                                        >  1:802052/1‑150 (MQ=255)
                                                                                                          cGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGAGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGAt                                       >  2:1780630/1‑150 (MQ=255)
                                                                                                              cTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGc                                   <  2:160710/150‑1 (MQ=255)
                                                                                                                  gCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGGGCTAGGCCTAGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgc                               <  1:211669/150‑1 (MQ=255)
                                                                                                                  gCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGTCCTGGTGATTGGCTCGAc                                                                                              <  2:1627430/88‑1 (MQ=255)
                                                                                                                  gCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGAc                                                                                              >  1:1627430/1‑88 (MQ=255)
                                                                                                                   cccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgcg                              >  1:197939/1‑150 (MQ=255)
                                                                                                                    ccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                               >  1:376993/1‑132 (MQ=255)
                                                                                                                    ccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt                                               <  2:376993/132‑1 (MQ=255)
                                                                                                                             tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGAc                                                                                              <  2:1315715/77‑1 (MQ=255)
                                                                                                                             tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGAc                                                                                              >  1:1315715/1‑77 (MQ=255)
                                                                                                                             tggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaa                    <  1:287032/150‑1 (MQ=255)
                                                                                                                              ggctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAg                   >  2:1039979/1‑150 (MQ=255)
                                                                                                                                ctTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTc                 >  2:1888022/1‑150 (MQ=255)
                                                                                                                                 ttGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGaacaag                >  2:226/1‑146 (MQ=255)
                                                                                                                                 ttGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcg                >  2:124113/1‑150 (MQ=255)
                                                                                                                                  tGGGCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCTCCCCCTGGAAGTcgc               <  1:794293/150‑1 (MQ=255)
                                                                                                                                    gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg             >  2:1222914/1‑150 (MQ=255)
                                                                                                                                    gTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccg             >  2:1473739/1‑150 (MQ=255)
                                                                                                                                        cactgcccAGGCCTAGCGCCTGGGTAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTAc         >  2:874051/6‑150 (MQ=255)
                                                                                                                                         actgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                     >  1:891084/5‑89 (MQ=38)
                                                                                                                                         actgcccAGGCCTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtggt                                                                     <  2:891084/85‑1 (MQ=38)
                                                                                                                                          ctgcccAGGCCTAGCGCCTGGGTAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActc       >  2:1516598/4‑150 (MQ=255)
                                                                                                                                               cAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGa  <  1:1428147/150‑1 (MQ=255)
                                                                                                                                               cAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGGGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTCGCGCGATCGGCCAGCGCCCCCTGGAACTCGCCGATACTCTGCgg  >  2:862216/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGA  >  NC_002947/4354942‑4355234

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