Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,161,034 G→C 5.4% intergenic (+16/‑193) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,161,0340GC 5.4% 557.6 / 16.0 166intergenic (+16/‑193)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/118);  new base C (3/6);  total (42/124)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 173 total aligned reads displayed.

GGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGG  >  NC_002947/4160886‑4161183
                                                                                                                                                    |                                                                                                                                                      
gggCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGC                                                                                                                                                       <  2:1666282‑M1/150‑1 (MQ=255)
 ggCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCT                                                                                                                                                      <  2:4486297‑M1/150‑1 (MQ=255)
   cATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGC                                                                                                                                                    <  2:3982822‑M1/150‑1 (MQ=255)
      ggAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCggggggggg                                                                                                                                                 >  2:2361747/1‑143 (MQ=255)
      ggAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCgggggcgtg                                                                                                                                                 >  1:4346380/1‑143 (MQ=255)
       gagaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCG                                                                                                                                                <  1:3041798‑M1/150‑1 (MQ=255)
        agaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCCACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGC                                                                                                                                               <  2:5267999‑M1/150‑1 (MQ=255)
        agaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGC                                                                                                                                               <  1:2237450‑M1/150‑1 (MQ=255)
          aCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCTGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGG                                                                                                                                            <  2:4680320‑M1/150‑1 (MQ=255)
          aCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGG                                                                                                                                            <  2:475034‑M1/150‑1 (MQ=255)
          aCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGG                                                                                                                                            <  2:1878799‑M1/150‑1 (MQ=255)
          aCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGG                                                                                                                                            >  1:2100004‑M1/1‑150 (MQ=255)
            ccGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCGGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCG                                                                                                                                             >  2:568044‑M1/1‑147 (MQ=255)
            ccGCTTCAGCCAGATTCACCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGT                                                                                                                                          >  1:4043238‑M1/1‑150 (MQ=255)
              gCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaa                                                                                                                                        >  1:1252142/1‑150 (MQ=255)
               cTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaa                                                                                                                                       <  1:4485207/150‑1 (MQ=255)
                   aGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCg                                                                                                                                   <  2:1252142/150‑1 (MQ=255)
                   aGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGG                                                                                                                                            >  2:2049937‑M1/1‑141 (MQ=255)
                    gCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGc                                                                                                                                  <  2:1463076/150‑1 (MQ=255)
                    gCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTgaacccgc                                                                                                                                  <  1:4696989‑M1/150‑9 (MQ=255)
                    gCCAGAGTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGc                                                                                                                                  <  2:3833255/150‑1 (MQ=255)
                     ccAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGG                                                                                                                                            >  1:2209229‑M1/1‑139 (MQ=255)
                      cAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTc                                                                                                                                <  1:1768377/150‑1 (MQ=255)
                       aGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGGTGGTTGGGGGGTttaaccgctcc                                                                                                                               >  2:3748392‑M1/1‑139 (MQ=255)
                         aTTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTa                                                                                                                             <  1:2956682/150‑1 (MQ=255)
                          ttCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGGGTCGCGGGGTAAAccccgcccac                                                                                                                             >  1:1745915/1‑143 (MQ=255)
                             gCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                         >  2:193620/1‑150 (MQ=255)
                              ccAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggt                                                                                                                        <  2:2984688/150‑1 (MQ=255)
                                aGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                      <  2:4834379/150‑1 (MQ=255)
                                aGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                      <  1:60559/150‑1 (MQ=255)
                                    tGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGtgc                                                                                                                  <  1:1618591/150‑1 (MQ=255)
                                     gggCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGtgct                                                                                                                 <  1:5121555/150‑1 (MQ=255)
                                        ccACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGcc                                                                                                              <  1:3946255/150‑1 (MQ=255)
                                        ccACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGcc                                                                                                              <  1:4251492/150‑1 (MQ=255)
                                          aCGGGCGCTTGCCACCGGTGCGATCCTCAACCTTTGGCTGGTTGGCGGGGGGCTGGCGTTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTATACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                            <  2:1552180/150‑1 (MQ=255)
                                           cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTg                                                                                                           <  1:698137/150‑1 (MQ=255)
                                            gggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTgaacccgctcctacaggggccttgccgtcttc                                                                                                          >  2:1503225‑M1/1‑118 (MQ=255)
                                            gggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGATGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGa                                                                                                          <  1:2964658/150‑1 (MQ=255)
                                             ggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                         <  2:1416775/150‑1 (MQ=255)
                                             ggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGATGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                         <  2:370871/150‑1 (MQ=255)
                                              gcgcTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgctcctacaggggccttgccgtcttcag                                                                                                        <  2:3631208‑M1/150‑35 (MQ=255)
                                                 cTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCgg                                                                                                                                                        <  1:2555603/100‑1 (MQ=255)
                                                  tGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGC                                                                                                                                               >  2:773999‑M1/1‑108 (MQ=255)
                                                  tGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGC                                                                                                                                               <  1:773999‑M1/108‑1 (MQ=255)
                                                     ccACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                 <  2:62603/150‑1 (MQ=255)
                                                      cACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGCTGCCTGCTTCGCGGGTgaacccgctcctacaggggccttgccgtcttcagggcagcac                                                                                                >  1:2844539‑M1/1‑108 (MQ=255)
                                                       aCCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGCTGCCTGCTTCGCGGGTgaacccgctcctacaggggccttgccgtcttcagggcagcacc                                                                                               <  2:2844539‑M1/150‑44 (MQ=255)
                                                          ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgatcctacaggggccttgccgtcttcagggcagcacccac                                                                                            <  2:2100004‑M1/150‑47 (MQ=255)
                                                          ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                            <  2:3173127/150‑1 (MQ=255)
                                                          ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                            <  2:438960/150‑1 (MQ=255)
                                                          ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                            <  1:3205712/150‑1 (MQ=255)
                                                          ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                            <  1:4522838/150‑1 (MQ=255)
                                                           gTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGcccc                                                                                           <  1:5240640/150‑1 (MQ=255)
                                                           gTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGAGTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGcccc                                                                                           <  1:4640194/150‑1 (MQ=255)
                                                            tGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgctcctacaggggccttgccgtcttcagggcagcacccacct                                                                                          <  1:42111‑M1/150‑49 (MQ=255)
                                                               gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgta                                                                                       <  2:674856‑M1/150‑52 (MQ=255)
                                                               gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                       <  1:2836320/150‑1 (MQ=255)
                                                               gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTgaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgta                                                                                       <  1:1067169‑M1/150‑52 (MQ=255)
                                                                aTCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCAATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtag                                                                                      <  2:2723033‑M1/150‑53 (MQ=255)
                                                                  ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTTGAGACCGAGCTGCCTGCTTCGCGGGTgaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtagga                                                                                    <  1:3604854‑M1/150‑55 (MQ=255)
                                                                  ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtagga                                                                                    <  1:193620‑M1/150‑55 (MQ=255)
                                                                  ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtagga                                                                                    <  1:539158‑M1/150‑55 (MQ=255)
                                                                  ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                      >  1:2353181/1‑116 (MQ=255)
                                                                  ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                      <  2:2353181/116‑1 (MQ=255)
                                                                  ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATGACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa                                                                                    <  2:666487/150‑1 (MQ=255)
                                                                      caCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                <  2:4998516/149‑1 (MQ=255)
                                                                      aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAAccc                                                                                                                                    <  2:2540108/98‑1 (MQ=255)
                                                                      aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                <  1:2080072/150‑1 (MQ=255)
                                                                      aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                <  1:1333541/150‑1 (MQ=255)
                                                                      aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                <  2:4060348/150‑1 (MQ=255)
                                                                      aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                <  2:4068543/150‑1 (MQ=255)
                                                                       aCCTGTGGCTGGTTGGCGGGGGGCTGGGGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGTCTGGAGACCGTGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg                                                                               <  1:2475191/150‑1 (MQ=255)
                                                                           gtgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCa                                                                           >  1:323867‑M2/1‑149 (MQ=255)
                                                                            tgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                       <  2:509533/137‑1 (MQ=255)
                                                                            tgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                       >  1:509533/1‑137 (MQ=255)
                                                                              gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTgaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggtttaccc                                                                        <  1:118517‑M1/150‑67 (MQ=255)
                                                                                tGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTaaacccgctcctacaggggccttgccgtcttcagggcagcacccacctgtagg                                                                                     >  2:3620719‑M1/1‑82 (MQ=255)
                                                                                   ttGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTGAacccgctcc                                                                   >  2:2868123‑M2/1‑141 (MQ=255)
                                                                                   ttGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCacccgcgaa                                                                   >  2:1189408‑M2/1‑141 (MQ=255)
                                                                                                 gCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc                                                                                            <  2:544439/111‑1 (MQ=255)
                                                                                                  cGATGACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtc                                                    <  2:1112161/150‑1 (MQ=255)
                                                                                                   gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTAGCGGGTAAACCCCCTCCAACAGGTGGGTGCTGCCCAGAAGACGGCAAGGCCCCAGTAGGTGCGGGTTCAACCGCCAAGCAGGCAGCGCGGtgt                                                   >  2:379287/1‑148 (MQ=255)
                                                                                                     ttACCCTGGGGGTGCAGACGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCACGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa                                                 <  1:2361747/150‑1 (MQ=255)
                                                                                                         ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                <  1:918666/115‑1 (MQ=255)
                                                                                                         ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGGTAACCCGCTCCTACACGTGGGTGCGGCCCGGAAGAAGGGCAGGCCCCGGTAGGATCGGGGTCAACCCCGCAGCAGGCAGCGCGGTCTCAAGGt                                             >  2:5212635/1‑150 (MQ=255)
                                                                                                             ggggTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttttt                                         <  2:4560141/150‑1 (MQ=255)
                                                                                                             ggggTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttttt                                         <  1:3094478/150‑1 (MQ=255)
                                                                                                               ggTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGACGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTca                                       <  2:5032541/150‑1 (MQ=255)
                                                                                                                gTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTcac                                      <  2:4204793/150‑1 (MQ=255)
                                                                                                                 tGCAGACGGTGTGGGGCTGACTGGAGACCGAGTCGGCTGCTTCGGGGGTAAACCCGCTCCTACTGGTGGGTGCGGCCCTGAAGCCGGAAAGGCCCCTGTCGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTcaca                                     <  1:560596/150‑1 (MQ=255)
                                                                                                                      atgCTGTTGGGCTGACTGGCGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGCCGGCACGGCCCCTTTTGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTCACAGccccc                               <  1:55398/148‑1 (MQ=38)
                                                                                                                      aaGGTGTTGGGCTGCCTGGAGCCCGAGTCGGCTGCTTCGGGGTTAAACCCGCTCCTACCGGTGGGTGCGGCCCTGAAGACGGCCAGGCCCCTGTAGGCGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcccc                                <  1:1086384/150‑1 (MQ=255)
                                                                                                                                 cTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGTTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCg                     <  1:3677821/150‑1 (MQ=255)
                                                                                                                                  tGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGt                    <  2:4183850/150‑1 (MQ=255)
                                                                                                                                   gACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTc                   <  2:4346380/150‑1 (MQ=255)
                                                                                                                                     cTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCca                 <  1:4276405/150‑1 (MQ=255)
                                                                                                                                           aCCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                 <  1:2912224/64‑1 (MQ=255)
                                                                                                                                             cGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc         <  2:5346445/150‑1 (MQ=255)
                                                                                                                                                gTCGGCTGCTTCGCGGTTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGt      <  1:1595593/150‑1 (MQ=255)
                                                                                                                                                    gctgctTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg  <  1:379287/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGG  >  NC_002947/4160886‑4161183

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