Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,082,935 C→T 6.1% P279L (CCG→CTG)  PP_3592 → RpiR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,082,9350CT 6.1% 190.4 / 2.6 66P279L (CCG→CTG) PP_3592RpiR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

GGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTG  >  NC_002947/4082786‑4083084
                                                                                                                                                     |                                                                                                                                                     
ggTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGAcc                                                                                                                                                       <  1:1301355/150‑1 (MQ=255)
   gACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACAATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGa                                                                                                                                                    <  2:1029231/150‑1 (MQ=255)
      cGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGt                                                                                                                                                 >  2:43144/1‑150 (MQ=255)
        acacCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTg                                                                                                                                               <  1:867771/150‑1 (MQ=255)
          acCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCACACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGCTGGCAAACGAGGTCTGCAGGAAGCTCGGACCGGATGTTGaa                                                                                                                                             >  2:510105/1‑150 (MQ=255)
            cTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGAGAATTCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGACGCTCGCACCGTATGTTGaaaa                                                                                                                                           >  2:59625/1‑150 (MQ=255)
             ttCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGaaaaa                                                                                                                                          >  2:1156599/1‑150 (MQ=255)
                        ggCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGc                                                                                                                               <  1:1458994/150‑1 (MQ=255)
                           ggATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATACGCCTCGACGAGGCTGTGGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAAc                                                                                                                            <  1:441581/150‑1 (MQ=255)
                               gCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGcc                                                                                                                        >  1:378026/1‑150 (MQ=255)
                                ccAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACTGGATGTTGAAAAACGCTTGGAAGCAACTGCCg                                                                                                                       >  1:1271866/1‑150 (MQ=255)
                                ccAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACTGGATGTTGAAAAACGCTTGGAAGCAACTGCCg                                                                                                                       >  1:1271848/1‑150 (MQ=255)
                                 cAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGAAAGCAACTgaccc                                                                                                                      >  1:797710/1‑146 (MQ=255)
                                 cAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGAATGTTGAAACACGCTTGGAAGCAACTGCCGc                                                                                                                      >  2:831649/1‑150 (MQ=255)
                                    tGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGctct                                                                                                                   >  2:1882629/1‑150 (MQ=255)
                                    tGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAAAGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAACCGCTTGGAAGCAACTGCCGctct                                                                                                                   <  1:831649/150‑1 (MQ=255)
                                      ccGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTAc                                                                                                                 <  1:982387/150‑1 (MQ=255)
                                       cgacgaGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACa                                                                                                                <  1:831670/150‑1 (MQ=255)
                                        gacgaGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAACACGCCTCGACTGTGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACAt                                                                                                               >  1:1408104/1‑150 (MQ=255)
                                          cgaGGTGTTCCGACTCGCCATCGCGTCCCGACTCTTCTTCGAATCCGTCCCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATaa                                                                                                             <  2:990581/150‑1 (MQ=255)
                                           gaGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGTGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTGGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAAc                                                                                                            <  1:1792536/150‑1 (MQ=255)
                                              gtgtTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGaa                                                                                                         >  2:734214/1‑150 (MQ=255)
                                               tgtTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAAt                                                                                                        <  1:945433/150‑1 (MQ=255)
                                                    cGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATtcgtc                                                                                                   <  2:1646497/150‑1 (MQ=255)
                                                    cGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATtcgtc                                                                                                   <  1:1882629/150‑1 (MQ=255)
                                                     gCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATtcgtcg                                                                                                  <  2:330527/150‑1 (MQ=255)
                                                      cATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCgg                                                                                                 <  1:43144/150‑1 (MQ=255)
                                                          cccACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGATGCACCTGCCGCTCTACATAACGAATTCGTCGGCTac                                                                                             >  1:408359/1‑150 (MQ=255)
                                                           ccACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATtcgtc                                                                                                   <  2:695163/143‑1 (MQ=255)
                                                           ccACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATtcgtc                                                                                                   >  1:695163/1‑143 (MQ=255)
                                                           ccACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTaca                                                                                            >  2:1490616/1‑150 (MQ=255)
                                                                   ttCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCTGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTgg                                                                                    >  2:1522270/1‑150 (MQ=255)
                                                                       aCCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATtcgtc                                                                                                   >  2:385244/1‑131 (MQ=255)
                                                                       aCCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATtcgtc                                                                                                   <  1:385244/131‑1 (MQ=255)
                                                                        ccTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGc                                                                                                                      <  1:1363786/111‑1 (MQ=255)
                                                                        ccTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGc                                                                                                                      >  2:1363786/1‑111 (MQ=255)
                                                                          tGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACgg                                                                             <  2:782138/150‑1 (MQ=255)
                                                                           gTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACggg                                                                            <  1:734214/150‑1 (MQ=255)
                                                                                  gAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAaacaa                                                                     <  1:1156599/150‑1 (MQ=255)
                                                                                   aaTCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAaacaac                                                                    <  1:1908557/150‑1 (MQ=255)
                                                                                    aTCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAacaaca                                                                   <  1:583624/150‑1 (MQ=255)
                                                                                        aCCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTacac                                                                                           >  2:1564951/1‑122 (MQ=255)
                                                                                        aCCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTacac                                                                                           <  1:1564951/122‑1 (MQ=255)
                                                                                        aCCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGc                                                                           >  2:619553/1‑138 (MQ=255)
                                                                                        aCCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGc                                                                           <  1:619553/138‑1 (MQ=255)
                                                                                                          tGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAag                                                           <  1:1737767/136‑1 (MQ=255)
                                                                                                          tGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAag                                                           >  2:1737767/1‑136 (MQ=255)
                                                                                                             gTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCga                                          >  2:546691/1‑150 (MQ=255)
                                                                                                              tGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAg                                                >  1:1681480/1‑143 (MQ=255)
                                                                                                              tGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAg                                                <  2:1681480/143‑1 (MQ=255)
                                                                                                               gCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACTGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCgaga                                        <  2:1271848/150‑1 (MQ=255)
                                                                                                               gCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACTGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCgaga                                        <  2:1271866/150‑1 (MQ=255)
                                                                                                                 aTCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATg                                      <  2:378026/150‑1 (MQ=255)
                                                                                                                    tGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAAc                                   >  1:1030896/1‑150 (MQ=255)
                                                                                                                      aTGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATGGCAAGAGTGCAATGAGGTGCGAGATGAACaa                                 <  1:510105/150‑1 (MQ=255)
                                                                                                                      aTGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACaa                                 >  1:950297/1‑150 (MQ=255)
                                                                                                                            aaCGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCAt                           >  1:1175347/1‑150 (MQ=255)
                                                                                                                               gAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGc                                                              >  2:576500/1‑112 (MQ=255)
                                                                                                                               gAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGc                                                              <  1:576500/112‑1 (MQ=255)
                                                                                                                                     tGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCt                  >  1:334232/1‑150 (MQ=255)
                                                                                                                                        aagaagCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGc                                                                           >  1:1443379/1‑90 (MQ=255)
                                                                                                                                        aagaagCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGc                                                                           <  2:1443379/90‑1 (MQ=255)
                                                                                                                                          gaagCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAaacaa                                                                     >  2:1197967/1‑94 (MQ=255)
                                                                                                                                          gaagCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTg                                             >  1:807241/1‑118 (MQ=255)
                                                                                                                                          gaagCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTg                                             <  2:807241/118‑1 (MQ=255)
                                                                                                                                          caagCTCGTACCGGATGTTGCAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAaacaa                                                                     <  1:1197967/93‑1 (MQ=255)
                                                                                                                                                 ggACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACAACTCGTGGTGGACGGGCAAACAACAATAGCACGAGTGCAATGAGCTGCGAGATGAACAAGTCCATGCCCTTGCTAGGCGCGTgcgc      >  2:580784/1‑150 (MQ=255)
                                                                                                                                                     cGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACg                                                                              <  1:134170/74‑1 (MQ=255)
                                                                                                                                                     cGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACg                                                                              >  2:134170/1‑74 (MQ=255)
                                                                                                                                                     cGGATGTTGAAAAACACTTGGCAGCAACTGCAGCACTAACTGACGAAGTCGCCGGCTACACCTCGTGCTCGACGGGCAAACAACACTAGCAATAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCctg  <  2:408359/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTG  >  NC_002947/4082786‑4083084

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: