Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,805 T→C 11.5% intergenic (+19/‑87) 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,8050TC11.5% 367.2 / 30.7 145intergenic (+19/‑87)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (45/79);  new base C (8/8);  total (58/87)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 145 total aligned reads displayed.

GAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAC  >  NC_002947/181663‑181951
                                                                                                                                              |                                                                                                                                                  
  aTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGACATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATg                                                                                                                                           >  1:2678115/1‑150 (MQ=17)
   tttCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGc                                                                                                                                          <  2:2107000/150‑1 (MQ=32)
      cTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAgg                                                                                                                                       <  2:1600566/150‑1 (MQ=32)
       ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt                                                                                                                                      >  1:1301010/1‑150 (MQ=32)
        tGACGACCAGCGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTa                                                                                                                                     <  2:3031438/150‑1 (MQ=17)
         gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCGAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                    <  2:1973853/150‑1 (MQ=17)
          acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAgg                                                                                                                                   <  1:3140430/150‑1 (MQ=32)
          acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAgg                                                                                                                                   <  2:2653809/150‑1 (MQ=32)
          acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAgg                                                                                                                                   <  2:2653897/150‑1 (MQ=32)
           cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCatcttctg                                                                                                                                  >  2:3401346/1‑143 (MQ=17)
           cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCatcttctg                                                                                                                                  >  2:3401196/1‑143 (MQ=17)
           cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                  <  2:4982716/150‑1 (MQ=32)
            gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAgggg                                                                                                                                 <  1:4798406/150‑1 (MQ=255)
              ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtt                                                                                                                               <  1:4464662/150‑1 (MQ=32)
               cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTTGTGAAACGATGCATCGCCGATGGGAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAGTTCGCTATACCCGTGTCTGCACGTGCGGGGGGGGGtgt                                                                                                                              >  2:2058755/1‑148 (MQ=11)
               cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                  <  1:1995392/146‑1 (MQ=32)
               cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                  >  2:1995392/1‑146 (MQ=32)
                   gagCATTGGAACCTCCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTc                                                                                                                          <  1:791397/150‑1 (MQ=17)
                   gagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTc                                                                                                                          <  1:3439046/150‑1 (MQ=32)
                   gagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGCCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTc                                                                                                                          <  1:916304/150‑1 (MQ=17)
                      cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTgcgaaa                                                                                                                       <  2:2711087/150‑6 (MQ=11)
                      cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTgcgaaa                                                                                                                       <  1:3016864/150‑6 (MQ=11)
                       aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTgcgaaat                                                                                                                      <  1:5207743/150‑7 (MQ=11)
                           gAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCCCAAAACCCATAAATCCAAATGCACGTAGGGGTTGTGTATTTaataa                                                                                                                  >  2:3406173/1‑147 (MQ=11)
                            aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAg                                                                                                                 <  1:4884493/150‑1 (MQ=32)
                            aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAg                                                                                                                 <  1:497641/150‑1 (MQ=32)
                             accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAgg                                                                                                                <  1:4794788/150‑1 (MQ=32)
                             accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAgg                                                                                                                <  2:4548451/150‑1 (MQ=32)
                               tccgagctcccatcccGAACTCAGTTGTGAAACGTTTCATCGCCTTTTGTTGTGTGGGGTTTCCCCATGTGAGTGTAGGTCATCGTCAAGTTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGTGTTTTGTCTTTCAGTAGGaa                                                                                                              <  2:2088961/143‑1 (MQ=11)
                               caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGaa                                                                                                              <  2:1591406/150‑1 (MQ=34)
                               caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGaa                                                                                                              <  1:4948196/150‑1 (MQ=34)
                                     atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTCTCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACag                                                                                                        <  2:2901204/150‑1 (MQ=255)
                                     atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACag                                                                                                        <  2:2901088/150‑1 (MQ=39)
                                      tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                       >  1:2064579/1‑150 (MQ=40)
                                      tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                       <  2:970728/150‑1 (MQ=40)
                                      tcccatcccGAACTCAGTAGTGAAACGATACATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                       <  2:970766/150‑1 (MQ=255)
                                           tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGATAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTc                                                                                                  <  1:5252524/150‑1 (MQ=255)
                                              cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa                                                                                               >  1:2891487/1‑150 (MQ=255)
                                               gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAg                                                                                              <  1:2813746/150‑1 (MQ=255)
                                                      gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTc                                                                                       <  1:2933692/150‑1 (MQ=255)
                                                          tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtttt                                                                                                                             >  2:3268682/1‑108 (MQ=32)
                                                          tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtttt                                                                                                                             <  1:3268682/108‑1 (MQ=32)
                                                            aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAATTCGCAAAACCCCTAACTGCAAATGCAGGGAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAggg                                                                                 >  2:3880081/1‑150 (MQ=255)
                                                             aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGa                                                                                <  1:3627479/150‑1 (MQ=255)
                                                                gATGCATCTCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAATATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGTGACgg                                                                             <  2:27527/150‑1 (MQ=255)
                                                                gATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                             <  1:2855843/150‑1 (MQ=255)
                                                                  tGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGc                                                                           <  1:3291734/150‑1 (MQ=255)
                                                                      tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTa                                                                                                                  >  1:3444335/1‑107 (MQ=32)
                                                                      tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTa                                                                                                                  <  2:3444335/107‑1 (MQ=32)
                                                                      tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGc                                                                       >  1:4925311/1‑150 (MQ=255)
                                                                        gCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCTGGCCCGTCCCTAGGACGTGCCTGcga                                                                     <  2:1878530/150‑3 (MQ=17)
                                                                        gCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTCTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCCCGTCCCTAGGACGTGCCTGcga                                                                     <  2:216528/150‑3 (MQ=17)
                                                                              ggTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGaa                                                               <  1:1589850/150‑1 (MQ=255)
                                                                                tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAtt                                                             >  1:3824172/1‑150 (MQ=255)
                                                                                 aGTGTGGGGGGTCCCCATGTGAGAGGAGGTCATCGTCAAGATTCATTTCGCAACACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTCAGGAGGCACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAttt                                                            <  2:4231575/150‑1 (MQ=21)
                                                                                   tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTAAATGCACAAGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGATATACGCAGGAACGTCCTAGGGACGGGCCAGCACACATAATTTCt                                                          >  1:4968474/1‑150 (MQ=21)
                                                                                      ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                       <  2:2662292/150‑1 (MQ=255)
                                                                                           ttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGAcc                                                  <  1:5413301/150‑1 (MQ=255)
                                                                                              cccATGTGAGAGTAGGTCATCGTCAAGATTCAATACCAAATACCCCTATCTGCACTAGCAGGTAGGGGTTTTGTCTTTAAGTATGAAATACAGATATTAGCAGGCACGTCCTAGGGACGGGCCACAACACAGAATATCTTGACGACCATa                                               >  1:1987471/1‑150 (MQ=21)
                                                                                              cccATGTGAGAGTAGGTCATCGACAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGAATTGTCTTTAAGTACGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                               <  2:836161/150‑1 (MQ=255)
                                                                                               ccATGTGAGAGTAGGTCATCGTCAAGATTCATTACGCAAAACCCCTAACTGCACCGGCAGGTAGGGGGGGTGTCTTGAAGTCGGACCTACAGAGTTTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTACTTGACGACCATag                                              >  2:3966316/1‑150 (MQ=25)
                                                                                                cATGTGAGAGTAGGTCATCGTCAAGATTCATGTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGAGTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGCCTTTCTTGACGACCATaga                                             <  2:3601961/150‑1 (MQ=255)
                                                                                                 aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCCGGTCGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATagag                                            >  1:4431608/1‑150 (MQ=255)
                                                                                                     gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                        <  1:504167/150‑1 (MQ=255)
                                                                                                     gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                        <  1:503835/150‑1 (MQ=255)
                                                                                                     gagaGTAGGTCATCGTCAAGATTCATTTCGAAAAACCCATATATGCACATGCAGTTATGGGTTTTTTCTTTAAGTAGGAACTACAGATATTCGCAGGCACGTCCTATGGACGGGACAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                        >  2:194007/1‑150 (MQ=25)
                                                                                                      agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTAGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTg                                       >  2:605052/1‑150 (MQ=255)
                                                                                                             gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAGTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGTCGACCATAGAGCATTGGAaccag                                >  1:4513580/1‑149 (MQ=255)
                                                                                                             gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccac                                >  1:2065641/1‑150 (MQ=255)
                                                                                                             gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCAAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccac                                >  1:4513637/1‑150 (MQ=255)
                                                                                                                  cGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                           <  2:5319882/150‑1 (MQ=255)
                                                                                                                  cGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                           <  1:1345400/150‑1 (MQ=255)
                                                                                                                   gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGaa                                                               >  1:3098497/1‑113 (MQ=255)
                                                                                                                   gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGaa                                                               <  2:3098497/113‑1 (MQ=255)
                                                                                                                   gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                          >  2:15714/1‑150 (MQ=255)
                                                                                                                   gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                          <  1:1687425/150‑1 (MQ=255)
                                                                                                                    tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                         >  1:1900227/1‑150 (MQ=255)
                                                                                                                    tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                         >  2:4838859/1‑150 (MQ=255)
                                                                                                                      aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatccca                       >  1:2767697/1‑150 (MQ=255)
                                                                                                                      aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatccca                       >  1:41017/1‑150 (MQ=255)
                                                                                                                         aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatcc                    >  1:113654/1‑150 (MQ=255)
                                                                                                                         aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatcc                    <  2:113654/150‑1 (MQ=255)
                                                                                                                          ttCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatccc                   <  2:2465631/150‑1 (MQ=255)
                                                                                                                            cATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGa                 >  2:3809942/1‑150 (MQ=255)
                                                                                                                              tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGaa                                                               <  1:2019740/102‑1 (MQ=255)
                                                                                                                                  gCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag           >  1:3693152/1‑150 (MQ=255)
                                                                                                                                    aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca         >  2:2054926/1‑150 (MQ=255)
                                                                                                                                    aaaaCCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca         >  1:1064034/1‑150 (MQ=255)
                                                                                                                                        ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatcc                    <  2:4172347/135‑1 (MQ=255)
                                                                                                                                        ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatcc                    >  1:4172347/1‑135 (MQ=255)
                                                                                                                                         cccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACag                                                                                                        <  2:3996485/50‑1 (MQ=255)
                                                                                                                                         cccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACag                                                                                                        >  1:3996485/1‑50 (MQ=255)
                                                                                                                                         cccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaa    >  2:3819930/1‑150 (MQ=255)
                                                                                                                                         cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGaa                                    >  1:111884/1‑118 (MQ=39)
                                                                                                                                         cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGaa                                    <  2:111884/118‑1 (MQ=39)
                                                                                                                                         cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaa    >  2:1690604/1‑150 (MQ=255)
                                                                                                                                          ccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                >  2:1732589/1‑137 (MQ=255)
                                                                                                                                          ccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                <  1:1732589/137‑1 (MQ=255)
                                                                                                                                           cTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc  >  2:3197409/1‑150 (MQ=255)
                                                                                                                                           cTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc  >  1:4236996/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAC  >  NC_002947/181663‑181951

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