Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,805 T→C 9.0% 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,8050TC 9.0% 893.4 / 55.7 312intergenic (+19/‑87)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (95/188);  new base C (9/19);  total (104/208)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

Only 100 of 315 total aligned reads displayed.

AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181662‑181954
                                                                                                                                               |                                                                                                                                                     
aGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGcat                                                                                                                                                 >  1:3672969/1‑149 (MQ=11)
   aTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGcatgcg                                                                                                                                              <  1:2506098/150‑5 (MQ=11)
      tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAg                                                                                                                                           <  2:7714413/150‑1 (MQ=32)
      tCTTGACGACCAGAGAGCATTGGAACCACCTGAGCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGcatgcgcag                                                                                                                                           <  2:3609193/150‑8 (MQ=11)
       cTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAgg                                                                                                                                          <  2:4103910/150‑1 (MQ=32)
        ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAgat                                                                                                                                         >  1:2374579/1‑148 (MQ=9)
         tGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTa                                                                                                                                        <  1:2007610/150‑1 (MQ=32)
         tGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTa                                                                                                                                        >  2:4784484/1‑150 (MQ=32)
          gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                       >  1:4676989/1‑150 (MQ=32)
          gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                       <  1:2849270/150‑1 (MQ=32)
          gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                       <  2:5027218/150‑1 (MQ=32)
           acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAgg                                                                                                                                      <  2:1337180/150‑1 (MQ=255)
            cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                     <  1:5106352/150‑1 (MQ=32)
            cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                     <  2:7119718/150‑1 (MQ=255)
            cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                     <  2:3503175/150‑1 (MQ=255)
            cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                     <  1:5014414/150‑1 (MQ=11)
            cgacCATAGAGCATTGGAAACACCTGATCCCATCCCGATCTCAGTTGTGAAACGATGAATGGCCGTTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTGATTTCTCAAAACGCCTACCTGCATGTGCAGATAggg                                                                                                                                     <  2:3059856/150‑1 (MQ=11)
             gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAgggg                                                                                                                                    <  2:4394804/150‑1 (MQ=32)
             gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAgggg                                                                                                                                    <  2:1166527/150‑1 (MQ=255)
              acCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGt                                                                                                                                   <  1:409119/150‑1 (MQ=32)
                cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGttt                                                                                                                                 <  2:6146284/150‑1 (MQ=32)
                 aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtttt                                                                                                                                <  2:5346491/150‑1 (MQ=32)
                 aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtttt                                                                                                                                <  1:8432402/150‑1 (MQ=32)
                 aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtttt                                                                                                                                <  2:266864/150‑1 (MQ=32)
                      gCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCtt                                                                                                                           <  2:3176885/150‑1 (MQ=32)
                      gCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCtt                                                                                                                           <  1:1163380/150‑1 (MQ=32)
                       gaTTTGAACCTCGTGATCCCATCCAGAACTCAGTAGTGAAACGATGCTTCGCCGATGGTAGTGTGGTGTTTCCCCATGTGAGAGTAGGTCATCGTCTAGATTCATTTCGCTAAACCCCTCTCGGCACATGCAGGTTGGGGTTTTGTCttt                                                                                                                          <  2:194659/149‑1 (MQ=11)
                        aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGACACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTGGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTa                                                                                                                         <  2:586692/150‑1 (MQ=17)
                        aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTa                                                                                                                         <  2:2761472/150‑1 (MQ=32)
                         ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        <  1:7655742/150‑1 (MQ=32)
                         ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        <  1:5819117/150‑1 (MQ=32)
                         ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        <  1:418265/150‑1 (MQ=32)
                         ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        <  1:2200919/150‑1 (MQ=32)
                                 accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTCGGGGTTTTGTCTTTAAGTAGGAAc                                                                                                                >  1:1859353/1‑150 (MQ=25)
                                 accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAc                                                                                                                <  1:621747/150‑1 (MQ=35)
                                      atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTGTTGTCTTTAAGTAGGAACTACag                                                                                                           >  2:8258283/1‑150 (MQ=255)
                                       tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                          <  1:7276781/150‑1 (MQ=40)
                                       tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                          <  2:372376/150‑1 (MQ=40)
                                         ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAgaga                                                                                                        >  2:8144621/1‑150 (MQ=255)
                                            tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTc                                                                                                     <  2:7436173/150‑1 (MQ=255)
                                             cccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCACAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCg                                                                                                    <  2:6180982/150‑1 (MQ=255)
                                                  aCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGc                                                                                               <  2:2313884/150‑1 (MQ=255)
                                                   cTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTa                                                                                                                     >  2:2470957/1‑127 (MQ=32)
                                                    tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCAc                                                                                             >  1:822788/1‑150 (MQ=255)
                                                          gtGAAACGATGCATCGCCTATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTa                                                                                       <  1:925881/150‑1 (MQ=255)
                                                          gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTa                                                                                       <  2:6450969/150‑1 (MQ=255)
                                                          gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCATATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTccga                                                                                       >  2:2618498/1‑148 (MQ=255)
                                                           tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACGGGCAGGGCGGGGTTGTGTCTTGAAGTAGGAACTACAGAGATGCGCAGGCACggcggag                                                                                      >  2:6608451/1‑144 (MQ=17)
                                                           tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAg                                                                                      <  2:7397238/150‑1 (MQ=255)
                                                            gAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAgg                                                                                     <  2:42894/150‑1 (MQ=255)
                                                              aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGa                                                                                   <  1:4388804/150‑1 (MQ=255)
                                                              aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGa                                                                                   <  1:1621594/150‑1 (MQ=255)
                                                                 gATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                <  1:4899470/150‑1 (MQ=255)
                                                                 gATGCATCGCCGATGGTAATGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                <  2:3973869/150‑1 (MQ=255)
                                                                      aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCCCAGGCCGGTAGGGGTTTTGTCTTTAAGTCGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg                                                                           >  2:8654096/1‑150 (MQ=255)
                                                                      aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCCCGTCCCTAGGACGTGcctg                                                                           <  1:6532515/150‑3 (MQ=17)
                                                                      aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg                                                                           <  2:4709520/150‑1 (MQ=255)
                                                                            gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacaca                                                                     <  1:7036650/150‑1 (MQ=255)
                                                                                gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAt                                                                 <  1:1234576/150‑1 (MQ=255)
                                                                                 tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAtt                                                                >  1:1323051/1‑150 (MQ=255)
                                                                                     gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                            <  1:4506490/150‑1 (MQ=255)
                                                                                      tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTg                                                           >  2:5931907/1‑150 (MQ=255)
                                                                                       ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                          <  1:7812406/150‑1 (MQ=255)
                                                                                       ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                          >  2:5156162/1‑150 (MQ=255)
                                                                                       ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                          <  2:8379413/150‑1 (MQ=255)
                                                                                       ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGACACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                          <  1:8674109/150‑1 (MQ=255)
                                                                                           tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacgac                                                      <  2:6967924/150‑1 (MQ=255)
                                                                                               cccATGTGAGAGTAGGTGATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                                  <  2:8481191/150‑1 (MQ=255)
                                                                                               cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                                  <  1:5567138/150‑1 (MQ=255)
                                                                                                ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATag                                                 >  2:5680570/1‑150 (MQ=255)
                                                                                                ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATag                                                 >  2:1431154/1‑150 (MQ=255)
                                                                                                 cATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATaga                                                >  1:8641556/1‑150 (MQ=255)
                                                                                                      gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                           <  2:4875204/150‑1 (MQ=255)
                                                                                                       agagTAGGTCATCGTCAAGATTCAGTTCGCAAAACCCCTATCTGCACATGCCGGGAGGGGGTTTGTCTTTAAGTAGGAACTACTGAGATTCGCAGGCACGTCGTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCGTAGAGCATTg                                          >  1:3987600/1‑150 (MQ=25)
                                                                                                        gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTgg                                         <  2:4676989/150‑1 (MQ=255)
                                                                                                         agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGa                                        <  2:6164953/150‑1 (MQ=255)
                                                                                                         agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGa                                        <  1:3145929/150‑1 (MQ=255)
                                                                                                               tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccacc                                  >  1:7951739/1‑150 (MQ=255)
                                                                                                                    gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGACTTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                             >  2:4493578/1‑150 (MQ=255)
                                                                                                                    gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                             >  1:6911391/1‑150 (MQ=255)
                                                                                                                    gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGAAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                             >  1:2907518/1‑150 (MQ=255)
                                                                                                                     tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCTGGCCCGTCCTAGGGACGGGCCa                                                                            >  2:8058901/1‑102 (MQ=255)
                                                                                                                     tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCTGGCCCGTCCTAGGGACGGGCCa                                                                            <  1:8058901/102‑1 (MQ=255)
                                                                                                                     tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                            >  1:5701865/1‑150 (MQ=255)
                                                                                                                     tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                            >  2:1583737/1‑150 (MQ=255)
                                                                                                                       aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatccca                          >  1:5097281/1‑150 (MQ=255)
                                                                                                                         gATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatc                        <  2:1859353/150‑1 (MQ=255)
                                                                                                                            tCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCg                     >  2:3815546/1‑150 (MQ=255)
                                                                                                                            tCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCg                     >  1:3380810/1‑150 (MQ=255)
                                                                                                                             ttaaagacAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTGCCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGa                    >  1:1008046/8‑150 (MQ=37)
                                                                                                                               tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAAc                  >  1:8536633/1‑150 (MQ=255)
                                                                                                                                  gacAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTca               >  2:6548756/3‑150 (MQ=255)
                                                                                                                                  cGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTca               <  2:5791530/150‑1 (MQ=255)
                                                                                                                                   gCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag              >  2:766218/1‑150 (MQ=255)
                                                                                                                                    cAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagc             >  1:7472345/1‑150 (MQ=255)
                                                                                                                                          cccTATCTGCACCTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaa       >  2:6368082/1‑150 (MQ=255)
                                                                                                                                            cTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGa                                        >  1:6003013/1‑115 (MQ=255)
                                                                                                                                            cTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc     >  1:6669985/1‑150 (MQ=255)
                                                                                                                                             tATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACg    >  2:6881303/1‑150 (MQ=255)
                                                                                                                                               tCTTCACATGCAGGTAGGGGTTTTGTCTTTAATTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt  <  2:5482085/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181662‑181954

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