Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,971,713 T→A 5.2% Y520F (TAC→TTC) ‡ PP_5236 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,971,7130TA 5.2% 453.2 / 9.6 135Y520F (TAC→TTC) ‡PP_5236hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (45/83);  new base A (4/3);  total (49/86)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 135 total aligned reads displayed.

CGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGC  >  NC_002947/5971566‑5971861
                                                                                                                                                   |                                                                                                                                                    
cGGCGTAACAAAATCCCCAGGTGTAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTAc                                                                                                                                                    <  1:4246984/150‑1 (MQ=255)
    gTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                                <  1:2259347/150‑1 (MQ=255)
    gTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                                >  1:215226/1‑150 (MQ=255)
    gTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTgg                                                                                                                                                <  2:4265998/150‑1 (MQ=255)
     tAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgc                                                                                                                                               >  2:1594124/1‑150 (MQ=255)
     tAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgc                                                                                                                                               <  2:3827994/150‑1 (MQ=255)
      aaCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcc                                                                                                                                              <  1:3604358/150‑1 (MQ=255)
      aaCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcc                                                                                                                                              <  1:3604308/150‑1 (MQ=255)
       aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcct                                                                                                                                             <  1:1214873/150‑1 (MQ=255)
       aCAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAAGTACAACCTGgcct                                                                                                                                             <  2:4643828/150‑1 (MQ=255)
        cAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctg                                                                                                                                            >  2:4946620/1‑150 (MQ=255)
         aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                           <  1:1820639/150‑1 (MQ=255)
         aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                           <  1:2174557/150‑1 (MQ=255)
         aaaaTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGgcctgc                                                                                                                                           <  1:4077540/150‑1 (MQ=255)
            aTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGcctgcctg                                                                                                                                        <  2:4689744/150‑1 (MQ=255)
            aTCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGcctgcctg                                                                                                                                        <  2:4056803/150‑1 (MQ=255)
             tcaccAGGTGCAACCGCGCACACGCAATGGACTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGTCAGGCCAGGTTGTACTTGGCctgcctgc                                                                                                                                       >  1:3781565/4‑150 (MQ=255)
                ccAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGTGCAGGCAGGCCAAGTACAACCTGGCctgcctgcgca                                                                                                                                    <  2:2000609/150‑4 (MQ=39)
                       cAACCGCGCACACGAATTGGTCTGCCCGTTGCCCCGCGCGTCGACACGCCGGCGCCTGCCCGCCCGCAGGGCACTGCGCTCGTTCAGATCCTCAGCGCCGGCCATGCGCAGGCAGGCCTGGTTGTACTGGGCCTGCCTGCTGCCCGCCTc                                                                                                                             <  1:4394494/150‑1 (MQ=255)
                        aaCCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTcc                                                                                                                            >  2:1759623/1‑150 (MQ=255)
                        aaCCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTcc                                                                                                                            >  2:4989051/1‑150 (MQ=255)
                         aCCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCACCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCa                                                                                                                           <  2:2264747/150‑1 (MQ=255)
                            gcgcACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCg                                                                                                                        <  1:2883135/150‑1 (MQ=255)
                                acacGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACg                                                                                                                    <  2:2399421/150‑1 (MQ=255)
                                acacGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACg                                                                                                                    <  2:2399404/150‑1 (MQ=255)
                                     cATTGGTCTGCCCGTTGCCCAGCGCGTCGAAAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCTGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCa                                                                                                               <  1:3087372/150‑1 (MQ=255)
                                      aTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAg                                                                                                              <  2:2651390/150‑1 (MQ=255)
                                      aTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAg                                                                                                              >  2:2382007/1‑150 (MQ=255)
                                      aTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAg                                                                                                              >  1:921378/1‑150 (MQ=255)
                                      aTTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAg                                                                                                              >  1:559019/1‑150 (MQ=255)
                                       ttGGTCGGCCCGTTGCCCAGCGCGTCGACCCGCAGGCGCCTGACCGCCAGCAGGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCAGGCGCAGGCAGGCCTGGTTGTACTGGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGc                                                                                                             <  1:2157202/150‑1 (MQ=255)
                                         ggTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCa                                                                                                           >  2:364082/1‑150 (MQ=255)
                                         ggTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCa                                                                                                           >  2:2921381/1‑150 (MQ=255)
                                          gTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAg                                                                                                          <  1:4387630/150‑1 (MQ=255)
                                          gTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAg                                                                                                          >  1:1209044/1‑150 (MQ=255)
                                             tGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAc                                                                                                       <  1:4815206/150‑1 (MQ=255)
                                               cccGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTAccc                                                                                                     >  1:3365982/1‑150 (MQ=255)
                                                cggTTGCCCAGCGCGTCGACAAGCTGGCGCGTGACCGCCAGCTTGGCACTGCGCTCGTTCTTATCCTCATCGCGGTCCATGCGCTGGCTGGCCAGTTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCg                                                                                                    <  1:2134934/148‑1 (MQ=255)
                                                ccctggCCCACCCCGTCGACACGCCGGCGCCTGACCGCCCGCAGGGCACTGCGCTCGTTCAGCTCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCg                                                                                                    <  1:4408365/145‑1 (MQ=255)
                                                        cAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCCGGCCGGCCAGGTTGTAGTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTc                                                                                            <  1:1830480/150‑1 (MQ=255)
                                                           cgcgTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTc                                                                                         <  2:1993859/150‑1 (MQ=255)
                                                              gTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGGACCCGGCTGCCTCGTCATc                                                                                      <  2:1025936/150‑1 (MQ=255)
                                                                cGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCa                                                                                    <  1:1669730/150‑1 (MQ=255)
                                                                 gACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAg                                                                                   <  2:3390350/150‑1 (MQ=255)
                                                                 gACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAg                                                                                   >  1:1432715/1‑150 (MQ=255)
                                                                 gACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAg                                                                                   >  1:570186/1‑150 (MQ=255)
                                                                 gACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGAAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCTGCTGCCTCGTCATCCAg                                                                                   <  2:817931/150‑1 (MQ=255)
                                                                  aCAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAgg                                                                                  >  2:4968995/1‑150 (MQ=255)
                                                                     aGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCCCCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTaa                                                                               <  2:517288/150‑1 (MQ=255)
                                                                     aGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTaa                                                                               <  2:2455805/150‑1 (MQ=255)
                                                                      gccggCGCCGGACCGCCAGCATGGCACGGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCCGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                              <  1:1445429/147‑1 (MQ=255)
                                                                      gCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                              >  1:4276076/1‑150 (MQ=255)
                                                                      gCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                              <  1:3695628/150‑1 (MQ=255)
                                                                      gCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAAt                                                                              >  2:1680814/1‑150 (MQ=255)
                                                                          gcgcCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATgcgc                                                                          >  1:577290/1‑150 (MQ=255)
                                                                          gcgcCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATgcgc                                                                          >  1:4353311/1‑150 (MQ=255)
                                                                          gcgcCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATgcgc                                                                          >  1:1774569/1‑150 (MQ=255)
                                                                              cTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCggg                                                                      <  1:2377617/150‑1 (MQ=255)
                                                                                gACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGgt                                                                    <  2:1040867/150‑1 (MQ=255)
                                                                                 aCCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGgtg                                                                   <  1:3693420/150‑1 (MQ=255)
                                                                                      cAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGc                                                              <  2:1521595/150‑1 (MQ=255)
                                                                                         cATGGCAATGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAg                                                           <  2:188715/150‑1 (MQ=255)
                                                                                            ggCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgcc                                                        >  1:4619301/1‑150 (MQ=255)
                                                                                            ggCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgcc                                                        >  2:2316642/1‑150 (MQ=255)
                                                                                             gCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccg                                                       >  1:2982196/1‑150 (MQ=255)
                                                                                             gCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCATGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccg                                                       >  1:4792766/1‑150 (MQ=255)
                                                                                              cACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgc                                                      <  2:2351129/150‑1 (MQ=255)
                                                                                              cACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgc                                                      <  2:1247224/150‑1 (MQ=255)
                                                                                              cACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgc                                                      <  1:2587578/150‑1 (MQ=255)
                                                                                              cACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCGGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgc                                                      <  1:205592/150‑1 (MQ=255)
                                                                                               aCTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGgccgcc                                                     >  2:2449418/1‑150 (MQ=255)
                                                                                                  gcgcTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATg                                                  <  2:1432715/150‑1 (MQ=255)
                                                                                                      tCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGc                                              >  2:4591611/1‑150 (MQ=255)
                                                                                                      tCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGc                                              >  1:225258/1‑150 (MQ=255)
                                                                                                        gTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAg                                            >  1:4809684/1‑150 (MQ=255)
                                                                                                            aGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCa                                        >  2:3642672/1‑150 (MQ=255)
                                                                                                             gATCCTCATCGCCGTCCATGCGCAGGCAGGCCTGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCAt                                       <  1:3630112/150‑1 (MQ=255)
                                                                                                                 cTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACaa                                   <  2:4276076/150‑1 (MQ=255)
                                                                                                                 cTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACaa                                   <  2:2543634/150‑1 (MQ=255)
                                                                                                                  tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                                  >  2:3310702/1‑150 (MQ=255)
                                                                                                                  tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                                  <  1:1594124/150‑1 (MQ=255)
                                                                                                                  tCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAg                                  >  2:4392942/1‑150 (MQ=255)
                                                                                                                   cATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAgg                                 >  1:2698007/1‑150 (MQ=255)
                                                                                                                    aTCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGc                                >  2:4418015/1‑150 (MQ=255)
                                                                                                                     tCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTACTGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGcc                               >  1:3036741/1‑150 (MQ=255)
                                                                                                                            ccATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgcag                        >  1:130391/1‑150 (MQ=255)
                                                                                                                            ccATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGgcag                        <  1:1467987/150‑1 (MQ=255)
                                                                                                                                 cgcAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCAccc                   <  1:2905584/150‑1 (MQ=255)
                                                                                                                                 agcaggcagGCCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCAccc                   >  1:3208185/2‑150 (MQ=255)
                                                                                                                                  gaaggcaggcCAGCTTGTACTTGGCCTGCCTTCTGACCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCTCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCg                  >  2:4393728/3‑150 (MQ=255)
                                                                                                                                   caggcaggCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCgg                 <  2:4141991/150‑1 (MQ=255)
                                                                                                                                    aggcaggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGc                >  1:1379142/1‑150 (MQ=255)
                                                                                                                                    aggcaggcCAAGTACAACCTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGc                >  1:4403082/1‑150 (MQ=255)
                                                                                                                                      gcaggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAg              <  1:1212038/150‑1 (MQ=255)
                                                                                                                                      gcaggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAg              <  2:1229958/150‑1 (MQ=255)
                                                                                                                                      gcaggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAg              <  2:4123917/150‑1 (MQ=255)
                                                                                                                                        aggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGt            >  2:3104279/1‑150 (MQ=255)
                                                                                                                                        aggcCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGt            >  2:3104482/1‑150 (MQ=255)
                                                                                                                                           ccAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCa         <  2:4243311/150‑1 (MQ=255)
                                                                                                                                                  gTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACCAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGgc  <  1:3697323/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGCGTAACAAAATCCCCAGGTGCAACCGCGCACACGCATTGGTCTGCCCGTTGCCCAGCGCGTCGACAAGCAGGCGCCTGACCGCCAGCATGGCACTGCGCTCGTTCAGATCCTCATCGCCGTCCATGCGCAGGCAGGCCAGGTTGTACTTGGCCTGCCTGCTGCCCGCCTCCACCGCACGCTGCAGCCAGTACCCGGCTGCCTCGTCATCCAGGTAATGCGCCGGGGTGTTGGCCAGGCCGCCATGGTGCAGCTCATACAAGGCCAGGGCAGCACCCGGCAGGTTCAGTTCGGC  >  NC_002947/5971566‑5971861

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