Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,818 G→A 6.8% intergenic (+32/‑74) PP_r01 → / → PP_5SB 5S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947181,8180GA 6.8% 282.2 / 8.9 103intergenic (+32/‑74)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (49/46);  new base A (5/2);  total (54/49)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 103 total aligned reads displayed.

TGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGA  >  NC_002947/181671‑181964
                                                                                                                                                   |                                                                                                                                                  
tGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTa                                                                                                                                                  <  2:541084/150‑1 (MQ=32)
   cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                               <  1:320901/150‑1 (MQ=32)
       cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGttt                                                                                                                                           >  1:1241503/1‑150 (MQ=32)
        aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtttt                                                                                                                                          <  1:1364594/150‑1 (MQ=32)
          agagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGt                                                                                                                                        <  2:2214057/150‑1 (MQ=32)
                ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                                  <  1:668204/150‑1 (MQ=32)
                ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                                  <  2:1116809/150‑1 (MQ=32)
                ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                                  >  2:259554/1‑150 (MQ=32)
                ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                                  <  2:1241503/150‑1 (MQ=32)
                   gAACCACCTGATCCCAGCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTa                                                                                                                               >  2:101951/1‑150 (MQ=17)
                          cTGATCCCATCCCGAACTCAGTTGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTa                                                                                                                        >  2:2221022/1‑150 (MQ=38)
                          cTGATCCCATCCCGAACTCAGTTGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTa                                                                                                                        >  2:2220519/1‑150 (MQ=38)
                          cTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGTTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTa                                                                                                                        >  1:1453280/1‑150 (MQ=37)
                            gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACa                                                                                                                      <  1:101951/150‑1 (MQ=39)
                                ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAgaga                                                                                                                  <  2:1151548/150‑1 (MQ=255)
                                 catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGAt                                                                                                                 >  1:1721254/1‑150 (MQ=255)
                                     ccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGGAGGTCATCGTCAAGATTCATTCCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGc                                                                                                             >  2:458912/1‑150 (MQ=255)
                                      cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACag                                                                                                                     >  2:870343/1‑141 (MQ=40)
                                      cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACag                                                                                                                     <  1:870343/141‑1 (MQ=40)
                                               tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTcc                                                                                                   >  2:2249241/1‑150 (MQ=255)
                                                     aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTGCAGATTCATTTCACAAAACCCCTATCTGAACATGCAGGTAGGGGTTTTGTCTTTAGGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGa                                                                                             >  1:1035712/1‑150 (MQ=18)
                                                     aaCGATGCATCGCCGATGGTAGTGTGGGGATTCCCCATGTGAGAGTAGGTCATCGTCAGGGTTCATTTCGCAAAACCTCTATCTGCACATGCAGGTAGTGGATTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGa                                                                                             >  2:2048788/1‑150 (MQ=18)
                                                      aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGAc                                                                                            <  1:961067/150‑1 (MQ=255)
                                                            cATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCa                                                                                      <  2:1005809/150‑1 (MQ=255)
                                                            cATCGACGCTGGTAGTGTGGCGCATACCCTTGTGAGAGTACGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCa                                                                                      <  2:1721254/150‑1 (MQ=18)
                                                              tCGCCGATGGTAGTGTGGGGTTTCCCTATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGaa                                                                                                                           <  2:1595144/111‑1 (MQ=25)
                                                              tCGCCGATGGTAGTGTGGGGTTTCCCTATGTGAGAGTAGGTCATCGTCAAGATTCAGTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGaa                                                                                                                           >  1:1595144/1‑111 (MQ=255)
                                                                  cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacac                                                                                <  2:1901704/150‑1 (MQ=255)
                                                                  cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacac                                                                                <  1:2221022/150‑1 (MQ=255)
                                                                  cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacac                                                                                <  1:2220519/150‑1 (MQ=255)
                                                                     tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAgg                                                                                                                             >  1:1514034/1‑102 (MQ=35)
                                                                     tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAgg                                                                                                                             <  2:1514034/102‑1 (MQ=35)
                                                                       gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCCTCGTATCACTCCAAGGCCAAAAACCCCTATCTGCCCAGTCATGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTGCTAGGGACGGGCCAGCACACAGAAt                                                                           >  1:1022613/1‑150 (MQ=21)
                                                                        tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAtt                                                                          >  2:1437405/1‑150 (MQ=255)
                                                                         aGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCTCTATCTGCACCTGGATTTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAttt                                                                         <  2:1771400/150‑1 (MQ=21)
                                                                          gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg                                                                                     <  2:1773976/137‑1 (MQ=255)
                                                                          gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg                                                                                     >  1:1773976/1‑137 (MQ=255)
                                                                          gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTc                                                                        <  1:1673759/150‑1 (MQ=255)
                                                                             tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTGGGGACGGGCCAGCACACAGAATTTCTTg                                                                     >  1:707126/1‑150 (MQ=255)
                                                                             tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTg                                                                     >  2:1453464/1‑150 (MQ=255)
                                                                              ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                                    <  1:2112315/150‑1 (MQ=255)
                                                                                ggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATTTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGTGACGGGCGAGCACACAGAATTTATTgaca                                                                  >  2:1820665/1‑149 (MQ=255)
                                                                                 gTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacga                                                                 >  2:2352533/1‑150 (MQ=255)
                                                                                 gTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacga                                                                 >  2:1973931/1‑150 (MQ=255)
                                                                                 gTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGCCCTAGGGACGGGCCAGCACACAGAATTTCTTgacga                                                                 >  1:14168/1‑150 (MQ=255)
                                                                                   ctCCCCATGTGACAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGAGTCCGATTAGCACACAGAATTTCTTGACGAcc                                                               <  1:675236/149‑1 (MQ=21)
                                                                                    tCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCNTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGAACGTCCTAGGGAGGGGCAAGCACACAGAATTTCTTGACGAGCa                                                              >  1:1189510/1‑150 (MQ=255)
                                                                                      cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                                 >  2:663994/1‑65 (MQ=37)
                                                                                      cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                                 <  1:663994/65‑1 (MQ=37)
                                                                                      cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACCGAGAAGGGAAAGCGTGTGCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCaag                                                            >  2:675236/1‑148 (MQ=18)
                                                                                      cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                                            >  1:1428079/1‑150 (MQ=255)
                                                                                      cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                                            <  2:1028121/150‑1 (MQ=255)
                                                                                      cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                                            >  2:641197/1‑150 (MQ=255)
                                                                                      cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACACCTATGTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTATGAACTACAGAGATTCGCAGGCACGTACTAGGGAGGGGCTAGTACACAGAATTTCTTGACGACCATa                                                            >  2:1187585/1‑150 (MQ=21)
                                                                                       ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCa                                                                                      >  2:620316/1‑123 (MQ=255)
                                                                                       ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCa                                                                                      <  1:620316/123‑1 (MQ=255)
                                                                                         aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAc                                                                                                                          <  2:1125870/85‑1 (MQ=40)
                                                                                         aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAc                                                                                                                          >  1:1125870/1‑85 (MQ=40)
                                                                                             gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                     <  1:2352533/150‑1 (MQ=255)
                                                                                                agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGa                                                  <  2:1428079/150‑1 (MQ=255)
                                                                                                     gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccac                                             <  1:1973931/150‑1 (MQ=255)
                                                                                                      tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa                                                                                                            >  2:1268610/1‑86 (MQ=255)
                                                                                                      tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa                                                                                                            <  1:1268610/86‑1 (MQ=255)
                                                                                                      tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAATTCTTGACGAACATAGAGCATTGGAaccacc                                            >  2:483109/1‑150 (MQ=255)
                                                                                                       cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                     >  2:1448352/1‑140 (MQ=255)
                                                                                                       cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                     >  2:1493999/1‑140 (MQ=255)
                                                                                                       cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                     <  1:1493999/140‑1 (MQ=255)
                                                                                                       cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                     <  1:1448352/140‑1 (MQ=255)
                                                                                                       cATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCt                                           >  2:228187/1‑150 (MQ=255)
                                                                                                           gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                                       <  2:1022613/150‑1 (MQ=255)
                                                                                                            tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                                      >  2:2280000/1‑118 (MQ=255)
                                                                                                            tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                                      <  1:2280000/118‑1 (MQ=255)
                                                                                                            tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                                      >  2:2191780/1‑150 (MQ=255)
                                                                                                               aGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                                      <  2:691168/115‑1 (MQ=255)
                                                                                                               aGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                                      >  1:691168/1‑115 (MQ=255)
                                                                                                                 aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt                                                                       >  2:1161842/1‑112 (MQ=255)
                                                                                                                 aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt                                                                       <  1:1161842/112‑1 (MQ=255)
                                                                                                                 aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccac                                             >  2:1240191/1‑138 (MQ=255)
                                                                                                                 aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccac                                             <  1:1240191/138‑1 (MQ=255)
                                                                                                                 aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatcc                                 >  2:1485493/1‑150 (MQ=255)
                                                                                                                     aTTCCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGCAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                             <  1:458912/150‑1 (MQ=255)
                                                                                                                      aaagacAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAacc                                               <  2:776149/126‑1 (MQ=37)
                                                                                                                      aaagacAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAacc                                               >  1:776149/6‑131 (MQ=38)
                                                                                                                        tCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGAcc                                                               >  2:2288039/1‑113 (MQ=255)
                                                                                                                        tCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGAcc                                                               <  1:2288039/113‑1 (MQ=255)
                                                                                                                          gCAAAACCCCTATCTGCACATGGAGGTAGGGGTTTGGTCGTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag                        <  2:2300451/150‑1 (MQ=255)
                                                                                                                          gCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag                        >  2:214020/1‑150 (MQ=255)
                                                                                                                            aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                      >  1:1663503/1‑150 (MQ=255)
                                                                                                                            aaaaCCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAAc                            >  2:2128027/1‑144 (MQ=37)
                                                                                                                            aaaaCCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAAc                            <  1:2128027/144‑1 (MQ=37)
                                                                                                                                 cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaa                 >  2:2213761/1‑150 (MQ=255)
                                                                                                                                 cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaa                 >  1:1200341/1‑150 (MQ=37)
                                                                                                                                   cTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTc                                                                                                    >  2:2141285/1‑65 (MQ=255)
                                                                                                                                   cTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTc                                                                                                    <  1:2141285/65‑1 (MQ=255)
                                                                                                                                   cTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaaa               >  1:754518/1‑149 (MQ=255)
                                                                                                                                      tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt            <  1:214020/150‑1 (MQ=255)
                                                                                                                                      tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt            <  1:1453464/150‑1 (MQ=255)
                                                                                                                                        tGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGc          >  1:587720/1‑150 (MQ=255)
                                                                                                                                           acaTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATc       >  1:1482174/1‑150 (MQ=255)
                                                                                                                                                cAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa  <  1:1485493/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGA  >  NC_002947/181671‑181964

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