Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 496,099 T→A 5.4% intergenic (‑63/‑33) PP_0408 ← / → PP_0409 hypothetical protein/two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947496,0990TA 5.4% 287.2 / 6.5 93intergenic (‑63/‑33)PP_0408/PP_0409hypothetical protein/two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (45/43);  new base A (3/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

TTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCTG  >  NC_002947/495953‑496248
                                                                                                                                                  |                                                                                                                                                     
ttCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCa                                                                                                                                                    <  1:1014202/150‑1 (MQ=255)
       tCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGt                                                                                                                                             <  1:1103795/150‑1 (MQ=255)
        cggcggCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTc                                                                                                                                            <  1:1765822/150‑1 (MQ=255)
         ggcggcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCt                                                                                                                                           >  2:261174/1‑150 (MQ=255)
             gcCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGc                                                                                                                                       <  2:1401571/150‑1 (MQ=255)
              ccAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGcc                                                                                                                                      >  2:1027993/1‑150 (MQ=255)
                     gcCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAg                                                                                                                                   <  2:1718082/146‑1 (MQ=255)
                     gcCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAg                                                                                                                                   >  1:1718082/1‑146 (MQ=255)
                          gTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGcc                                                                                                                                               >  2:527452/1‑129 (MQ=255)
                          gTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGca                                                                                                                                               <  1:527452/129‑2 (MQ=255)
                           tGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACTGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGtat                                                                                                                         >  2:1731708/1‑150 (MQ=255)
                             gAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAg                                                                                                                       >  2:801831/1‑150 (MQ=255)
                             gAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAg                                                                                                                       >  1:154683/1‑150 (MQ=255)
                                ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCt                                                                                                                                           <  1:362311/127‑1 (MQ=255)
                                ggACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCt                                                                                                                                           >  2:362311/1‑127 (MQ=255)
                                  aCCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGACGAGCGTATAGATGaa                                                                                                                  <  1:269662/150‑1 (MQ=255)
                                   ccGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAg                                                                                                                 >  1:1686820/1‑150 (MQ=255)
                                     gAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCg                                                                                                               >  1:672741/1‑150 (MQ=255)
                                      aGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGa                                                                                                              <  1:1273982/150‑1 (MQ=255)
                                       ggCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAg                                                                                                             >  1:635078/1‑150 (MQ=255)
                                        gCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGt                                                                                                            <  1:1651063/150‑1 (MQ=255)
                                        gCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGt                                                                                                            >  2:561223/1‑150 (MQ=255)
                                               agcaCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGc                                                                                                     <  2:922437/150‑1 (MQ=255)
                                               agcaCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCAAACACCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGc                                                                                                     >  1:1788512/1‑150 (MQ=255)
                                                gcaCATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTACTAGTATTTGCCCGTGACGGCCATGCTGCATTGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGcc                                                                                                    >  2:687865/1‑150 (MQ=255)
                                                  acaTTGCCAGCGTCGTGCGAGAAAGTGTGGACATGCGTCAACCCAAGCCCAGAATGAAACCGGCAGCCAAAAAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCCAGCGTATAGATGAAGCGAGTGACCCGCCTg                                                                                                  >  2:339862/1‑150 (MQ=255)
                                                     ttGCCAGCGTCGTGCGAGAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTg                                                                                               <  2:1788512/150‑1 (MQ=255)
                                                       gCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGt                                                                                             >  2:1603171/1‑150 (MQ=255)
                                                        ccAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTgggt                                                                                            >  2:435667/1‑148 (MQ=255)
                                                        ccAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGtt                                                                                            <  2:672741/150‑1 (MQ=255)
                                                         cAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTa                                                                                           >  1:845542/1‑150 (MQ=255)
                                                         cAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTa                                                                                           >  2:619625/1‑150 (MQ=255)
                                                            cgtcgtGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCg                                                                                        <  1:1150079/150‑1 (MQ=255)
                                                            cgtcgtGCGAGAAAGTGTGGACATGAGTCAACCCAAGGCAAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCg                                                                                        <  2:1536925/150‑1 (MQ=255)
                                                              tcgtGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGgc                                                                                      >  2:1498434/1‑150 (MQ=255)
                                                                 tGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGgctgc                                                                                   >  2:921944/1‑150 (MQ=255)
                                                                  gCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGGGCGTCGCCTGTTGGTTATCGgctgct                                                                                  <  1:619625/150‑1 (MQ=255)
                                                                      tAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcacgg                                                                                                                                        >  2:1515564/1‑88 (MQ=255)
                                                                      tAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCAGCATGGCCGTcacgg                                                                                                                                        <  1:1515564/92‑5 (MQ=255)
                                                                      tAAAGTGTGGACATGAGTCAACCCAAGGCCAGAAGGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTtgtg                                                                              <  2:697560/150‑1 (MQ=255)
                                                                        aaGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGc                                                                            >  2:425244/1‑150 (MQ=255)
                                                                         aGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCt                                                                           >  1:879962/1‑150 (MQ=255)
                                                                         aGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCt                                                                           >  1:1253727/1‑150 (MQ=255)
                                                                               gACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCtt                                                                     <  1:561223/150‑1 (MQ=255)
                                                                                 cATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                             <  1:660818/92‑1 (MQ=255)
                                                                                 cATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                             >  2:660818/1‑92 (MQ=255)
                                                                                    gAGTCGACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAg                                                                                                             <  2:996429/105‑1 (MQ=255)
                                                                                    gAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAg                                                                                                             >  1:996429/1‑105 (MQ=255)
                                                                                       gcAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCGATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCt                                                             <  1:822846/149‑1 (MQ=255)
                                                                                         aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCCTGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCTCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTgg                                                           >  2:253310/1‑150 (MQ=255)
                                                                                         aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                             <  2:437979/84‑1 (MQ=255)
                                                                                         aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGc                                                                                                                             >  1:437979/1‑84 (MQ=255)
                                                                                         aaCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTTATTGCCTgg                                                           >  2:272183/1‑150 (MQ=255)
                                                                                          aCCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTggg                                                          <  2:1253727/150‑1 (MQ=255)
                                                                                           cccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGAt                                                                                                                     <  1:1857811/90‑1 (MQ=255)
                                                                                           cccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGAt                                                                                                                     >  2:1857811/1‑90 (MQ=255)
                                                                                            ccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGcc                                                        >  1:1665932/1‑150 (MQ=255)
                                                                                            ccAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGcc                                                        <  1:261174/150‑1 (MQ=255)
                                                                                               aGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAg                                                     <  1:1028135/150‑1 (MQ=255)
                                                                                                 gCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGcc                                                   >  1:467487/1‑150 (MQ=255)
                                                                                                 gCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGGCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGTCTGCTTGTAGCTGCCCTTGATTGCCGGGGCCAAccc                                                   >  2:98889/1‑150 (MQ=255)
                                                                                                     gAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGt                                                                                                        <  2:1123666/93‑1 (MQ=255)
                                                                                                     gAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGt                                                                                                        >  1:1123666/1‑93 (MQ=255)
                                                                                                      aaTGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAgg                                              >  1:1223524/1‑150 (MQ=255)
                                                                                                      aaTGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAgg                                              >  2:139618/1‑150 (MQ=255)
                                                                                                        tGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGa                                                <  2:666183/146‑1 (MQ=255)
                                                                                                        tGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGa                                                >  1:666183/1‑146 (MQ=255)
                                                                                                          aaaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCt                                          <  1:139618/150‑1 (MQ=255)
                                                                                                           aaCCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGGATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTg                                         >  2:147999/1‑150 (MQ=255)
                                                                                                                gCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATAGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGa                                                                                                              >  2:1152438/1‑76 (MQ=255)
                                                                                                                gCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATAGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGa                                                                                                              <  1:1152438/76‑1 (MQ=255)
                                                                                                                 cAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCctct                                   <  2:91272/150‑1 (MQ=255)
                                                                                                                  aGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtg                                  <  2:1686820/150‑1 (MQ=255)
                                                                                                                    ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                                <  1:921944/150‑1 (MQ=255)
                                                                                                                    ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                                <  1:706969/150‑1 (MQ=255)
                                                                                                                    ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                                <  2:1389107/150‑1 (MQ=255)
                                                                                                                    ccTAATAGTATTTGCCCGTGACGGCCATGCTGCATGCCCGTCTTTGCGAAGTCGAGCTTATAGATGAAGCGCGTGCTGCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCtgtg                                <  1:425244/150‑1 (MQ=255)
                                                                                                                     cTAATGGTATTTGCCCCTGACGGCGATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTgtgg                               <  2:970277/150‑1 (MQ=255)
                                                                                                                         tAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCGTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCt                           >  1:1378063/1‑150 (MQ=255)
                                                                                                                            tATTTGCCCGTGCCGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGa                        <  1:861630/150‑1 (MQ=255)
                                                                                                                              tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                      >  2:1171854/1‑150 (MQ=255)
                                                                                                                              tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                      <  2:1665932/150‑1 (MQ=255)
                                                                                                                              tttGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGc                      <  1:1018117/150‑1 (MQ=255)
                                                                                                                                tGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCt                    <  1:1494497/150‑1 (MQ=255)
                                                                                                                                tGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCt                    <  2:1378063/150‑1 (MQ=255)
                                                                                                                                  gcaaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCCGCTGGTCATCGGCTGCTTGTCGCTGCCCTCGATGGCCTGGGCCAATCCCGAGGCGCCGCCCCCTCTGCTGCTGGACCCTGa                  >  2:1652532/6‑150 (MQ=255)
                                                                                                                                   ccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                 <  1:801831/150‑1 (MQ=255)
                                                                                                                                   ccGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                 >  1:765814/1‑150 (MQ=255)
                                                                                                                                   caaagACGGCCATGCAGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGAc                 >  1:1908559/5‑150 (MQ=255)
                                                                                                                                      tGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcag              <  1:1498434/150‑1 (MQ=255)
                                                                                                                                          ggCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcagc          >  2:1364637/1‑150 (MQ=255)
                                                                                                                                          ggCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACcagcagc          >  1:1830920/1‑150 (MQ=255)
                                                                                                                                              aTGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTg      >  1:567797/1‑150 (MQ=255)
                                                                                                                                                  tgcATGGCCGACTTTGCCAAGACGAACGAATAAATGAAGAGAGTGTGTCGCCTGTTGGTTATCGGATGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAACCCAAGGCGCTGCCCTCCGTGGGGCTGGAACATCACCAGCAGCAGTGGCTg  >  2:856825/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTCGGCATCGGCGGCCAGGGCGCCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTCTTTGCCAAGTCGAGCGTATAGATGAAGCGAGTGCGTCGCCTGTTGGTTATCGGCTGCTTGTGGCTGCCCTTGATTGCCTGGGCCAAGCCCGAGGCGCTGCCCTCTGTGGTGCTGGAGCCTGACCAGCAGCAGTGGCTG  >  NC_002947/495953‑496248

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