Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,943 C→T 22.4% 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,9430CT22.4% 174.1 / 48.0 98noncoding (52/120 nt)PP_5SB5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base T (13/9);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 107 total aligned reads displayed.

AAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGA  >  NC_002947/181795‑182092
                                                                                                                                                    |                                                                                                                                                      
aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                       <  1:1038483/150‑1 (MQ=255)
aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                       <  2:1704027/150‑1 (MQ=255)
aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                                                                                                                                                       <  2:1665384/150‑1 (MQ=255)
 aaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagcag                                                                                                                                                      >  1:1279059/1‑150 (MQ=255)
     cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaa                                                                                                                                                  >  2:31893/1‑150 (MQ=37)
      ccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa                                                                                                                                                 <  1:333974/150‑1 (MQ=255)
             gCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCa                                                                                                                                          >  2:436256/1‑150 (MQ=255)
                caTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCg                                                                                                                                       <  2:1279059/150‑1 (MQ=255)
                      ggTAGGGGTTTTGTGTTTAAGTAGGAACTAAAGAGATGTGCCGGCACGTACTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGAATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACTATGCATCGCCGATg                                                                                                                                 >  2:295755/1‑150 (MQ=21)
                       gTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATgg                                                                                                                                >  2:1756976/1‑150 (MQ=255)
                              ttttGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa                                                                                                                                   <  1:354374/140‑1 (MQ=255)
                              ttttGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa                                                                                                                                   >  2:354374/1‑140 (MQ=255)
                                   tCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGtt                                                                                                                    >  2:357640/1‑150 (MQ=255)
                                    cTTTAAGTAGGAACTACAGAGATTCGCATGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGttt                                                                                                                   >  1:461203/1‑150 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTc                                                                                                                  >  2:589909/1‑150 (MQ=255)
                                     tttAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCAAACAGAATTTCTTGAGGACCATAGAGCATTGGGACCACCTGATGCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGCTTc                                                                                                                  >  2:53047/1‑150 (MQ=21)
                                         aGTAGGAACTAAAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCa                                                                                                              >  2:767106/1‑150 (MQ=255)
                                          gTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAt                                                                                                             >  1:585063/1‑150 (MQ=255)
                                          gTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATGTCTGGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAt                                                                                                             >  2:1626548/1‑150 (MQ=255)
                                                   aCAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc                                                                                                                                                <  1:1040068/106‑1 (MQ=255)
                                                   aCAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc                                                                                                                                                >  2:1040068/1‑106 (MQ=255)
                                                   aCAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTa                                                                                                    >  1:359518/1‑150 (MQ=255)
                                                    cAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAg                                                                                                   <  2:1554399/150‑1 (MQ=255)
                                                    cAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCAACCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAg                                                                                                   <  2:1396062/150‑1 (MQ=255)
                                                    cAGAGATTCGCAGGCACGTCCTAGGGACCGGCCACCACCCAGCATTTCTTGCCGACCACAGAGCATTGGAACCACCTAATCTCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAg                                                                                                   <  1:1401038/150‑1 (MQ=17)
                                                     agagATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAgg                                                                                                  <  1:357640/150‑1 (MQ=255)
                                                      gagaTTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGCATTTCTTGACGACCATAGAGCATTGGAACCACCTCATCCCATCCCCAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTATGGGGTTTCCCCATGTGAGAGTAGGt                                                                                                 <  1:53047/150‑1 (MQ=255)
                                                               aGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCaa                                                                                        <  1:589909/150‑1 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCCTTGCAAGCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaa                                                                                 >  1:49872/1‑150 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa                                                                                                                                   <  1:564829/100‑1 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa                                                                                                                                   >  2:564829/1‑100 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGagag                                                                                                      >  2:591099/1‑129 (MQ=255)
                                                                      tCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaa                                                                                 >  1:77141/1‑150 (MQ=255)
                                                                            ggACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCaa                                                                           >  1:606000/1‑150 (MQ=255)
                                                                                gggCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAATGcc                                                                       <  2:585063/150‑1 (MQ=255)
                                                                                 ggCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGccc                                                                      <  2:606000/150‑1 (MQ=255)
                                                                                 ggCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGccc                                                                      <  2:461203/150‑1 (MQ=255)
                                                                                    cAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTa                                                                                                    >  1:857272/1‑117 (MQ=40)
                                                                                    cAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTa                                                                                                    <  2:857272/117‑1 (MQ=40)
                                                                                     aGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCGCCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGCGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGt                                                                  >  1:1844670/1‑150 (MQ=255)
                                                                                      gCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATg                                                                                                                                            <  1:884193/75‑1 (MQ=255)
                                                                                      gCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATg                                                                                                                                            >  2:884193/1‑75 (MQ=255)
                                                                                         aacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGt                                                                                                                                                     <  1:802422/62‑1 (MQ=18)
                                                                                         aacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGt                                                                                                                                                     >  2:802422/2‑63 (MQ=18)
                                                                                          acaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAt                                                                                                             <  1:128817/102‑1 (MQ=33)
                                                                                          acaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTTTCTGCg                                                             >  2:1723191/1‑150 (MQ=255)
                                                                                           caGAATTTCTTGACGACCATAGCGAAAGGTAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGa                                                            >  1:1232943/1‑150 (MQ=17)
                                                                                            aGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGcatc                                                           >  2:1729906/1‑147 (MQ=255)
                                                                                              aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaaa                                                                                <  1:1513505/127‑1 (MQ=35)
                                                                                              aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCaaa                                                                                >  2:1513505/1‑127 (MQ=35)
                                                                                                   cTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACa                                                   <  1:8477/150‑1 (MQ=255)
                                                                                                    ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAg                                                  <  1:1729906/150‑1 (MQ=255)
                                                                                                     tGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgg                                                 <  2:359518/150‑1 (MQ=255)
                                                                                                     tGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                                 <  1:1756976/150‑1 (MQ=255)
                                                                                                       acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagg                                                >  2:1307742‑M2/1‑116 (MQ=255)
                                                                                                        cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                               >  1:1648093‑M2/1‑115 (MQ=255)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagggg                                              >  1:102419‑M2/1‑114 (MQ=255)
                                                                                                         gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagggg                                              >  1:1057524‑M2/1‑114 (MQ=255)
                                                                                                           ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                           >  1:540120/1‑150 (MQ=255)
                                                                                                           ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                           <  2:1844670/150‑1 (MQ=255)
                                                                                                            cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                               <  1:1941153‑M2/146‑36 (MQ=255)
                                                                                                            cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                               >  2:1941153‑M2/1‑111 (MQ=255)
                                                                                                            cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatag                                                 >  2:853885‑M2/1‑111 (MQ=255)
                                                                                                            cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatag                                                 <  1:853885‑M2/144‑34 (MQ=255)
                                                                                                            cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcacaaaacccctacctgcatgtgcagataggggttt                                           <  1:1761408‑M2/150‑40 (MQ=255)
                                                                                                                gagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGAtt                                                                                    <  2:1355921/105‑1 (MQ=33)
                                                                                                                gagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGAtt                                                                                    >  1:1355921/1‑105 (MQ=33)
                                                                                                                         aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggttt                                           <  1:1942398‑M2/137‑40 (MQ=255)
                                                                                                                         aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggggttt                                           >  2:1942398‑M2/1‑98 (MQ=255)
                                                                                                                          accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                               >  1:648731‑M2/1‑97 (MQ=255)
                                                                                                                          accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                               <  2:648731‑M2/132‑36 (MQ=255)
                                                                                                                           ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                                >  2:1615334‑M2/1‑96 (MQ=255)
                                                                                                                           ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagatagg                                                <  1:1615334‑M2/130‑35 (MQ=255)
                                                                                                                              ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggggtttt                                          <  2:1331481‑M2/133‑41 (MQ=255)
                                                                                                                              ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagataggggtttt                                          >  1:1331481‑M2/1‑93 (MQ=255)
                                                                                                                              ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACAggg                                                >  2:1051075/1‑126 (MQ=14)
                                                                                                                              ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACACAAAGCCCCTGTCTGCATCGCAGACAggg                                                <  1:1051075/126‑1 (MQ=17)
                                                                                                                              ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAATAGGTCATCGTCAAGATTCAAACCAGAAGGGACTGGCAGCGGTGTGGACAGGGGCTTTGTCTTTTCCGGGTTTCAat                        >  1:1098089/1‑150 (MQ=17)
                                                                                                                               cTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAtgt                                                                                                           >  1:1808449/1‑67 (MQ=37)
                                                                                                                                 gATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTga                                                                                                         <  1:1156669/67‑1 (MQ=35)
                                                                                                                                  atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                               >  2:269906‑M2/1‑89 (MQ=255)
                                                                                                                                  atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcatgcgcagataggg                                               <  1:269906‑M2/124‑36 (MQ=255)
                                                                                                                                     ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGGTGGTAGTGTGGGGTTTCCCCATGTGAGAGGAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAAACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGtt                 >  1:1564721/1‑150 (MQ=255)
                                                                                                                                     ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTcc                                                                                                                 >  1:33733/1‑55 (MQ=35)
                                                                                                                                     ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTcc                                                                                                                 <  2:33733/55‑1 (MQ=38)
                                                                                                                                     ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTATTTTCAGGGTTTCTATATCTGtt                 >  2:415328/1‑150 (MQ=255)
                                                                                                                                      catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagg                                                <  1:1001781‑M2/119‑35 (MQ=255)
                                                                                                                                      catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAtttcgcaaaacccctacctgcgtgcgcagatagg                                                >  2:1001781‑M2/1‑85 (MQ=255)
                                                                                                                                          ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAtgt                                                                                                           <  2:629544/56‑1 (MQ=38)
                                                                                                                                          ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCAtgt                                                                                                           >  1:629544/1‑56 (MQ=38)
                                                                                                                                           cGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGa                                                                                      <  1:1879288/76‑1 (MQ=35)
                                                                                                                                             aaCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa         <  1:1626548/150‑1 (MQ=255)
                                                                                                                                              aCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACt           <  2:376932/147‑1 (MQ=255)
                                                                                                                                              aCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACt           >  1:376932/1‑147 (MQ=255)
                                                                                                                                              aCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTCTCCCCATGTGAGAGTAGGTCATCGTCAAGCTTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGAt        <  1:436256/150‑1 (MQ=255)
                                                                                                                                                tCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAg      <  1:1084948/150‑1 (MQ=255)
                                                                                                                                                  agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCg                                 >  1:894396/1‑121 (MQ=255)
                                                                                                                                                  agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt    <  2:1474145/150‑1 (MQ=255)
                                                                                                                                                    cagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGa  >  2:1033944/1‑150 (MQ=255)
                                                                                                                                                    cagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCCGGTTTCAATATCGGTTAGAACTGATAGCTGa  <  2:540120/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGA  >  NC_002947/181795‑182092

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: