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

ATTCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAA  >  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)
 ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtgt                                                                                                                                              <  7:81198/148‑1 (MQ=255)
 ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtgt                                                                                                                                              <  5:119684/148‑1 (MQ=255)
 ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGgtgtgt                                                                                                                                              <  5:277431/148‑1 (MQ=255)
    ggcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTggg                                                                                                                                           <  5:19130/148‑1 (MQ=255)
    ggcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTgg                                                                                                                                            <  7:139181/147‑1 (MQ=255)
     gcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTgggg                                                                                                                                          >  3:180835/1‑148 (MQ=255)
       ggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCa                                                                                                                                        <  7:210002/148‑1 (MQ=255)
           aCGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACg                                                                                                                                     <  3:170619/147‑1 (MQ=255)
            cGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGt                                                                                                                                    >  8:341283/1‑147 (MQ=255)
              cTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGaa                                                                                                                                 <  6:74039/148‑1 (MQ=255)
              cTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGGGAATCATGATGACAAAAATGAGGCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGaa                                                                                                                                 <  3:37912/148‑1 (MQ=255)
                gTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAAcc                                                                                                                               >  7:168373/1‑148 (MQ=255)
                gTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAAc                                                                                                                                >  2:178890/1‑147 (MQ=255)
                 tACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTAAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAAcc                                                                                                                               <  3:259058/147‑1 (MQ=255)
                              aTTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTag                                                                                                                    >  5:151394/1‑145 (MQ=255)
                                 ttCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAgaagaa                                                                                                               <  4:148698/147‑1 (MQ=255)
                                  tCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTACGGGGTGAACCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              <  4:168679/147‑1 (MQ=255)
                                  tCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              <  4:180835/147‑1 (MQ=255)
                                     gggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGAGGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGaa                                                                                                          <  2:384205/148‑1 (MQ=255)
                                      ggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGaa                                                                                                          <  1:306074/147‑1 (MQ=255)
                                      ggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCGGGCGATGTTTCTCGCCGGGTGGGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGaa                                                                                                          <  1:248063/147‑1 (MQ=255)
                                      ggggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGGGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGa                                                                                                           <  7:118438/146‑1 (MQ=255)
                                       gggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGtgtgtg                                                                                                                                             <  3:347198/111‑1 (MQ=255)
                                       gggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGtgtgtg                                                                                                                                             >  4:347198/1‑111 (MQ=255)
                                         gcgcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCa                                                                                                      <  1:178890/148‑1 (MQ=255)
                                           gcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAAc                                                                                                                                <  3:265517/120‑1 (MQ=255)
                                           gcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAAc                                                                                                                                >  4:265517/1‑120 (MQ=255)
                                              caAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              <  7:104164/135‑1 (MQ=255)
                                              caAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGaagaag                                                                                                              >  8:104164/1‑135 (MQ=255)
                                                 aGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGa                                                                                              <  2:251660/148‑1 (MQ=255)
                                                      gACTTTATCTATATTTGTATGGGGTCAATCTTGATGACAAATATGAGTCGCTACGCCTTGATTACCGAGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGaacaac                                                                                         <  7:298168/148‑1 (MQ=255)
                                                        cTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGACGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCa                                                                                       <  8:168373/148‑1 (MQ=255)
                                                           tATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGcc                                                                                    >  2:390073/1‑148 (MQ=255)
                                                           tATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGcc                                                                                    >  3:238352/1‑148 (MQ=255)
                                                             tCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGTCg                                                                                   <  5:266377/147‑1 (MQ=255)
                                                             tCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCg                                                                                   <  1:325858/147‑1 (MQ=255)
                                                                 tatTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTTGGGGCAACGTGAACCTGCACCGGGAGAAGAAGTGAACCCAGCGCCGGAACAACCCGCCGAGc                                                                                >  3:236270/1‑147 (MQ=255)
                                                                 tatTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGcc                                                                               <  4:180538/147‑1 (MQ=255)
                                                                 tatTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACAGCGCTGGCGATGTTTCTAGCCGGGTGTGTTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGc                                                                                >  3:231429/1‑147 (MQ=255)
                                                                                   cTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCca                                                              >  7:360827/1‑146 (MQ=255)
                                                                                     tgatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             <  5:25462/145‑1 (MQ=255)
                                                                                     tgatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             >  6:25462/1‑145 (MQ=255)
                                                                                      gatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcaca                                                            <  7:341283/145‑1 (MQ=255)
                                                                                      gatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             <  1:390073/144‑1 (MQ=255)
                                                                                      gatgaCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCcac                                                             <  4:238352/144‑1 (MQ=255)
                                                                                          aCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCt                                                     >  8:244895/1‑148 (MQ=255)
                                                                                              aaaTGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGt                                                 >  5:340094/1‑148 (MQ=255)
                                                                                               aaTGAGGCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGTAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTg                                                <  8:154165/148‑1 (MQ=255)
                                                                                                       gCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGAt                                        <  2:212727/148‑1 (MQ=255)
                                                                                                         tACGCGTTGTTTAGCGGGCTGGCGATGTTTGTCGCCGGGTGTGTGGGGCAACGGGAACCTGCACCGGTAGAAGAAGTGATACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATcc                                      <  8:360827/148‑1 (MQ=255)
                                                                                                               ttGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCcaca                                                                            >  3:277540/1‑104 (MQ=255)
                                                                                                               ttGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCcaca                                                                            <  4:277540/104‑1 (MQ=255)
                                                                                                               ttGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCgcag                                 >  5:304696/1‑147 (MQ=255)
                                                                                                                  aTTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCagcagc                             >  3:16793/1‑148 (MQ=255)
                                                                                                                   ttACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCg                                                   >  7:255196/1‑125 (MQ=255)
                                                                                                                   ttACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCg                                                   <  8:255196/125‑1 (MQ=255)
                                                                                                                   ttACCGCGCTGGCGATGTTTCTCGCCGGGGGGGTGGGGCAACGGGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGcc                            <  6:340094/148‑1 (MQ=255)
                                                                                                                      ccGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAgg                         <  7:244895/148‑1 (MQ=255)
                                                                                                                             ggCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAAtt                  >  5:191373/1‑148 (MQ=255)
                                                                                                                               cGATGTTTCTCGCCGGGTGTGTGGGGCACCGTGAACCTGCACCGCTAGAAGAAGTGAAACAAGAGCCGGAACAACCAGCCGCGCCACAAAAGCCTGTCACCACAGTGCCCTCTGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTgc                >  6:197068/1‑147 (MQ=255)
                                                                                                                                gATGTTTCTCGCCGGGTGGGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAg               <  4:330294/148‑1 (MQ=255)
                                                                                                                                       ctcGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGaa         <  4:16793/147‑1 (MQ=255)
                                                                                                                                          gCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGa       <  7:252057/146‑1 (MQ=255)
                                                                                                                                            cGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCaa   <  6:191373/148‑1 (MQ=255)
                                                                                                                                             gcgtgtgtGGGGCAACGTGAACCTGCACCGGTAGAAGAAGAGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCACGCCCAAATGAGCACGAAGATCaaa  >  4:57017/3‑148 (MQ=255)
                                                                                                                                              ggTGTGTTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCaaa  <  4:236270/148‑1 (MQ=255)
                                                                                                                                              ggTGTGTTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCaaa  <  4:231429/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATTCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTCGCCGGGTGTGTGGGGCAACGTGAACCTGCACCGGTAGAAGAAGTGAAACCAGCGCCGGAACAACCAGCCGAGCCACAACAGCCTGTCCCCACAGTGCCCTCGGTGCCGACGATCCCGCAGCAGCCAGGCCCAATTGAGCACGAAGATCAAA  >  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: