Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,034 C→G 8.6% T44S (ACC→AGC)  yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,0340CG 8.6% 183.7 / 8.5 70T44S (ACC→AGC) yfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base C (37/26);  new base G (3/3);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

ATCAGGTGACGGACAACCATGACCGAATCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAATATGGCT  >  NC_000913/2356886‑2357179
                                                                                                                                                    |                                                                                                                                                 
atCAGGTGACGGACAACCATGACCGACTCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAAc                                                                                                                                                   <  1:672673/149‑1 (MQ=255)
   cggTGACGGACAACCATGACCGTATCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGTAAGGTTTGCAGGTTGCAAAGCtt                                                                                                                                                <  1:787604/148‑1 (MQ=255)
   aGGTGACGGACAACCATGACCGAATCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAGGTTGCAAAGCtt                                                                                                                                                <  2:358401/149‑1 (MQ=255)
   aGGTGACGGACAACCATGACCGAATCAACAACCTCCTCCCCGCATCATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAGGTTGCAAAGCtt                                                                                                                                                <  2:294013/149‑1 (MQ=255)
     gTGACGGACAACCATGACCGAATCACGAACCTCCTCCCCGCATCATGCGGAAATTAAAACCTATATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGCCCCACTTCTCAAGCTTGGCAACCTGCTAACCTTAc                                                                                                                                              >  1:694519/1‑149 (MQ=255)
      tGACGGACAACCATGTCCGAATCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCTCGGGATTTTCTCGATATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCATACCTTACg                                                                                                                                             >  1:475621/1‑149 (MQ=255)
       gACGGACAACAATGACCGAATCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCCCGGGAGTTTCTCGCAATACATTGACCAGAACCACTGCGCAATCTTTGAACCCTGCAAACCTTACGc                                                                                                                                            >  1:618629/1‑149 (MQ=255)
                         aaTCAACAACCTCCTCCCCGCATGATGCGGTATTTAACACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCCAACCTTACGCCTGGAACCCACTAGtttt                                                                                                                          >  2:810802/1‑149 (MQ=255)
                           tCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTAt                                                                                                                        >  1:899410/1‑149 (MQ=255)
                           tCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTAt                                                                                                                        >  1:903048/1‑149 (MQ=255)
                           tCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACGAGTTTTAt                                                                                                                        >  1:151028/1‑149 (MQ=255)
                               caacCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTATATGGaa                                                                                                                    >  1:862505/1‑149 (MQ=255)
                                aacCTCCTCCCCGCATGATGCGGTATTTAAAACCTATATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACGAGAACCACTGCGCAAGCTTTGCAACCTGCTAACCTTACGCCTGGCACCAACTAGTTTTATTGaaa                                                                                                                   >  2:200508/1‑149 (MQ=255)
                                 acCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGaaaa                                                                                                                  >  1:331279/1‑149 (MQ=255)
                                   ctcctcCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGaaaaaa                                                                                                                >  1:20821/1‑149 (MQ=255)
                                    tcctccCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAg                                                                                                               <  1:103072/149‑1 (MQ=255)
                                      ctccCCGCATGATGCGGTATTTAAAATCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGtt                                                                                                             <  2:86656/149‑1 (MQ=255)
                                      ctccCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGCATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGACCCCACTAGTTTTATTGAAAAAAGtt                                                                                                             >  2:250917/1‑149 (MQ=255)
                                              atgatgCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCt                                                                                                     <  1:988796/149‑1 (MQ=255)
                                                gatgCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTta                                                                                                   <  2:357305/149‑1 (MQ=255)
                                                gatgAGGTATTTAAAACCTTTATGTTCACACCCGTAACCGCACGGGATTTTCTCGAAAAACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTTCAATCCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTta                                                                                                   <  1:238175/149‑1 (MQ=255)
                                                       tATTTAAAACCTTTATGTTCACACCCGAAACCGCCCGGGATTTTCTCGAAATACATTTACCCGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCGGGAACCCACTCGTTTTATTGCAAACAGTTTACGCGCTTACTACTCg                                                                                            <  2:74404/149‑1 (MQ=255)
                                                          ttAAAAGATTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGAt                                                                                         <  1:329220/149‑1 (MQ=255)
                                                           tAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGCAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATg                                                                                        <  1:542655/149‑1 (MQ=255)
                                                             aaaCCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGtt                                                                                      >  2:13841/1‑149 (MQ=255)
                                                                 cTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTtgt                                                                                  >  1:151346/1‑149 (MQ=255)
                                                                  tttATGTTCACACCCGAAACCGCACGTGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGAGCTTACTACTCGGATGTTTtgtg                                                                                 <  2:57036/149‑1 (MQ=255)
                                                                       gTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCg                                                                            >  1:606391/1‑149 (MQ=255)
                                                                             acCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGACCCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTGTTGTGGTCCGTGGaaa                                                                      >  2:432689/1‑149 (MQ=255)
                                                                             acCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTAATCGGATGTTTTGTGGTCCGTGGaaa                                                                      >  1:659291/1‑149 (MQ=255)
                                                                               ccGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAAcc                                                                    >  1:580184/1‑149 (MQ=255)
                                                                               ccGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAAcc                                                                    >  2:488266/1‑149 (MQ=255)
                                                                               ccGAAACCGCACGGGATTTTCTCGAAATACATTAACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTAATTGAAAAACTTTTAAGCGCTAACTACACGGATGTTTTGTGGTCCGTGGAAAcc                                                                    >  1:945411/1‑149 (MQ=255)
                                                                                      cGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACg                                                             >  1:889670/1‑149 (MQ=255)
                                                                                      cGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACg                                                             <  2:331178/149‑1 (MQ=255)
                                                                                        cACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGt                                                           <  1:72383/149‑1 (MQ=255)
                                                                                           gggATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGAc                                                        <  2:919410/149‑1 (MQ=255)
                                                                                              aTTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGc                                                     <  2:296992/149‑1 (MQ=255)
                                                                                               ttttCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCt                                                    <  1:520039/149‑1 (MQ=255)
                                                                                                tttCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCta                                                   >  1:200094/1‑149 (MQ=255)
                                                                                                          tACATTTACCAGAACCACTGCGCAAGCTTTGCACACTGCAAAACTTACGCCTGGAACCCACTATGTTAATTGAAAAAAGTTTACTCGCTTACTACTCGTATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGTCTATATCTACTGc                                         <  2:618629/149‑1 (MQ=255)
                                                                                                             aTTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTg                                      >  1:289849/1‑149 (MQ=255)
                                                                                                              tttACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGAGGGTGACGGCTATATCTACTGAGTGa                                     >  1:336935/1‑149 (MQ=255)
                                                                                                               ttACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGAt                                    <  2:899410/149‑1 (MQ=255)
                                                                                                               ttACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCATGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGAt                                    <  2:903048/149‑1 (MQ=255)
                                                                                                                   cAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGaa                                >  2:500596/1‑149 (MQ=255)
                                                                                                                   cAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGaa                                <  1:409511/149‑1 (MQ=255)
                                                                                                                    aGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAAc                               >  1:52047/1‑149 (MQ=255)
                                                                                                                         cACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCaa                          >  2:7343/1‑149 (MQ=255)
                                                                                                                           atGCGCAATCTTTGCAACCTGCAAACCTTACGCCTGGTACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTAGTCGGATGTTTTGTGGTCCTTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAAGATCAAAg                        <  1:594326/148‑1 (MQ=255)
                                                                                                                               gtaaGGTTTGCAGGTTGCAAAGCTTGCGCCTGGAACCCACTTGATTTATTGAAAAAAGTTTACGCGCTTACTACTCGGAAGTTTTGTGGTCCGTGGAAACCAGCCACGGTGCCGGCTATATCTACTGCGTGATTGACCATCAAAGctcg                    >  2:318224/3‑148 (MQ=255)
                                                                                                                               gtaaGGTTTGCAGGTTGCAAAGCTTACGCCTGGATCCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGTCTATATCTACTGCGTGATTGAACATCAAAGctct                    >  1:540245/3‑149 (MQ=255)
                                                                                                                               gcAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATCTGTTGTGGTCCGTGGAAACCAGCGACGGAGACGGCTATATCTACTGCCTGATTTAACAACAAAGctct                    >  2:385423/1‑149 (MQ=255)
                                                                                                                                 aaGGTTTGCAGGTTGCAAAGCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGc                  >  2:639968/1‑149 (MQ=255)
                                                                                                                                 aaGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGc                  >  2:373213/1‑149 (MQ=255)
                                                                                                                                 aaGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGc                  >  1:638388/1‑149 (MQ=255)
                                                                                                                                  aGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCa                 <  1:242488/149‑1 (MQ=255)
                                                                                                                                  aGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCa                 <  2:889670/149‑1 (MQ=255)
                                                                                                                                    cTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGa               >  1:16257/1‑149 (MQ=255)
                                                                                                                                       tGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGaaaa            <  2:151346/149‑1 (MQ=255)
                                                                                                                                       tGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGaaaa            <  1:736469/149‑1 (MQ=255)
                                                                                                                                         cAACCTGCAAACCTTACGCCTGGACCCCACTAGTTTTAGTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGa          >  2:532/1‑149 (MQ=255)
                                                                                                                                          aaCCTGCAAACCTTACGCCTGGAACCCACTTGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTATTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGaa         <  2:973308/149‑1 (MQ=255)
                                                                                                                                           aCCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAat        <  2:221596/149‑1 (MQ=255)
                                                                                                                                             cTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGGGACGGCTATATCTACGGGGGGATGGAACATCAAAGCTCTGCAGAGAAGAatat      >  2:969685/1‑149 (MQ=255)
                                                                                                                                             cTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGCGGAACCCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAatat      >  2:989768/1‑149 (MQ=255)
                                                                                                                                               gCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAATATgg    <  1:514635/149‑1 (MQ=255)
                                                                                                                                               gCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAATATgg    >  2:238179/1‑149 (MQ=255)
                                                                                                                                               gCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAATATgg    <  2:475621/149‑1 (MQ=255)
                                                                                                                                               gCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGGGGAACCCCGCGACGGTGACCGCTTTACCTACTGCGTGATTGAACATCAAAGCTCTGCGGAAAAGAATATgg    >  2:492635/1‑149 (MQ=255)
                                                                                                                                                 aaaCCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAATATGGCt  >  2:502403/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCAGGTGACGGACAACCATGACCGAATCAACAACCTCCTCCCCGCATGATGCGGTATTTAAAACCTTTATGTTCACACCCGAAACCGCACGGGATTTTCTCGAAATACATTTACCAGAACCACTGCGCAAGCTTTGCAACCTGCAAACCTTACGCCTGGAACCCACTAGTTTTATTGAAAAAAGTTTACGCGCTTACTACTCGGATGTTTTGTGGTCCGTGGAAACCAGCGACGGTGACGGCTATATCTACTGCGTGATTGAACATCAAAGCTCTGCAGAAAAGAATATGGCT  >  NC_000913/2356886‑2357179

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