Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,943 C→T 23.8% noncoding (52/120 nt) PP_5SB → 5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947181,9430CT23.8% 438.5 / 154.4 249noncoding (52/120 nt)PP_5SB5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (101/88);  new base T (57/2);  total (158/91)
Fisher's exact test for biased strand distribution p-value = 3.32e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 259 total aligned reads displayed.

AAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGA  >  NC_002947/181795‑182092
                                                                                                                                                    |                                                                                                                                                     
aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCATCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                      >  2:2150124/1‑150 (MQ=255)
  aaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGt                                                                                                                                                    >  1:3571856/1‑150 (MQ=255)
      ccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa                                                                                                                                                >  1:433167/1‑150 (MQ=255)
          tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt                                                                                                                                            >  1:1318467/1‑150 (MQ=255)
          tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt                                                                                                                                            >  2:1380382/1‑150 (MQ=255)
              catgtGCAGATAGGGGTTTTGCCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCAt                                                                                                                                        >  1:465012/5‑150 (MQ=255)
                   gCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCg                                                                                                                                   >  2:1279410/1‑150 (MQ=255)
                       gTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATgg                                                                                                                               >  2:2474725/1‑150 (MQ=255)
                                   tCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCCTTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAGGGTAGTGTGGGGtt                                                                                                                   >  1:2461976/1‑150 (MQ=255)
                                   tCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGACACGATGCATCGCCGATGGTAGTGTGGGGtt                                                                                                                   >  1:2498468/1‑150 (MQ=255)
                                   tCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGtt                                                                                                                   >  2:203536/1‑150 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                 >  2:3376639/1‑150 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                 >  1:3553969/1‑150 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                 >  1:2319935/1‑150 (MQ=255)
                                                     agagATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAgg                                                                                                 >  1:4132236/1‑150 (MQ=255)
                                                     agagATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAgg                                                                                                 >  1:2415286/1‑150 (MQ=255)
                                                     agagATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAGTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAgg                                                                                                 >  2:2512983/1‑150 (MQ=255)
                                                          ttCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCAGAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATc                                                                                            <  1:2479417/150‑1 (MQ=255)
                                                           tCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCg                                                                                           >  1:3293933/1‑150 (MQ=255)
                                                                 gCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGa                                                                                     >  1:3722853/1‑150 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaa                                                                                >  2:2145775/1‑150 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaa                                                                                >  2:2599878/1‑150 (MQ=255)
                                                                        cTAGGGACGGGCCAGCACACAGAATTTCTTGACGACAATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTACACATGTGAGAGTAGGTCATCGTCAAGATACAAAc                                                                              >  2:2984743/1‑150 (MQ=255)
                                                                          aGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAAccc                                                                            >  1:2818768/1‑150 (MQ=255)
                                                                            ggACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCaa                                                                          >  2:3289958/1‑150 (MQ=255)
                                                                            ggACGGGCCAGCACACAGAATTTCTTGACGACCAAAGAGCATTGGAACCTCCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCaa                                                                          <  1:3907343/150‑1 (MQ=255)
                                                                                  gCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGcccc                                                                    >  1:1074565/1‑150 (MQ=255)
                                                                                  gCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGcccc                                                                    >  1:2125534/1‑150 (MQ=255)
                                                                                   ccAGCACACAGTATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCt                                                                   <  1:3964891/150‑1 (MQ=255)
                                                                                      gCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTc                                                                <  1:2913012/150‑1 (MQ=255)
                                                                                      gCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTc                                                                >  2:4275369/1‑150 (MQ=255)
                                                                                         aacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTg                                                                  <  2:2740574/144‑1 (MQ=255)
                                                                                         aacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGc                                                                                                                                     >  1:3668928/2‑78 (MQ=17)
                                                                                         aacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGc                                                                                                                                     <  2:3668928/77‑1 (MQ=17)
                                                                                          acaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAATGCCCCTGTCTGCg                                                            >  2:4185149/1‑150 (MQ=255)
                                                                                          acaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCg                                                            >  1:3278675/1‑150 (MQ=255)
                                                                                             gAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                         >  1:177876/1‑150 (MQ=255)
                                                                                                tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                      >  2:2764301/1‑86 (MQ=33)
                                                                                                tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                      <  1:2764370/86‑1 (MQ=33)
                                                                                                tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTggg                                                                                                                      <  1:2764301/86‑1 (MQ=33)
                                                                                                tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAg                                                      >  2:3851567/1‑150 (MQ=255)
                                                                                                  tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGAc                                                    <  2:3553969/150‑1 (MQ=255)
                                                                                                    ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                  >  1:1878123/1‑150 (MQ=255)
                                                                                                    ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                  <  1:1380382/150‑1 (MQ=255)
                                                                                                    ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAg                                                  >  2:3475010/1‑150 (MQ=255)
                                                                                                      gacgagCATAGCGCATTGGACCCACCTGTTCCCTTCCGGACCTCAGCAGTGTCACGATGCCACGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAggg                                                <  1:895834/150‑1 (MQ=17)
                                                                                                      gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatag                                                >  1:593120‑M2/1‑117 (MQ=255)
                                                                                                       acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagg                                               >  2:1065118‑M2/1‑116 (MQ=255)
                                                                                                       acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgtgcagatagg                                               >  2:2374559‑M2/1‑116 (MQ=255)
                                                                                                       acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgtgcagatagg                                               >  2:1933603‑M2/1‑116 (MQ=255)
                                                                                                       acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                               >  2:2765734‑M2/1‑116 (MQ=255)
                                                                                                       acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                               >  2:4390957‑M2/1‑116 (MQ=255)
                                                                                                        cgacCATCGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTGTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                              >  1:2743753‑M2/1‑115 (MQ=255)
                                                                                                        cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                              >  1:1195477‑M2/1‑115 (MQ=255)
                                                                                                        cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                              >  1:746763‑M2/1‑115 (MQ=255)
                                                                                                        cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggg                                              >  1:3817213‑M2/1‑115 (MQ=255)
                                                                                                        cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgtgcagataggg                                              >  2:1696228‑M2/1‑115 (MQ=255)
                                                                                                        cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgtgcagataggg                                              >  1:1405660‑M2/1‑115 (MQ=255)
                                                                                                        cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                              >  2:3596559‑M2/1‑115 (MQ=255)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagggg                                             >  1:4191449‑M2/1‑114 (MQ=255)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagggg                                             >  1:3790065‑M2/1‑114 (MQ=255)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagggg                                             >  1:1513348‑M2/1‑114 (MQ=255)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagggg                                             >  1:2599598‑M2/1‑114 (MQ=255)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATc                                                                                                                                       >  2:4015323/1‑60 (MQ=35)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATc                                                                                                                                       <  1:4015690/60‑1 (MQ=35)
                                                                                                          acCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggggt                                            >  2:1307233‑M2/1‑113 (MQ=255)
                                                                                                          acCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggt                                            >  1:1011411‑M2/1‑113 (MQ=255)
                                                                                                           ccAGAGCGCCTTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCAGCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                           <  2:919869/150‑1 (MQ=255)
                                                                                                             aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggtttt                                         >  2:4488698‑M2/1‑110 (MQ=255)
                                                                                                              tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggttttg                                        >  2:1961391‑M2/1‑109 (MQ=255)
                                                                                                              tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggttttg                                        >  1:698181‑M2/1‑109 (MQ=255)
                                                                                                               agagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaacccctacctgcgtgcgcagataggggttttgc                                       >  1:2576590‑M2/1‑108 (MQ=255)
                                                                                                               agagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaacccctacctgcgtgcgcagataggggttttgc                                       >  2:1990000‑M2/1‑108 (MQ=255)
                                                                                                               agagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAtgt                                                                                                          <  1:1974153/83‑1 (MQ=33)
                                                                                                                gagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt                                                                                                                                            <  1:1283003/48‑1 (MQ=38)
                                                                                                                   cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggttttgcgaaa                                   >  2:2207838‑M2/1‑104 (MQ=255)
                                                                                                                   cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaaccccgacctgcatgcgcagataggggttttgcgaaa                                   >  2:760760‑M2/1‑104 (MQ=255)
                                                                                                                   cATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa                                                                                                                                  >  2:3270557/1‑55 (MQ=35)
                                                                                                                       ggAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCg                                                                                           <  2:591641/90‑1 (MQ=33)
                                                                                                                        gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                              <  2:1110392/150‑1 (MQ=255)
                                                                                                                        gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                              <  2:4464342/150‑1 (MQ=255)
                                                                                                                        gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                              >  1:1652/1‑150 (MQ=255)
                                                                                                                         aaCCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTgg                                                                                                                       <  2:3239458/60‑1 (MQ=35)
                                                                                                                           ccaccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTa                                                                                                                             <  1:4058339/52‑1 (MQ=35)
                                                                                                                           ccaccTGACCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTc                           <  2:4094156/150‑1 (MQ=255)
                                                                                                                            caccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGGAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCa                          >  2:374334/1‑150 (MQ=255)
                                                                                                                             accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGAt                                                          >  1:420788/1‑117 (MQ=255)
                                                                                                                              ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAat                        <  2:465012/150‑1 (MQ=255)
                                                                                                                              ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCGTTGTCTTTTCCGGGTTTCAat                        >  1:2037264/1‑150 (MQ=255)
                                                                                                                                  atcccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCATTGTCTTTTCCGGGTTTCAATATCg                    <  2:2489700/150‑1 (MQ=255)
                                                                                                                                    cccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                  <  2:1603089/150‑1 (MQ=255)
                                                                                                                                     ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt                 >  2:599712/1‑150 (MQ=255)
                                                                                                                                               cTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGTTTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATa       <  2:223758/150‑1 (MQ=255)
                                                                                                                                               cTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATa       <  2:1853806/150‑1 (MQ=255)
                                                                                                                                               cTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATa       <  1:3179250/150‑1 (MQ=255)
                                                                                                                                                tCAGCAGTGCCACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTTGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCATTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAg      <  2:3005022/150‑1 (MQ=255)
                                                                                                                                                tCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAg      >  1:599447/1‑150 (MQ=255)
                                                                                                                                                  agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAATCCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt    >  1:2002934/1‑150 (MQ=255)
                                                                                                                                                  agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt    >  1:1737562/1‑150 (MQ=255)
                                                                                                                                                  agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt    >  2:710582/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGA  >  NC_002947/181795‑182092

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