Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,081,166 C→T 5.2% A464A (GCG→GCA)  PP_2691 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,081,1660CT 5.2% 294.3 / 4.2 96A464A (GCG→GCA) PP_2691hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (50/41);  new base T (2/3);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ATGTCCCGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTC  >  NC_002947/3081017‑3081315
                                                                                                                                                     |                                                                                                                                                     
aTGTCCCGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCagcagc                                                                                                                                                       <  1:933201/150‑1 (MQ=255)
   tCCCGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAgcgcg                                                                                                                                                    >  1:960473/1‑150 (MQ=255)
     ccGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTTTTTTTTTTTTTTGTCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCt                                                                                                                                                  >  2:962733/1‑150 (MQ=255)
       gggCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTAt                                                                                                                                                <  1:475571/150‑1 (MQ=255)
        ggCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATg                                                                                                                                               <  1:854216/150‑1 (MQ=255)
              ttcagtgctgccgatgcggagcacttcttctggtaccgctcagctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagacatgcagatcgataacgcccgcaataTCCAGCAGCGCGCTATGGCActg                                                                                                                                         <  2:1840364‑M1/23‑1 (MQ=255)
               agGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgc                                                                                                                                        <  1:907233/150‑1 (MQ=255)
                 agtgctgccgatgcggagcacttcttctggtaccgctcagctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActnctg                                                                                                                                      >  2:1501685‑M1/125‑150 (MQ=255)
                 agtgctgccgatgcggagcacttcttctggtaccgctcagctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtaccgcccgcaataTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      >  2:1969128‑M1/125‑150 (MQ=255)
                 agtgctgccgaggcggagcacttcttctggtaccgctccgcggaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      >  1:1840364‑M1/125‑150 (MQ=255)
                     gCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCGGCGCGCTATGGCACTGCTGGata                                                                                                                                  >  1:1459052/1‑150 (MQ=255)
                     ctgccgatgcggagcacttcttctggtaccgctcagctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCACTGCTGGata                                                                                                                                  >  1:843462‑M1/121‑150 (MQ=255)
                      tgtcgatgcggagcacttcttctggtaccgctcagctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      <  2:1161743‑M1/26‑1 (MQ=255)
                      tgtcgatgcggagcacttcttctggtaccgctcagctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      >  1:1161743‑M1/120‑145 (MQ=255)
                      cACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGatat                                                                                                                                 >  1:1827494/1‑150 (MQ=255)
                       aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      <  1:1473121/144‑1 (MQ=255)
                       aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      >  2:1473121/1‑144 (MQ=255)
                       aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATAtt                                                                                                                                >  2:914678/1‑150 (MQ=255)
                       aCGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTAATAGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATAtt                                                                                                                                >  1:124263/1‑150 (MQ=255)
                            tttCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCggg                                                                                                                           <  1:854002/150‑1 (MQ=255)
                                aCCGTAGTTCGTTTTTGCCATTCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCT‑GATATTGCGGGCGATa                                                                                                                      >  2:2094999/1‑150 (MQ=255)
                                aCCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGAt                                                                                                                       >  1:1837750/1‑150 (MQ=255)
                                   gTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAAAATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCGGCGCGCTATGGCACTGCTGGATATTGCGGGCGATAcc                                                                                                                    <  2:1459052/150‑1 (MQ=255)
                                    tAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATAccc                                                                                                                   >  2:191109/1‑150 (MQ=255)
                                          gTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGc                                                                                                             >  1:26352/1‑150 (MQ=255)
                                           tttttGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCgg                                                                                                                            >  1:726092/1‑134 (MQ=255)
                                           tttttGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCgg                                                                                                                            <  2:726092/134‑1 (MQ=255)
                                                 ccATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCtt                                                                                                      <  2:2245388/150‑1 (MQ=255)
                                                  cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGa                                                                                                                        >  2:1030944/1‑131 (MQ=255)
                                                  cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGa                                                                                                                        >  2:1030965/1‑131 (MQ=255)
                                                  cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGa                                                                                                                        <  1:1030944/131‑1 (MQ=255)
                                                  cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGa                                                                                                                        <  1:1030965/131‑1 (MQ=255)
                                                  cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCttt                                                                                                     <  2:26352/150‑1 (MQ=255)
                                                  cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCttt                                                                                                     >  1:1113528/1‑150 (MQ=255)
                                                  cATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCttt                                                                                                     >  2:1671201/1‑150 (MQ=255)
                                                   aTGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCtttt                                                                                                    >  2:325102/1‑150 (MQ=255)
                                                     gCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTc                                                                                                  <  1:357783/150‑1 (MQ=255)
                                                     gCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTc                                                                                                  <  1:359547/150‑1 (MQ=255)
                                                     gCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTc                                                                                                  >  1:2178007/1‑150 (MQ=255)
                                                     gCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTc                                                                                                  <  1:914678/150‑1 (MQ=255)
                                                        gctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActgct                                                                                                                                       >  1:1081364‑M1/86‑110 (MQ=255)
                                                        gctgaaaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActgct                                                                                                                                       <  2:1081364‑M1/25‑1 (MQ=255)
                                                          gCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCggctc                                                                                             <  1:325102/150‑1 (MQ=255)
                                                             ggcggcGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCggctcctg                                                                                          >  2:1646044/1‑150 (MQ=255)
                                                             aaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCACTGCTGGa                                                                                                                                    <  2:1848616‑M1/28‑1 (MQ=255)
                                                             aaaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCACTGCTGGa                                                                                                                                    >  1:1848616‑M1/81‑108 (MQ=255)
                                                              gcggcgTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCggctcctgg                                                                                         >  2:978678/1‑150 (MQ=255)
                                                              aaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      <  2:601180‑M1/26‑1 (MQ=255)
                                                              aaagaagagccgcgcctgggcattcgcgcagaggagccaggagccgaaaaagaaatgcagctgggtatcgcccgcaataTCCAGCAGCGCGCTATGGCActgctg                                                                                                                                      >  1:601180‑M1/80‑105 (MQ=255)
                                                                 gcgTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCggctcctggctc                                                                                      <  1:978678/150‑1 (MQ=255)
                                                                    tACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGctcctc                                                                                   >  1:5677/1‑150 (MQ=255)
                                                                     aCCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGCCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCctct                                                                                  >  1:508528/1‑150 (MQ=255)
                                                                       cgcgCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTgc                                                                                <  1:1015741/150‑1 (MQ=255)
                                                                         cgcCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTgcgc                                                                              <  1:407778/150‑1 (MQ=255)
                                                                             cTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAAt                                                                          <  2:1670476/150‑1 (MQ=255)
                                                                                       ggTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGgcgc                                                                >  2:1467811/1‑150 (MQ=255)
                                                                                        gTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGcgcg                                                               >  1:368404/1‑150 (MQ=255)
                                                                                        gTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGcgcg                                                               >  1:1966850/1‑150 (MQ=255)
                                                                                         ttGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCgg                                                              <  2:1827494/150‑1 (MQ=255)
                                                                                             gCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGctct                                                          >  2:2296556/1‑150 (MQ=255)
                                                                                              cGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtctt                                                         >  2:1162652/1‑150 (MQ=255)
                                                                                              cGATCAGGAAGTGAGCGGTCAGCCACTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtctt                                                         <  1:2501/150‑1 (MQ=255)
                                                                                                 tCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTcttctt                                                      >  1:2019602/1‑150 (MQ=255)
                                                                                                  cAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCTCGATTTCCCAGGCGCGGCTCTTCttt                                                     >  2:555138/1‑150 (MQ=255)
                                                                                                    ggAAGTGCGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGTGCGCTATGGCACTGCTGGATATTGCGCGCGATACCCAGCTGCATTTATTTTTCGGCTCCTGGCTCCTCTACGCGAATGCCCAGGCGCGGCTCTTCttttt                                                   <  2:287840/150‑1 (MQ=255)
                                                                                                    ggAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCttttt                                                                                                   <  1:2014928/102‑1 (MQ=255)
                                                                                                    ggAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCttttt                                                                                                   >  2:2014928/1‑102 (MQ=255)
                                                                                                    ggAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCCGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCttttt                                                   >  2:160620/1‑150 (MQ=255)
                                                                                                      aaGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCa                                                 >  2:1427910/1‑150 (MQ=255)
                                                                                                         tGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCt                                              >  1:581863/1‑150 (MQ=255)
                                                                                                         tGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCt                                              <  1:1162652/150‑1 (MQ=255)
                                                                                                            gCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAg                                           <  2:2178007/150‑1 (MQ=255)
                                                                                                               gTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCgg                                        >  1:2008736/1‑150 (MQ=255)
                                                                                                                tCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGt                                       >  1:1697673/1‑150 (MQ=255)
                                                                                                                tCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGt                                       <  1:962733/150‑1 (MQ=255)
                                                                                                                  aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCggctcct                                                                                           <  1:387902/96‑1 (MQ=255)
                                                                                                                  aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCggctcct                                                                                           >  2:387902/1‑96 (MQ=255)
                                                                                                                  aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttc                                                        >  1:640231/1‑131 (MQ=255)
                                                                                                                  aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttc                                                        <  2:640231/131‑1 (MQ=255)
                                                                                                                  aGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTAc                                     >  2:485184/1‑150 (MQ=255)
                                                                                                                          aGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAgaagaa                             <  2:368404/150‑1 (MQ=255)
                                                                                                                                cAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTc                       <  1:1467811/150‑1 (MQ=255)
                                                                                                                                 aGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTcc                      <  1:17520/150‑1 (MQ=255)
                                                                                                                                   ttGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTACATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGc                    >  1:1796976/1‑150 (MQ=255)
                                                                                                                                   gcccgcaataTCCAGCAGTGCCATAGCGCGCTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCa                                                 >  2:1759913‑M1/11‑121 (MQ=255)
                                                                                                                                   gcccgcaataTCCAGCAGTGCCATAGCGCGCTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCa                                                 <  1:1759913‑M1/111‑1 (MQ=255)
                                                                                                                                      cgcaataTCCAGCAGTGCCATAGCGCGCTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGc                                                             <  1:458388‑M1/99‑1 (MQ=255)
                                                                                                                                      cgcaataTCCAGCAGTGCCATAGCGCGCTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGc                                                             >  2:458388‑M1/8‑106 (MQ=255)
                                                                                                                                      cAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAgaagaa                             >  2:1639284/1‑138 (MQ=255)
                                                                                                                                      cAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAgaagaa                             <  1:1639284/138‑1 (MQ=255)
                                                                                                                                        ggTGCTCCAGCAGCGCGCTATGGCACTGCTGGACATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCgg               >  1:1559001/1‑150 (MQ=255)
                                                                                                                                         gTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttct                                                       >  2:2147490/1‑109 (MQ=255)
                                                                                                                                         gTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCtcttct                                                       <  1:2147490/109‑1 (MQ=255)
                                                                                                                                         gTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGgc              >  2:430834/1‑150 (MQ=255)
                                                                                                                                          tGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTCCCTCGCAGCGGTACCAGAAGAAGTGCTCCGCATCGgca             <  1:191109/150‑1 (MQ=255)
                                                                                                                                            cTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGgcagc           <  1:430834/150‑1 (MQ=255)
                                                                                                                                               CAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGATGCAATTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACt        <  1:1864979‑M1/150‑1 (MQ=255)
                                                                                                                                                   AGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAg    <  1:1165155‑M1/150‑1 (MQ=255)
                                                                                                                                                     cgcgcTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTc  >  2:836724/1‑150 (MQ=255)
                                                                                                                                                     cgcgcTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTc  <  2:1697673/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGTCCCGGGCCAGCAGGTTGGCACGGATTTCACCGTAGTTCGTTTTTGCCATGCCTTGCAGGCGGCGTACCGCGCCCTTGAACATCGGTTGGGCGATCAGGAAGTGAGCGGTCAGCCCCTGAGGCTGCAGTTGCAGGTGCTCCAGCAGCGCGCTATGGCACTGCTGGATATTGCGGGCGATACCCAGCTGCATTTCTTTTTCGGCTCCTGGCTCCTCTGCGCGAATGCCCAGGCGCGGCTCTTCTTTTTCAGCTGAGCGGTACCAGAAGAAGTGCTCCGCATCGGCAGCACTGAAGTC  >  NC_002947/3081017‑3081315

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