Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,155 C→A T11N (ACC→AAC)  fabR → transcriptional repressor of fabA and fabB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,1550CA97.1% 246.6 / NA 69T11N (ACC→AAC) fabRtranscriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (37/30);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTG  >  NC_000913/4161008‑4161297
                                                                                                                                                   |                                                                                                                                              
aTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAAccg                                                                                                                                             <  2:71996/151‑1 (MQ=255)
    cAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGt                                                                                                                                         >  2:79544/1‑151 (MQ=255)
       tAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCg                                                                                                                                      >  1:110405/1‑151 (MQ=255)
        aaaTATTTTTTTTTGCGTTACCGTTCTTTCACAATACTGGAGAATCCCAGTATTTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGAAAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAAAACCGCCGTTCGc                                                                                                                                     <  1:796849/151‑1 (MQ=255)
        aaaTATTTTTTATTGCTTTCCCTTTCATTCACAATACTGGAGCAACCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGAAAGAATACAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAGAAAAAAACCGCCGTTCg                                                                                                                                      <  1:151318/150‑1 (MQ=255)
          atatGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCtg                                                                                                                                   >  1:16273/1‑151 (MQ=255)
               ttttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCtggtgg                                                                                                                               >  1:52094/1‑150 (MQ=255)
               ttttATGGCGTTGCCTTTCATTCACATTACTGGAGCAACCCAGTATGTTCATTCTCTGGTAGGGTGCCAGCAGTCCTTTTGGCAGGGATTCAGACACCGTGATGGGCGTAAGAGGGCAACAAAAAGAAAAAAACCGCCGTTc                                                                                                                                       <  2:704873/142‑1 (MQ=255)
                tttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGa                                                                                                                              >  1:126667/1‑150 (MQ=255)
                    ttGCGTTACCTTTCATTCACAATACTGGAGCAATCCAGTATTTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGcc                                                                                                                          <  2:212042/150‑1 (MQ=255)
                         ttACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCAt                                                                                                                      <  2:811186/149‑1 (MQ=255)
                          tACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTTCTTTGGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                                    <  2:500139/150‑1 (MQ=255)
                           aCCGTTCATTCACAATACTGGAGCAATCCAGTATGTCCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTAAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCAttt                                                                                                                    >  1:112724/1‑149 (MQ=255)
                            ccGTTCATTCACAATACTGGAGCAACCCAGTATTTTCATTCTCTGGTATAGTGCCAGCAGTCCTTTTGGAAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTa                                                                                                                   <  1:743875/149‑1 (MQ=255)
                             cGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc                                                                                                                >  1:86064/1‑151 (MQ=255)
                                 cattcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                                                                                                            >  2:376795/1‑151 (MQ=255)
                                   ttcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                                                                                                            >  2:751482/1‑149 (MQ=255)
                                   tccacaATACTGGAGCAACCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa                                                                                                           <  1:267704/148‑1 (MQ=255)
                                    tcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTa                                                                                                           >  1:309431/1‑149 (MQ=255)
                                         aTACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGGGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                    >  2:109010/1‑151 (MQ=255)
                                         aTACGGAAGCAATCCAGTATGTTCATTCTCGGTTATAGTGCCAGAATTACTTTGGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCACAAAAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                    <  1:379573/151‑1 (MQ=255)
                                           aCTGGAGAAATCCAGTATGTTCATTCTCTGGTATAGGGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGaa                                                                                                  <  2:980096/151‑1 (MQ=255)
                                              ggAGCAATCCAGTATTTTCATTCTCTGGTATAGTGCCAGAAGTATTTTTGGCAAGGATTCAGACATCGTGATGGGCTTAAGAGCGCAACAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                                 <  1:677903/149‑1 (MQ=255)
                                              ggAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc                                                                                               <  2:463030/151‑1 (MQ=255)
                                               gagaAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGAATTATTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGGGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa                                                                                              <  1:682395/151‑1 (MQ=255)
                                                aGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGc                                                                                               >  1:945308/1‑149 (MQ=255)
                                                  cAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc                                                                                            <  2:63329/150‑1 (MQ=255)
                                                   aaTCCAGTATGTTCATTCTCTGGTATAGTGCCAGAAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCtt                                                                                          <  1:690155/151‑1 (MQ=255)
                                                    aTCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTc                                                                                         >  1:315489/1‑151 (MQ=255)
                                                      ccAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGc                                                                                       >  1:976450/1‑151 (MQ=255)
                                                      ccAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTGGGAAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCg                                                                                        <  2:915837/150‑1 (MQ=255)
                                                        aGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGAAAGGATTCAGACATCGTGATGGGCGAAAGAGCGAAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCa                                                                                     <  2:782272/151‑1 (MQ=255)
                                                            ttttCATTCTCTGGTTTAGTGCCAGAAGTACTTTTGGCAAGGATTCAGACAGCGTGATGGGCTTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGcc                                                                                  <  1:209880/148‑1 (MQ=255)
                                                            tGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCt                                                                                 >  1:146658/1‑151 (MQ=255)
                                                                      ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGt                                                                       >  2:434291/1‑151 (MQ=255)
                                                                       tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                                      >  1:1003184/1‑151 (MQ=255)
                                                                       tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                                      >  2:244503/1‑151 (MQ=255)
                                                                         gTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAACAGAAAAAAAACGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTg                                                                      >  1:67944/1‑149 (MQ=255)
                                                                           ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGt                                                                  >  2:504922/1‑151 (MQ=255)
                                                                           ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGt                                                                  >  1:765003/1‑151 (MQ=255)
                                                                              gTGCCAGCAGTCCTTTTGGAAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg                                                                 <  1:376783/149‑1 (MQ=255)
                                                                               tGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcg                                                              >  2:480696/1‑151 (MQ=255)
                                                                                 ccAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGcgcg                                                            <  1:16625/151‑1 (MQ=255)
                                                                                  cagcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGt                                                           >  1:426831/1‑151 (MQ=255)
                                                                                   agcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTg                                                          >  1:280993/1‑151 (MQ=255)
                                                                                    gcatTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTg                                                          <  1:38671/150‑1 (MQ=255)
                                                                                    gcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTg                                                          <  2:418569/150‑1 (MQ=255)
                                                                                        tACTTTTGGAAAGGATTCAGACGTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGa                                                         <  1:811743/147‑1 (MQ=255)
                                                                                            tttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGc                                                  <  2:530932/150‑1 (MQ=255)
                                                                                              tGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAtt                                               >  2:576956/1‑151 (MQ=255)
                                                                                              tGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAtt                                               >  1:884754/1‑151 (MQ=255)
                                                                                              tGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAt                                                <  2:357127/150‑1 (MQ=255)
                                                                                               ggCAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAACAAAAAAAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAt                                                <  1:285860/149‑1 (MQ=255)
                                                                                                gCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGc                                             >  1:113408/1‑151 (MQ=255)
                                                                                                                cGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTa                              <  2:424772/150‑1 (MQ=255)
                                                                                                                 gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATc                            >  1:919542/1‑151 (MQ=255)
                                                                                                                    aTGGGCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGc                         >  1:717055/1‑151 (MQ=255)
                                                                                                                        gCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAttt                     >  2:475324/1‑151 (MQ=255)
                                                                                                                        gCGTAAGAGCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAttt                     >  2:515329/1‑151 (MQ=255)
                                                                                                                             agagCGCAACAAAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgcg                >  2:582296/1‑151 (MQ=255)
                                                                                                                                   caacaaAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt            >  2:352593/1‑149 (MQ=255)
                                                                                                                                    aaaaaaaaaaaaaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGa         <  2:827524/148‑1 (MQ=255)
                                                                                                                                     acaaAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc        >  1:393620/1‑151 (MQ=255)
                                                                                                                                      caaAAAGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACg       >  1:532225/1‑151 (MQ=255)
                                                                                                                                       aaaaaGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACg       >  1:600489/1‑150 (MQ=255)
                                                                                                                                       aaaaaGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGa      <  2:765028/151‑1 (MQ=255)
                                                                                                                                        aaaaGAAAAAAACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGa      >  1:656147/1‑150 (MQ=255)
                                                                                                                                         aaaaaaaaaaACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAc    <  1:989380/151‑1 (MQ=255)
                                                                                                                                             aaaaaaACCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg  <  1:396398/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTG  >  NC_000913/4161008‑4161297

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: