Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 976,769 G→A 6.1% L148L (CTG→CTA)  mukB → chromosome condensin MukBEF, ATPase and DNA‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913976,7690GA 6.1% 207.6 / 3.1 65L148L (CTG→CTA) mukBchromosome condensin MukBEF, ATPase and DNA‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base A (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTGCGCTCCGCATCTGACCGTA  >  NC_000913/976622‑976908
                                                                                                                                                   |                                                                                                                                           
tgtggtgGTCGGTGTGGGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCTTTTGCCATTCTTGGACTGCCGATGTCTGTGCAGCCGACACAGCTGGTGACCGACGCCCTGAACGAACGCAAGGCGCGCGtgctgc                                                                                                                                            <  4:429519/148‑1 (MQ=255)
     tggtCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCtt                                                                                                                                       <  4:245077/149‑1 (MQ=255)
        tCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAAc                                                                                                                                    <  1:612083/149‑1 (MQ=255)
             gtgCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCt                                                                                                                               <  3:824716/149‑1 (MQ=255)
             gtgCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCt                                                                                                                               <  1:672637/149‑1 (MQ=255)
               gCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGa                                                                                                                             >  2:26201/1‑149 (MQ=255)
                 gTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCCTGCTGCCGCTTAACGAGCTGaaa                                                                                                                           <  3:681425/149‑1 (MQ=255)
                      cAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACaa                                                                                                                      >  2:192467/1‑149 (MQ=255)
                      cAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACaa                                                                                                                      <  1:521/149‑1 (MQ=255)
                       aaCAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAg                                                                                                                     >  2:548820/1‑149 (MQ=255)
                       aaCAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTACCGCTTAACGAGCTGAAAGACAAg                                                                                                                     >  3:852826/1‑149 (MQ=255)
                       aaCAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTACCGCTTAACGAGCTGAAAGACAAg                                                                                                                     >  3:757001/1‑149 (MQ=255)
                             ttGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAg                                                                                                               <  1:226614/149‑1 (MQ=255)
                                 cGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGa                                                                                                           <  4:327841/149‑1 (MQ=255)
                                        gATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGaa                                                                                                      <  1:234552/147‑1 (MQ=255)
                                        gATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAgg                                                                                                    >  4:188418/1‑149 (MQ=255)
                                        gATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGATGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTCAACGAGCTGAAAGACAAGCTCGAGGCGATGGaa                                                                                                      >  2:234552/1‑147 (MQ=255)
                                         aTCGTAAAGTGTATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTCACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAgcc                                                                                                   <  4:634295/149‑3 (MQ=255)
                                          tCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCg                                                                                                  >  3:128172/1‑149 (MQ=255)
                                            gTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTg                                                                                                >  1:439235/1‑149 (MQ=255)
                                            gTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTg                                                                                                >  2:198627/1‑149 (MQ=255)
                                               aaGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAg                                                                                             >  4:858305/1‑149 (MQ=255)
                                                  tGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGttt                                                                                          <  1:211261/149‑1 (MQ=255)
                                                     atatCAAGCCGTTTGCCATTCATGGACTGCCGATGTCGTTGTAGCCGTCACAGCTGGTGGCCGAAACCCTGAACGAACGCCAGGCGGGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTaaa                                                                                       <  2:404326/149‑1 (MQ=255)
                                                     atatCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTaaa                                                                                       <  3:489417/149‑1 (MQ=255)
                                                       atCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACa                                                                                     <  4:211954/149‑1 (MQ=255)
                                                        tCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAg                                                                                    >  4:6643/1‑149 (MQ=255)
                                                          aaGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGATTAAACAGtt                                                                                  >  2:377144/1‑149 (MQ=255)
                                                           aGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGa                                                                                                                                  >  1:884249/1‑100 (MQ=255)
                                                           aGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGa                                                                                                                                  <  2:884249/100‑1 (MQ=255)
                                                              cGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAGACAGTTCAAc                                                                              <  3:55724/149‑1 (MQ=255)
                                                              cGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAAc                                                                              <  1:198627/149‑1 (MQ=255)
                                                                tttGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTc                                                                            >  2:624365/1‑149 (MQ=255)
                                                                      aTTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTAc                                                                      >  3:216354/1‑149 (MQ=255)
                                                                      aTTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTAc                                                                      >  3:784890/1‑149 (MQ=255)
                                                                      aTTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTAc                                                                      >  1:272147/1‑149 (MQ=255)
                                                                       ttCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTa                                                                       <  4:10704/147‑1 (MQ=255)
                                                                       ttCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTa                                                                       >  3:10704/1‑147 (MQ=255)
                                                                       ttCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACt                                                                     <  3:199201/149‑1 (MQ=255)
                                                                          aGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGAt                                                                  >  3:155577/1‑149 (MQ=255)
                                                                           gggACTGCCGATGTCGGTGCAGCCGACACAGTTGGTGACCGAAACCCTGAACGAGCGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGAtt                                                                 <  4:376874/149‑1 (MQ=255)
                                                                           gggACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGAtt                                                                 <  3:858305/149‑1 (MQ=255)
                                                                            ggACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTa                                                                >  2:175882/1‑149 (MQ=255)
                                                                                    gATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTACCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGc                                                        <  4:757001/149‑1 (MQ=255)
                                                                                    gATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTACCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGc                                                        <  4:852826/149‑1 (MQ=255)
                                                                                         cGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATg                                                   <  2:191878/149‑1 (MQ=255)
                                                                                             gCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGATTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCg                                               <  1:377144/149‑1 (MQ=255)
                                                                                                gCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATc                                            >  1:223056/1‑149 (MQ=255)
                                                                                                  cGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTg                                          >  3:312726/1‑149 (MQ=255)
                                                                                                   gACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTgg                                         <  1:110197/149‑1 (MQ=255)
                                                                                                   gACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTgg                                         <  4:312726/149‑1 (MQ=255)
                                                                                                   gACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTgg                                         <  2:660678/149‑1 (MQ=255)
                                                                                                     cacaGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGc                                       <  4:216354/149‑1 (MQ=255)
                                                                                                      acaGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGca                                      >  3:29194/1‑149 (MQ=255)
                                                                                                               gACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATcgcgc                             <  4:384583/149‑1 (MQ=255)
                                                                                                                   gAAACCCTGAACGAACGCCAGGAGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGACGTTCGATCTGGGCATCATCGCGcgtcg                         >  2:47649/1‑149 (MQ=255)
                                                                                                                      aCCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTcc                                                                           <  3:759825/96‑1 (MQ=255)
                                                                                                                      aCCCTGAACAAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTcc                                                                           >  4:759825/1‑96 (MQ=255)
                                                                                                                        ccTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTgc                    >  2:360427/1‑149 (MQ=255)
                                                                                                                         cTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTgcg                   >  1:485582/1‑149 (MQ=255)
                                                                                                                            aacgaacgCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTGCGCTc                <  2:471109/149‑1 (MQ=255)
                                                                                                                               gaacgCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTGCGCTCCGc             >  3:165475/1‑149 (MQ=255)
                                                                                                                                  cgCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTGCGCTCCGCATc          <  4:192964/149‑1 (MQ=255)
                                                                                                                                      aGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTGCGCTCCGCATCTGAc      <  4:13882/149‑1 (MQ=255)
                                                                                                                                          gcgcgTGCTGCCGCTTAACGAGCTGAAAGACCAGCTCGAGGCGATGGAAGGCGTGCAGTTTAACCAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCCTCACGCGTCGTCTGCGCTCCGCATCTGACCGTa  <  4:387064/149‑1 (MQ=255)
                                                                                                                                          gcgcgTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTGCGCTCCGCATCTGACCGTa  <  3:447954/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTGCCGATGTCGGTGCAGCCGACACAGCTGGTGACCGAAACCCTGAACGAACGCCAGGCGCGCGTGCTGCCGCTTAACGAGCTGAAAGACAAGCTCGAGGCGATGGAAGGCGTGCAGTTTAAACAGTTCAACTCCATTACTGATTACCACTCGCTGATGTTCGATCTGGGCATCATCGCGCGTCGTCTGCGCTCCGCATCTGACCGTA  >  NC_000913/976622‑976908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑