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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947181,8180GA 8.4% 194.0 / 7.2 71intergenic (+32/‑74)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (4/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAT  >  NC_002947/181672‑181965
                                                                                                                                                    |                                                                                                                                                      
gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTTTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt                                                                                                                                                       <  5:296988/148‑1 (MQ=17)
   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                                                                                                                               <  3:334/137‑1 (MQ=32)
                                catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAg                                                                                                                               >  4:334/1‑137 (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                                                                                                         >  8:283958/1‑140 (MQ=255)
                                                   aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCa                                                                                                         <  7:283958/140‑1 (MQ=255)
                                                        aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                           <  5:110062/149‑1 (MQ=255)
                                                        aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                           <  2:4969/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                                                      <  1:368911/149‑1 (MQ=255)
                                                                                             agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                      <  2:207257/149‑1 (MQ=255)
                                                                                               agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                               <  7:304841/138‑1 (MQ=255)
                                                                                               agTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                               >  8:304841/1‑138 (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                              >  3:120857/1‑149 (MQ=255)
                                                                                                                      tttCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                              >  5:175730/1‑149 (MQ=255)
                                                                                                                        tCGCAAAACCCCAATCTGCACATGCAGGTAGT‑‑‑GGATTTGCCGTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTATTGACGACCATAGAGCATTGGAACCTCCGGATCCCATCCCGAACt                            >  3:39365/1‑149 (MQ=37)
                                                                                                                        agacAAAACCCCTACCTGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                            >  8:219585/4‑149 (MQ=37)
                                                                                                                        agacAAAACCCCTACCTGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                            <  7:219585/146‑1 (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                    >  8:21918/1‑148 (MQ=255)
                                                                                                                                 cccTACCTGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg                    >  3:211416/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)
                                                                                                                                                gCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGcc     >  3:161171/1‑149 (MQ=255)
                                                                                                                                                 cAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGcc     <  5:257030/148‑1 (MQ=255)
                                                                                                                                                  aGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGa   <  1:19361/149‑1 (MQ=255)
                                                                                                                                                   ggTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAt  <  8:260100/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAT  >  NC_002947/181672‑181965

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: