Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,070 C→A 39.9% A436E (GCG→GAG)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947543,0700CA39.9% 93.2 / 88.3 97A436E (GCG→GAG) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base A (20/19);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTGGTCTGGCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGA  >  NC_002947/542923‑543216
                                                                                                                                                   |                                                                                                                                                  
tGCGTGGTCTGGCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGt                                                                                                                                                  <  2:447467/150‑1 (MQ=255)
 gCGTGGTCTGGCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGtt                                                                                                                                                 <  2:457012/150‑1 (MQ=255)
      gTCTGGCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTAGTGCGGAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAAc                                                                                                                                            >  1:142419/1‑150 (MQ=255)
         tGGCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCg                                                                                                                                         >  1:542613/1‑150 (MQ=255)
                caccaTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTTGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGAt                                                                                                                                  >  2:1181183/1‑150 (MQ=255)
                 accaTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATc                                                                                                                                 >  2:2290554/1‑150 (MQ=255)
                   caTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGa                                                                                                                               >  1:915284/1‑150 (MQ=255)
                     tCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAg                                                                                                                             <  1:170250/150‑1 (MQ=255)
                     tCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAg                                                                                                                             <  1:1208854/150‑1 (MQ=255)
                       aaGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGt                                                                                                                           >  1:1465598/1‑150 (MQ=255)
                           cTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTTGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAgg                                                                                                                       <  1:1181183/150‑1 (MQ=255)
                             gCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTAGTGATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCt                                                                                                                     >  1:2135775/1‑150 (MQ=255)
                               cAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTAt                                                                                                                   >  1:85873/1‑150 (MQ=255)
                                 agaagaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTc                                                                                                                 <  2:915284/150‑1 (MQ=255)
                                  gaagaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCa                                                                                                                <  1:115389/150‑1 (MQ=255)
                                    agaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGc                                                                                                              >  1:1371935/1‑150 (MQ=255)
                                          tCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCAcc                                                                                                        >  2:1734879/1‑150 (MQ=255)
                                               cgcgAGCTGCCAGAGGTGTGGGACGTTCTAGCTGGCGAATTAAGGCTGCACGGTGTTCTGCTCAACCGTGAACCCATCGTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCATCTGCACCCGCTg                                                                                                   <  1:522843/150‑1 (MQ=255)
                                                        ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCAc                                                                                                         <  2:54657/135‑1 (MQ=255)
                                                        ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCAc                                                                                                         >  1:54657/1‑135 (MQ=255)
                                                        ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgct                                                                                          >  2:866828/1‑150 (MQ=255)
                                                        ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTa                                                                                                                    >  2:1588517/1‑124 (MQ=255)
                                                        ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTa                                                                                                                    <  1:1588517/124‑1 (MQ=255)
                                                        ccAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgct                                                                                          >  1:764373/1‑150 (MQ=255)
                                                         cAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCt                                                                                               <  1:1605671/144‑1 (MQ=255)
                                                         cAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCt                                                                                               >  2:1605671/1‑144 (MQ=255)
                                                          agagGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGAt                                                                                                                                  >  1:875590/1‑108 (MQ=255)
                                                          agagGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGAt                                                                                                                                  <  2:875590/108‑1 (MQ=255)
                                                          agagGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                        >  1:1640517/1‑150 (MQ=255)
                                                              gtgtGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTAAACCCGTACATGGCGAAGATAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTAc                                                                                    >  1:1687192/1‑150 (MQ=255)
                                                                gtgGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACaa                                                                                  >  2:1522673/1‑150 (MQ=255)
                                                                   ggACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGt                                                                                                                                                  <  1:2123927/83‑1 (MQ=255)
                                                                   ggACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGt                                                                                                                                                  >  2:2123927/1‑83 (MQ=255)
                                                                   ggACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGc                                                                               >  2:2171701/1‑150 (MQ=255)
                                                                   ggACGTGCGCGCTTTAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGc                                                                               >  2:1203033/1‑150 (MQ=255)
                                                                    gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCa                                                                                                                >  2:1087449/1‑116 (MQ=255)
                                                                    gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCa                                                                                                                <  1:1087449/116‑1 (MQ=255)
                                                                       gTTCTCGCTGACGTGATTCGCAAACACCCCGGACTGGTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGAc                                                                           >  2:1311408/1‑150 (MQ=255)
                                                                       gTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGAc                                                                           <  2:542613/150‑1 (MQ=255)
                                                                       gTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGAc                                                                           >  2:438238/1‑150 (MQ=255)
                                                                       gTTCTCGCTGAAGTGATTCGCGAACAACCCTTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGAc                                                                           >  1:429366/1‑150 (MQ=255)
                                                                         tctcGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACtt                                                                         <  1:1522673/150‑1 (MQ=255)
                                                                             gCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGAc                                                                     >  2:1750059/1‑150 (MQ=255)
                                                                              cTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACg                                                                    <  2:1371935/150‑1 (MQ=255)
                                                                                 aaGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGa                                                                                                                               >  2:1336942/1‑88 (MQ=255)
                                                                                 aaGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGa                                                                                                                               <  1:1336942/88‑1 (MQ=255)
                                                                                 aaGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTg                                                                 <  1:1750059/150‑1 (MQ=255)
                                                                                 aaGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCt                                                                                                             <  1:1050892/106‑1 (MQ=255)
                                                                                 aaGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCt                                                                                                             >  2:1050892/1‑106 (MQ=255)
                                                                                 aaGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTg                                                                 >  1:1116677/1‑150 (MQ=255)
                                                                                     gATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCAGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCa                                                             >  2:302946/1‑150 (MQ=255)
                                                                                        tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTTACCAGAt                                                          <  1:1298624/150‑1 (MQ=255)
                                                                                          gcgAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATgg                                                        <  2:1116677/150‑1 (MQ=255)
                                                                                           cgAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                        >  1:1909654/1‑117 (MQ=255)
                                                                                           cgAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                        <  2:1909654/117‑1 (MQ=255)
                                                                                                cccccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCACGCCTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTc                                                  <  1:2273706/149‑1 (MQ=255)
                                                                                                  cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGc                                                                                                     >  1:829623/1‑97 (MQ=255)
                                                                                                  cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGc                                                                                                     <  2:829623/97‑1 (MQ=255)
                                                                                                  cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                              >  1:2205896/1‑120 (MQ=255)
                                                                                                  cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                              <  2:2205896/120‑1 (MQ=255)
                                                                                                  cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACgg                                                                   <  2:439159/131‑1 (MQ=255)
                                                                                                  cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACgg                                                                   >  1:439159/1‑131 (MQ=255)
                                                                                                   ccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTg                                                                                                   >  2:872048/1‑98 (MQ=255)
                                                                                                   ccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTg                                                                                                   <  1:872048/98‑1 (MQ=255)
                                                                                                    cGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGt                                              <  2:2037938/150‑1 (MQ=255)
                                                                                                      tACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGc                                            <  1:438238/150‑1 (MQ=255)
                                                                                                         tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGc                                         <  1:866828/150‑1 (MQ=255)
                                                                                                         tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTa                                                                                     >  1:2334473/1‑106 (MQ=255)
                                                                                                         tGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTa                                                                                     <  2:2334473/106‑1 (MQ=255)
                                                                                                           cTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTCCTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTg                                       >  2:81230/1‑150 (MQ=255)
                                                                                                            tCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGGCGGTGACCAGATGGCCGTTCACGTGCCGCTGa                                      >  2:7420/1‑150 (MQ=255)
                                                                                                            tCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGa                                      <  1:401052/150‑1 (MQ=255)
                                                                                                              aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCa                                                                                                          >  1:774256/1‑80 (MQ=255)
                                                                                                              aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCa                                                                                                          <  2:774256/80‑1 (MQ=255)
                                                                                                              aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATgg                                                        <  2:539168/130‑1 (MQ=255)
                                                                                                              aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATgg                                                        >  1:539168/1‑130 (MQ=255)
                                                                                                                     cACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAg                             <  2:1465598/150‑1 (MQ=255)
                                                                                                                      aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                        <  1:1175090/90‑1 (MQ=255)
                                                                                                                      aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                        >  2:1175090/1‑90 (MQ=255)
                                                                                                                      aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                            >  2:961281/1‑150 (MQ=255)
                                                                                                                        cGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGccc                          >  2:459426/1‑150 (MQ=255)
                                                                                                                        cGACCCTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGGAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGccc                          >  1:1876207/1‑150 (MQ=255)
                                                                                                                         gACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGGCGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCa                         <  1:7420/150‑1 (MQ=255)
                                                                                                                         gACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTc                                                  <  1:818536/125‑1 (MQ=255)
                                                                                                                         gACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTc                                                  >  2:818536/1‑125 (MQ=255)
                                                                                                                            ccTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCt                      >  2:2272671/1‑150 (MQ=255)
                                                                                                                             cTTCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTc                     >  2:1386412/1‑150 (MQ=255)
                                                                                                                              ttCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCa                         <  1:887949/145‑1 (MQ=255)
                                                                                                                              ttCACCGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCa                         >  2:887949/1‑145 (MQ=255)
                                                                                                                                  ccGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgc                <  2:85873/150‑1 (MQ=255)
                                                                                                                                  ccGTCTGGGTATCCAGGAGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgc                <  1:2182940/150‑1 (MQ=255)
                                                                                                                                      cTGGGTATCCAGGCGTTTGAACCCTTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCGCAGCCGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGt            <  1:1311408/150‑1 (MQ=255)
                                                                                                                                             tCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtga     <  1:2290554/150‑1 (MQ=255)
                                                                                                                                             tCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtga     >  1:1760697/1‑150 (MQ=255)
                                                                                                                                              ccAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgat    <  1:2272671/150‑1 (MQ=255)
                                                                                                                                                aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga  >  2:330154/1‑150 (MQ=255)
                                                                                                                                                aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCACCTGAGTCCTCGGGTCTGTCGTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga  >  2:80845/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCGTGGTCTGGCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGA  >  NC_002947/542923‑543216

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: