Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,503,125 C→A 7.7% intergenic (‑186/‑73) tas ← / → rplM NAD(P)‑linked reductase/50S ribosomal protein L13

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,503,1250CA 7.7% 266.7 / 9.3 87intergenic (‑186/‑73)tas/rplMNAD(P)‑linked reductase/50S ribosomal protein L13
Reads supporting (aligned to +/- strand):  ref base C (42/38);  new base A (4/3);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CAGGTTGAAATTTTTCGCGCAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGACCCTG  >  NC_002947/1502977‑1503272
                                                                                                                                                    |                                                                                                                                                   
cAGGTTGTAATTTTTCGCGGAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGGTATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGGCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGccc                                                                                                                                                     <  7:46639/149‑1 (MQ=255)
cAGGTTGAAATTTTTCGCGCAATCTGGATAATTCAGCCCATCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGACCGCCGTCGAACATTGGCGTGCCCGACCGGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGccc                                                                                                                                                     <  4:332050/149‑1 (MQ=255)
   gTTGAAATTTTTCGCGCAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTg                                                                                                                                                  <  2:142839/149‑1 (MQ=255)
   gTTGAAATTTTTCGCGCAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTg                                                                                                                                                  >  2:72106/1‑149 (MQ=255)
      gAAATTTTTCGCGCAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGc                                                                                                                                                >  3:175915/1‑148 (MQ=255)
          tttttCGCGCAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCACCAAAAGAGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCt                                                                                                                                           >  6:30873/1‑149 (MQ=255)
             ttCGCGCAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCAAGGGCAAGcgcg                                                                                                                                         <  8:31035/148‑4 (MQ=255)
                   cAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTAt                                                                                                                                  <  2:252760/149‑1 (MQ=255)
                    aaTCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATc                                                                                                                                 <  5:319502/149‑1 (MQ=255)
                        tGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGt                                                                                                                             >  2:156375/1‑148 (MQ=255)
                          gATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCa                                                                                                                               <  2:119081/145‑1 (MQ=255)
                          gATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGt                                                                                                                             >  5:297193/1‑147 (MQ=255)
                             aaTTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGAt                                                                                                                        >  5:12463/1‑149 (MQ=255)
                               ttCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTc                                                                                                                      >  8:171057/1‑149 (MQ=255)
                                tCAGCCCCTCTCCGCGTGGCGGAAGTGAAGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGc                                                                                                                                                <  4:12380/122‑1 (MQ=255)
                                tCAGCCCCTCTCCGCGTGGCGGAAGTGAAGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGc                                                                                                                                                >  3:12380/1‑122 (MQ=255)
                                  aGCCCCTCTCCGCGTGGCGGAAGTGATGCACCCTCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCAGACCCGCGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGcc                                                                                                                   <  6:283208/149‑1 (MQ=255)
                                     cccTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGCCCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCCAGCTGTGCACTATCCAGTAAGATTCGCCGTc                                                                                                                >  4:56287/1‑149 (MQ=255)
                                     cccTCTCCGCCTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCAGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAACCGTCTCTGTTCGGCTCCGCGCTTGCCCTTGGCAACCTGTGCACTATCCACTAAGATTCGCCGt                                                                                                                 >  8:82848/1‑148 (MQ=255)
                                       ctctCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTa                                                                                                              >  1:26345/1‑149 (MQ=255)
                                         ctcCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTAtt                                                                                                            >  2:304759/1‑149 (MQ=255)
                                            cgcgTGGCGGAAGTGATGCTCCATGACTCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACGCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTTGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTAttttt                                                                                                         <  6:12463/149‑1 (MQ=255)
                                            cgcgTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTAttttt                                                                                                         >  7:335507/1‑149 (MQ=255)
                                            cgcgTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTAttttt                                                                                                         >  1:90150/1‑149 (MQ=255)
                                            cgcgTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTAtt                                                                                                            <  4:175915/146‑1 (MQ=255)
                                                             gCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATc                                                                                                                                 >  3:89662/1‑108 (MQ=255)
                                                             gCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATc                                                                                                                                 <  4:89662/108‑1 (MQ=255)
                                                               aCCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAgg                                                                                      >  6:4451/1‑149 (MQ=255)
                                                                 cATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCt                                                                                    >  2:62141/1‑149 (MQ=255)
                                                                 cATAACGAAGGAAGAAACCCGTCGAAAATTGGCGTGCACGACCCGTGACCCCAAAAGCGTCTGTTTTCGGCTGCGCGCTTTCACTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGTTTGGGCGCCTCAGAGGCt                                                                                    >  2:133072/1‑149 (MQ=255)
                                                                  aTCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCtt                                                                                   >  2:70611/1‑149 (MQ=255)
                                                                  aTCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCtt                                                                                   >  8:230515/1‑149 (MQ=255)
                                                                     aCGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGGCCCCAAAAGCGTCTGTTTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTa                                                                                 <  7:82848/148‑1 (MQ=255)
                                                                     aCGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGggc                                                                                                 >  7:328294/1‑132 (MQ=255)
                                                                     aCGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGggc                                                                                                 <  8:328294/132‑1 (MQ=255)
                                                                     aCGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTag                                                                                <  7:230515/149‑1 (MQ=255)
                                                                     aCGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTa                                                                                 <  1:304759/148‑1 (MQ=255)
                                                                      cGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTa                                                                                 >  2:347540/1‑147 (MQ=255)
                                                                       gCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTa                                                                                 >  6:77817/1‑146 (MQ=255)
                                                                         aGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGaa                                                                            <  8:39279/149‑1 (MQ=255)
                                                                             agaaCCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGaa                                                                        <  5:30873/149‑1 (MQ=255)
                                                                                aCCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAAc                                                                     >  6:19348/1‑149 (MQ=255)
                                                                                 ccccGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACt                                                                    >  1:110336/1‑149 (MQ=255)
                                                                                 ccccGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACt                                                                    <  5:77817/149‑1 (MQ=255)
                                                                                  cccGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTGTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACtt                                                                   >  2:188155/1‑149 (MQ=255)
                                                                                  cccGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGGCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACtt                                                                   <  3:331997/149‑1 (MQ=255)
                                                                                   ccGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACttt                                                                  >  3:317702/1‑149 (MQ=255)
                                                                                    cGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTCTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACtttt                                                                 <  1:188155/149‑1 (MQ=255)
                                                                                     gTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTa                                                                <  1:72106/149‑1 (MQ=255)
                                                                                           cATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAg                                                                                       <  7:231282/120‑1 (MQ=255)
                                                                                           cATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAg                                                                                       >  8:231282/1‑120 (MQ=255)
                                                                                               ggCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAAc                                                       <  5:19348/148‑1 (MQ=255)
                                                                                                gCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAAc                                                                     >  6:42438/1‑133 (MQ=255)
                                                                                                gCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAAc                                                                     <  5:42438/133‑1 (MQ=255)
                                                                                                  gTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGa                                                   <  8:335507/149‑1 (MQ=255)
                                                                                                    gCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGGTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGaaa                                                 <  6:86121/149‑1 (MQ=255)
                                                                                                    gCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTGCTGCTAAACCGGaaa                                                 <  2:110336/149‑1 (MQ=255)
                                                                                                     cccgacccgcGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGGAAGCTGTGCACGATCCAGTAAGATTCGCCGTCAAGTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGCATGACAACGTTTACTGCTAAACCGGaaca                                                >  1:285690/1‑147 (MQ=255)
                                                                                                     cccgacccgaGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGaaa                                                 <  1:347540/148‑1 (MQ=255)
                                                                                                     cccgacccgaGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGaaa                                                 <  1:62141/148‑1 (MQ=255)
                                                                                                     cccgacccgaGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGaaa                                                 >  6:151286/1‑148 (MQ=255)
                                                                                                     cccgacccgaGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGAGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGaaa                                                 >  6:199219/1‑148 (MQ=255)
                                                                                                      ccgacccgaGCCCCTAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACa                                               <  3:165202/149‑1 (MQ=255)
                                                                                                      ccgacccgaGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGGCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCTTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACa                                               <  7:171057/149‑1 (MQ=255)
                                                                                                      ccgacccgaGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCgg                                                    >  1:110603/1‑144 (MQ=255)
                                                                                                      ccgacccgaGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCgg                                                    <  2:110603/144‑1 (MQ=255)
                                                                                                           ccgaGCCCCCAACAGCGTGTGGGTTCGGCTGCGCGCGTGCCCTTGGCAAGCTGTGCACTATCCAGTCATATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGCTAATGAAAACTTTTACTGCTAAACCGGAACCAGTaaa                                          <  2:285690/149‑1 (MQ=255)
                                                                                                               gCCCCCAAAAGCGGCTGTGTTCGGCTGCGCGCATGCCCTTGGCAAGCTGGGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCGGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAgcgc                                      <  2:26345/149‑1 (MQ=255)
                                                                                                                 cgcctaaaGCGTCTGGGTTCGGCTGCGCGATTGCCCTTGGCATGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGATAACTTTTACTGCTAAACCGGAAACATTAAATCGCGa                                    <  4:17962/144‑1 (MQ=255)
                                                                                                                 ccccAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGa                                    <  3:163211/149‑1 (MQ=255)
                                                                                                                  cccAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAg                                   >  6:122975/1‑149 (MQ=255)
                                                                                                                    cAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTg                                 >  8:38404/1‑149 (MQ=255)
                                                                                                                       aaGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGt                               >  4:288960/1‑148 (MQ=255)
                                                                                                                            tCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAg                         <  6:297193/149‑1 (MQ=255)
                                                                                                                              tgtgTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTaga                                                                               >  5:28412/1‑93 (MQ=255)
                                                                                                                              tgtgTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTaga                                                                               <  6:28412/93‑1 (MQ=255)
                                                                                                                                tgtTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCg                      <  5:4451/148‑1 (MQ=255)
                                                                                                                                        tGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTg             >  1:141747/1‑149 (MQ=255)
                                                                                                                                          acagcTTGCCAAGGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTa                          >  6:350721/4‑134 (MQ=255)
                                                                                                                                          acagcTTGCCAAGGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTa                          <  5:350721/131‑1 (MQ=255)
                                                                                                                                          acagcTTGCCAAGGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTc                       >  4:137151/4‑137 (MQ=255)
                                                                                                                                          acagcTTGCCAAGGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTc                       <  3:137151/134‑1 (MQ=255)
                                                                                                                                          acagcTTGCCAAGGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGc           >  2:309104/4‑149 (MQ=255)
                                                                                                                                            cgcTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCa         <  1:156375/149‑1 (MQ=255)
                                                                                                                                            agcTTGCCAAGGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAACGCGCGAGTGGTTCGTAGTCGACGCCGCTGGc           >  2:140598/2‑147 (MQ=255)
                                                                                                                                             gcTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAg        >  1:175378/1‑149 (MQ=255)
                                                                                                                                               ttGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGAc      >  5:130176/1‑149 (MQ=255)
                                                                                                                                               ttGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGAc      >  5:203567/1‑149 (MQ=255)
                                                                                                                                                 gCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGaa                                                  >  4:364591/1‑103 (MQ=255)
                                                                                                                                                  cccTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGACCCt   <  1:70611/149‑1 (MQ=255)
                                                                                                                                                   ccTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGAcgccg                >  3:121589/1‑135 (MQ=255)
                                                                                                                                                   ccTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGAcgccg                <  4:121589/135‑1 (MQ=255)
                                                                                                                                                   ccTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGACCCTg  >  7:89833/1‑149 (MQ=255)
                                                                                                                                                   ccTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGACCCTg  <  2:90150/149‑1 (MQ=255)
                                                                                                                                                   ccTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGACCCTg  >  5:359291/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGGTTGAAATTTTTCGCGCAATCTGGATAATTCAGCCCCTCTCCGCGTGGCGGAAGTGATGCACCATCACGCAGGAAGAACCCCGTCGAACATTGGCGTGCCCGACCCGAGCCCCCAAAAGCGTCTGTGTTCGGCTGCGCGCTTGCCCTTGGCAAGCTGTGCACTATCCAGTAAGATTCGCCGTCTATTTTTCGCTGGGCGGCCTCTGAGGCTTTAGAGAATGAAAACTTTTACTGCTAAACCGGAAACAGTAAAGCGCGAGTGGTTCGTAGTCGACGCCGCTGGCCAGACCCTG  >  NC_002947/1502977‑1503272

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