Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,003,385 Δ1 bp coding (186/681 nt) ycbF → putative chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,003,3840A.94.3% 769.2 / 31.7 194coding (185/681 nt)ycbFputative chaperone
Reads supporting (aligned to +/- strand):  ref base A (2/0);  major base . (122/59);  minor base G (5/6);  total (129/65)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 213 total aligned reads displayed.

GAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGA  >  NC_000913/1003293‑1003472
                                                                                           |                                                                                        
gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGag                                                                                         <  2:5162813/93‑2 (MQ=255)
gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGag                                                                                         >  2:1458431/1‑92 (MQ=255)
  aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGGTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGagcg                                                                                       >  1:162525/1‑90 (MQ=255)
       tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGAc                                                                                 <  2:1004548/93‑1 (MQ=255)
            aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGc                                                                            >  1:881583/1‑93 (MQ=255)
             ttttACTGGCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGcc                                                                           <  1:967127/93‑1 (MQ=255)
               ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAg                                                                         >  1:835822/1‑93 (MQ=255)
                tACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGc                                                                        >  2:2716247/1‑93 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCAcc                                                                     <  2:1550778/93‑1 (MQ=255)
                   tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGCGGGGGAGCGTCGACAGTGCCAGCAcc                                                                     >  1:3116703/1‑93 (MQ=255)
                    gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCa                                                                    >  2:2984671/1‑93 (MQ=255)
                     aCTAATACCAGTCAGGAATGCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGTCAGCAccgt                                                                   >  1:6564190/1‑91 (MQ=255)
                     aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCAt                                                                   >  1:231421/1‑93 (MQ=255)
                      cTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCAtt                                                                  >  1:4587025/1‑93 (MQ=255)
                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTAc                                                                <  1:4387813/93‑1 (MQ=255)
                        aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTAc                                                                >  2:3605374/1‑93 (MQ=255)
                         aTACGAGGCAGGAATCCTGTCTTATTAACAGGAAAATCAACCGCCCAACGGGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACt                                                               <  2:724472/93‑1 (MQ=255)
                         aTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTCCt                                                               >  2:4017307/1‑93 (MQ=255)
                         aTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACt                                                               >  2:1171257/1‑93 (MQ=255)
                         aTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACt                                                               >  1:5874588/1‑93 (MQ=255)
                          tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTg                                                              >  2:1318059/1‑93 (MQ=255)
                          aaCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTg                                                              <  2:5206442/92‑1 (MQ=255)
                            ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGc                                                            >  2:5406188/1‑93 (MQ=255)
                              aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCg                                                          >  2:1748673/1‑93 (MQ=255)
                                 cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTc                                                       >  1:3484891/1‑93 (MQ=255)
                                  aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                      <  1:1041930/93‑1 (MQ=255)
                                  aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                      <  1:5427779/93‑1 (MQ=255)
                                  aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                      <  1:4218759/93‑1 (MQ=255)
                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                     >  2:281113/1‑93 (MQ=255)
                                   ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                     >  1:3786426/1‑93 (MQ=255)
                                     aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                   >  1:4711496/1‑93 (MQ=255)
                                     aaTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                   >  2:4868599/1‑93 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                  >  1:737569/1‑93 (MQ=255)
                                      aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCa                                                  >  1:1595307/1‑93 (MQ=255)
                                       tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                 <  1:2177420/93‑1 (MQ=255)
                                       tCCTGGATTATTAACAGTAAAATCAACCGCCCACCGCTTGGGCCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                 <  1:1866642/93‑1 (MQ=255)
                                        ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                <  1:3382823/93‑1 (MQ=255)
                                         cTGGTTTATCAACAGAAAAATCACCCGCCAAACGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                               <  2:1616561/93‑1 (MQ=255)
                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                               <  1:3256282/93‑1 (MQ=255)
                                         cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                               >  1:6640552/1‑93 (MQ=255)
                                          tgttttattaAAAGAAAAACCACCCGCCCAACGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                              <  1:5165959/89‑1 (MQ=255)
                                          tGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                              >  2:730055/1‑93 (MQ=255)
                                          tGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                              >  1:1778206/1‑93 (MQ=255)
                                           gggttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                             >  1:5511125/4‑93 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                             >  2:814579/1‑93 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                             <  1:6694697/93‑1 (MQ=255)
                                           ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                             >  2:6595449/1‑93 (MQ=255)
                                            gCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAtt                                            <  2:1837325/93‑1 (MQ=255)
                                             cTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtc                                           >  2:6351200/1‑93 (MQ=255)
                                               tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                          >  2:2293889/1‑93 (MQ=255)
                                               taataACAAAAAAATCACCCGCCAACCGGTTGGGCCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                         <  1:417273/90‑1 (MQ=255)
                                                attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGCGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                        >  1:1907431/1‑93 (MQ=255)
                                                 ttaCCAGGAAAACCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                       <  1:2141700/93‑1 (MQ=255)
                                                 ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                       >  2:3294302/1‑93 (MQ=255)
                                                 ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                       >  2:5924672/1‑93 (MQ=255)
                                                  taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                      >  2:1504107/1‑93 (MQ=255)
                                                   aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                     >  1:1091642/1‑93 (MQ=255)
                                                   aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                     >  2:4591456/1‑93 (MQ=255)
                                                   aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                     <  2:6257044/93‑1 (MQ=255)
                                                   aaCAGTAAAATCAACCCCCCACCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                     <  2:5377254/93‑1 (MQ=255)
                                                    aCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                    >  2:6565775/1‑93 (MQ=255)
                                                     cAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAg                                   >  2:1128500/1‑93 (MQ=255)
                                                     cAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTCATTCTCCTGAAg                                   >  1:5648384/1‑93 (MQ=255)
                                                      aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                  >  1:6368141/1‑93 (MQ=255)
                                                       gaaaaaTCAACCGCCCAACGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                 <  1:4330481/91‑1 (MQ=255)
                                                       gTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                                 >  1:2513491/1‑93 (MQ=255)
                                                       aaaaaaTCACCCGCCCAACGCTTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTACTGAAGcc                                 <  2:6133648/91‑1 (MQ=255)
                                                        tAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                <  2:5070212/93‑1 (MQ=255)
                                                        tAAAATCAACCGCCCAACGCGTTGGGCCTGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                <  2:852565/93‑1 (MQ=255)
                                                        tAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                <  1:4212736/93‑1 (MQ=255)
                                                          aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                              >  1:3332003/1‑93 (MQ=255)
                                                             tCAACCGCCCACCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                           >  2:6718213/1‑93 (MQ=255)
                                                             tCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                           >  1:2845897/1‑93 (MQ=255)
                                                               aaCCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACgc                         >  2:3352517/1‑92 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                        <  1:2043264/93‑1 (MQ=255)
                                                                aCCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                        >  2:3037174/1‑93 (MQ=255)
                                                                  cGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                      >  2:2462721/1‑93 (MQ=255)
                                                                  cGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                      <  2:773850/93‑1 (MQ=255)
                                                                  cGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                      >  1:4157108/1‑93 (MQ=255)
                                                                  cGCCCAACGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                      <  1:6297241/93‑1 (MQ=255)
                                                                   gCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                     >  1:4640193/1‑93 (MQ=255)
                                                                     ccACCGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                   <  1:687887/93‑1 (MQ=255)
                                                                     ccAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                   >  2:2111748/1‑93 (MQ=255)
                                                                      cAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCAc                  >  2:6691689/1‑93 (MQ=255)
                                                                        aCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                <  2:4284864/93‑1 (MQ=255)
                                                                        aCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                <  2:199769/93‑1 (MQ=255)
                                                                        aCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                <  1:1971759/93‑1 (MQ=255)
                                                                           gggcgggCCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg              <  1:1169280/89‑1 (MQ=255)
                                                                           ctgggggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc             <  1:1440355/89‑1 (MQ=255)
                                                                              tGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt          <  2:2797489/93‑1 (MQ=255)
                                                                              tGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt          <  1:1069533/93‑1 (MQ=255)
                                                                              ggggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt          <  1:3664409/92‑1 (MQ=255)
                                                                               gggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt         <  1:1644506/93‑1 (MQ=255)
                                                                               gggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt         <  2:5467603/93‑1 (MQ=255)
                                                                               gggCCAGGGGGGAGTGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt         <  1:5043748/93‑1 (MQ=255)
                                                                                ggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg        >  2:6437004/1‑93 (MQ=255)
                                                                                   ccggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt      >  2:6579322/3‑93 (MQ=255)
                                                                                     ggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga    >  1:4006806/1‑93 (MQ=255)
                                                                                      ggggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgag   >  1:5281566/1‑93 (MQ=255)
                                                                                       ggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgaga  >  1:271320/1‑93 (MQ=255)
                                                                                           |                                                                                        
GAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGA  >  NC_000913/1003293‑1003472

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