Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,162,698:1 +T 6.6% coding (61/642 nt) lpoB → OM lipoprotein stimulator of MrcB transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,162,6981.T 6.6% 233.9 / 11.7 61coding (61/642 nt)lpoBOM lipoprotein stimulator of MrcB transpeptidase
Reads supporting (aligned to +/- strand):  ref base . (19/38);  new base T (2/2);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATTCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAA  >  NC_000913/1162551‑1162839
                                                                                                                                                    |                                                                                                                                             
atacggcggTAACGATGTACGGCAGCGCCAATTTTCTGGGCGCTCACACAGTCAGCCTTACTCTATATTTGTAAGGGGTGAATCTTGATTACAAAAATTAGTCTATACGCCTAGATTACCGCGCTTGATATGTTTCTCGCCGGgtgtg                                                                                                                                                <  4:220051/145‑1 (MQ=255)
aTTCGGCGGTAAGGCTGTACGGCAGCGGCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGCGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtg                                                                                                                                                <  7:73099/148‑1 (MQ=255)
aTTCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtg                                                                                                                                                >  3:148698/1‑148 (MQ=255)
 ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtg‑t                                                                                                                                              <  7:81198/148‑1 (MQ=255)
 ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtg‑t                                                                                                                                              <  5:119684/148‑1 (MQ=255)
 ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtg‑t                                                                                                                                              <  5:277431/148‑1 (MQ=255)
    ggcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑Tggg                                                                                                                                           <  5:19130/148‑1 (MQ=255)
    ggcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑Tgg                                                                                                                                            <  7:139181/147‑1 (MQ=255)
     gcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑Tgggg                                                                                                                                          >  3:180835/1‑148 (MQ=255)
       ggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCa                                                                                                                                        <  7:210002/148‑1 (MQ=255)
           aCGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACg                                                                                                                                     <  3:170619/147‑1 (MQ=255)
            cGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGt                                                                                                                                    >  8:341283/1‑147 (MQ=255)
              cTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGaa                                                                                                                                 <  6:74039/148‑1 (MQ=255)
              cTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGGGAATCATGATGACAAAAATGAGGCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGaa                                                                                                                                 <  3:37912/148‑1 (MQ=255)
                gTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAAcc                                                                                                                               >  7:168373/1‑148 (MQ=255)
                gTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAAc                                                                                                                                >  2:178890/1‑147 (MQ=255)
                 tACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTAAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAAcc                                                                                                                               <  3:259058/147‑1 (MQ=255)
                              aTTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTag                                                                                                                    >  5:151394/1‑145 (MQ=255)
                                 ttCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAgaagaa                                                                                                               <  4:148698/147‑1 (MQ=255)
                                  tCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTACGGGGTGAACCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              <  4:168679/147‑1 (MQ=255)
                                  tCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              <  4:180835/147‑1 (MQ=255)
                                     gggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGAGGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGaa                                                                                                          <  2:384205/148‑1 (MQ=255)
                                      ggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGaa                                                                                                          <  1:306074/147‑1 (MQ=255)
                                      ggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCGGGCGATGTTTCTCGCCGGGTGGG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGaa                                                                                                          <  1:248063/147‑1 (MQ=255)
                                      ggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGGGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGa                                                                                                           <  7:118438/146‑1 (MQ=255)
                                       gggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGtgtg‑tg                                                                                                                                             <  3:347198/111‑1 (MQ=255)
                                       gggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGtgtg‑tg                                                                                                                                             >  4:347198/1‑111 (MQ=255)
                                         gcgcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCa                                                                                                      <  1:178890/148‑1 (MQ=255)
                                           gcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAAc                                                                                                                                <  3:265517/120‑1 (MQ=255)
                                           gcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAAc                                                                                                                                >  4:265517/1‑120 (MQ=255)
                                              caAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              <  7:104164/135‑1 (MQ=255)
                                              caAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              >  8:104164/1‑135 (MQ=255)
                                                 aGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGa                                                                                              <  2:251660/148‑1 (MQ=255)
                                                      gACTTTATCTATATTTGTATGGGGTCAATCTTGATGACAAATATGAGTCGCTACGCCTTGATTACCGAGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGaacaac                                                                                         <  7:298168/148‑1 (MQ=255)
                                                        cTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGACGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCa                                                                                       <  8:168373/148‑1 (MQ=255)
                                                           tATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGcc                                                                                    >  2:390073/1‑148 (MQ=255)
                                                           tATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGcc                                                                                    >  3:238352/1‑148 (MQ=255)
                                                             tCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGTCg                                                                                   <  5:266377/147‑1 (MQ=255)
                                                             tCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCg                                                                                   <  1:325858/147‑1 (MQ=255)
                                                                 tatTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTTGGGGCAACGTGAACCTGCACCGGGAGAAGAAGTGAACCCAGCGCCGGAACAACCCGCCGAGc                                                                                >  3:236270/1‑147 (MQ=255)
                                                                 tatTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGcc                                                                               <  4:180538/147‑1 (MQ=255)
                                                                 tatTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACAGCGCTGGCGATGTTTCTAGCCGGGTGTGTTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGc                                                                                >  3:231429/1‑147 (MQ=255)
                                                                                   cTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCca                                                              >  7:360827/1‑146 (MQ=255)
                                                                                     tgatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             <  5:25462/145‑1 (MQ=255)
                                                                                     tgatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             >  6:25462/1‑145 (MQ=255)
                                                                                      gatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcaca                                                            <  7:341283/145‑1 (MQ=255)
                                                                                      gatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             <  1:390073/144‑1 (MQ=255)
                                                                                      gatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             <  4:238352/144‑1 (MQ=255)
                                                                                          aCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCt                                                     >  8:244895/1‑148 (MQ=255)
                                                                                              aaaTGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGt                                                 >  5:340094/1‑148 (MQ=255)
                                                                                               aaTGAGGCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGTAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTg                                                <  8:154165/148‑1 (MQ=255)
                                                                                                       gCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGAt                                        <  2:212727/148‑1 (MQ=255)
                                                                                                         tACGCGTTGTTTAGCGGGCTGGCGATGTTTGTCGCCGGGTGTG‑TGGGGCAACGGGAACCTGCACCGGTAGAAGAAGTGATACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATcc                                      <  8:360827/148‑1 (MQ=255)
                                                                                                               ttGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCcaca                                                                            >  3:277540/1‑104 (MQ=255)
                                                                                                               ttGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCcaca                                                                            <  4:277540/104‑1 (MQ=255)
                                                                                                               ttGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCgcag                                 >  5:304696/1‑147 (MQ=255)
                                                                                                                  aTTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCagcagc                             >  3:16793/1‑148 (MQ=255)
                                                                                                                   ttACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCg                                                   >  7:255196/1‑125 (MQ=255)
                                                                                                                   ttACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCg                                                   <  8:255196/125‑1 (MQ=255)
                                                                                                                   ttACCGCGCTGGCGATGTTTCTCGCCGGGGGGG‑TGGGGCAACGGGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGcc                            <  6:340094/148‑1 (MQ=255)
                                                                                                                      ccGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAgg                         <  7:244895/148‑1 (MQ=255)
                                                                                                                             ggCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAAtt                  >  5:191373/1‑148 (MQ=255)
                                                                                                                               cGATGTTTCTCGCCGGGTGTG‑TGGGGCACCGTGAACCTGCACCGCTAGAAGAAGTGAAACAAGAGCCGGAACAACCAGCCGCGCCACAAAAGCCTGTCACCACAGTGCCCTCTGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTgc                >  6:197068/1‑147 (MQ=255)
                                                                                                                                gATGTTTCTCGCCGGGTGGG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAg               <  4:330294/148‑1 (MQ=255)
                                                                                                                                       ctcGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGaa         <  4:16793/147‑1 (MQ=255)
                                                                                                                                          gCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGa       <  7:252057/146‑1 (MQ=255)
                                                                                                                                            cGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCaa   <  6:191373/148‑1 (MQ=255)
                                                                                                                                             gcgtgtg‑tGGGGCAACGTGAACCTGCACCGGTAGAAGAAGAGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCACGCCCAAATGAGCACGAAGATCaaa  >  4:57017/3‑148 (MQ=255)
                                                                                                                                              ggTGTGTTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCaaa  <  4:236270/148‑1 (MQ=255)
                                                                                                                                              ggTGTGTTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCaaa  <  4:231429/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATTCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTG‑TGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAA  >  NC_000913/1162551‑1162839

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑