Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,812 A→G 8.6% intergenic (+26/‑80) 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,8120AG 8.6% 188.1 / 7.2 70intergenic (+26/‑80)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base G (4/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

TCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181668‑181954
                                                                                                                                                   |                                                                                                                                                 
tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAAGAGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATg                                                                                                                                                 >  1:294450/1‑148 (MQ=12)
    gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTTTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt                                                                                                                                            <  5:296988/148‑1 (MQ=17)
       gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATg                                                                                                                                                 <  3:69029/140‑1 (MQ=32)
       gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATg                                                                                                                                                 >  4:69029/1‑140 (MQ=32)
       gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAgg‑‑‑g                                                                                                                                     <  8:167752/149‑1 (MQ=32)
         ccATAGAGCATGGTAACCACCTGATCCCATCACGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTTTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGt                                                                                                                                   <  1:157070/149‑1 (MQ=255)
             agagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTg                                                                                                                               <  2:251867/149‑1 (MQ=32)
                gCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTc                                                                                                                             <  1:284462/149‑1 (MQ=17)
                gCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCt                                                                                                                            <  2:110548/149‑1 (MQ=32)
                  aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCtt                                                                                                                           >  6:56452/1‑148 (MQ=32)
                   ttGGAACCACCTGATCCCATCCCGAACGCAGTAGTGAAACGATGCATAGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCttt                                                                                                                          >  1:207257/1‑148 (MQ=17)
                     ggAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTaa                                                                                                                        <  6:91099/148‑1 (MQ=32)
                      gAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAg                                                                                                                       >  6:145563/1‑148 (MQ=32)
                       aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTa                                                                                                                     <  7:261551/149‑1 (MQ=32)
                          caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAgg                                                                                                                   >  4:214587/1‑148 (MQ=32)
                          caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAgg                                                                                                                   >  5:163727/1‑148 (MQ=32)
                          caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGa                                                                                                                  <  7:84101/149‑1 (MQ=33)
                           accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGaa                                                                                                                 <  7:2884/149‑1 (MQ=34)
                            ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAAc                                                                                                                >  1:195717/1‑149 (MQ=35)
                            ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAAc                                                                                                                >  5:6411/1‑149 (MQ=35)
                                   ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACagag                                                                                                         <  3:285928/149‑1 (MQ=40)
                                   ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACa                                                                                                            >  2:368911/1‑146 (MQ=38)
                                   ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAgg‑‑‑g                                                                                                                                     <  7:92421/121‑1 (MQ=255)
                                   ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAgg‑‑‑g                                                                                                                                     >  8:92421/1‑121 (MQ=255)
                                    catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAg                                                                                                                    >  4:334/1‑137 (MQ=32)
                                    catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAg                                                                                                                    <  3:334/137‑1 (MQ=32)
                                         cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCg                                                                                                    >  1:4969/1‑148 (MQ=255)
                                          gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGTTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGTAACTACAGAGATTCGCa                                                                                                  <  3:214587/149‑1 (MQ=255)
                                              tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCa                                                                                              >  6:110062/1‑149 (MQ=255)
                                                  tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTc                                                                                          <  1:18424/149‑1 (MQ=255)
                                                   agtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTcc                                                                                         >  5:255513/1‑149 (MQ=255)
                                                       aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCa                                                                                              <  7:283958/140‑1 (MQ=255)
                                                       aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCa                                                                                              >  8:283958/1‑140 (MQ=255)
                                                             aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                <  2:4969/149‑1 (MQ=255)
                                                             aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                <  5:110062/149‑1 (MQ=255)
                                                                cATCGCCGATGGTAGTGTGGGGTATCACCATGTGAGAGTACGTCATCGTCAAGATTCATTTCGCAAAACCCATATCTGCACATGCAGGAAGG‑‑‑GGTTTTCTGTTTAAGTAg                                                                                                                    >  2:35657/1‑110 (MQ=11)
                                                                 aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCCCGTCCCTAGGACGTGcct                                                                            <  6:6411/149‑2 (MQ=17)
                                                                  tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg                                                                           >  2:19361/1‑149 (MQ=255)
                                                                        aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacaca                                                                     <  5:45076/149‑1 (MQ=255)
                                                                           gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGaa                                                                  <  2:195717/149‑1 (MQ=255)
                                                                              gtgtGGGGTCTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAtt                                                                <  5:230260/148‑1 (MQ=255)
                                                                                 tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt                                                             <  5:145563/148‑1 (MQ=255)
                                                                                   gggTTTGCCCATGTGCCTTTAGGGCATCGTCAACATTCATTTCCCAAAACCCCTTTCTTCCACATGCAGGTAGGCTCCCTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTc                                                              <  1:116707/149‑1 (MQ=255)
                                                                                    ggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgac                                                         <  6:163727/149‑1 (MQ=255)
                                                                                           ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                    >  4:113365/1‑147 (MQ=255)
                                                                                           ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                                  >  7:159040/1‑149 (MQ=255)
                                                                                                  agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                           <  2:207257/149‑1 (MQ=255)
                                                                                                  agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                           <  1:368911/149‑1 (MQ=255)
                                                                                                    agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                    >  8:304841/1‑138 (MQ=255)
                                                                                                    agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                    <  7:304841/138‑1 (MQ=255)
                                                                                                     gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTgg                                         <  6:255513/148‑1 (MQ=255)
                                                                                                     gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGa                                        <  5:56452/149‑1 (MQ=255)
                                                                                                                gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                              >  6:257030/1‑149 (MQ=255)
                                                                                                                 tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                             >  6:57026/1‑149 (MQ=255)
                                                                                                                           tttCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                   >  5:175730/1‑149 (MQ=255)
                                                                                                                           tttCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                   >  3:120857/1‑149 (MQ=255)
                                                                                                                             tCGCAAAACCCCAATCTGCACATGCAGGTAGT‑‑‑GGATTTGCCGTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTATTGACGACCATAGAGCATTGGAACCTCCGGATCCCATCCCGAACt                 >  3:39365/1‑149 (MQ=37)
                                                                                                                             agacAAAACCCCTACCTGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                 <  7:219585/146‑1 (MQ=37)
                                                                                                                             agacAAAACCCCTACCTGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                 >  8:219585/4‑149 (MQ=37)
                                                                                                                               gCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                 >  6:191077/1‑147 (MQ=255)
                                                                                                                                 aaaaCCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag              >  3:97950/1‑148 (MQ=255)
                                                                                                                                  aaaCCCCAAAATCCACATGCGGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAAGTACAGAGATACCCCGGCACGACCTAGGACCGGCCCAGCACACAGAATTTCTTGACGACCTTAACGCATTGTAACCACCTCATCCAAGCCCGaaagaccac            >  2:107912/1‑142 (MQ=255)
                                                                                                                                      cccTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  2:121374/1‑148 (MQ=255)
                                                                                                                                      cccTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  7:260100/1‑148 (MQ=255)
                                                                                                                                      cccTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  6:188622/1‑148 (MQ=255)
                                                                                                                                      cccTACCTGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  3:211416/1‑148 (MQ=255)
                                                                                                                                      cccTACCTGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  8:21918/1‑148 (MQ=255)
                                                                                                                                        caaTCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa      >  5:90275/3‑149 (MQ=255)
                                                                                                                                        cTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa      >  7:271746/1‑149 (MQ=255)
                                                                                                                                        cTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa      >  1:105929/1‑149 (MQ=255)
                                                                                                                                                cacaTGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt  <  6:175730/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181668‑181954

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