Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,811 C→T 8.5% intergenic (+25/‑81) 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,8110CT 8.5% 189.4 / 6.3 70intergenic (+25/‑81)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base T (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 = 7.09e-01

GAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181663‑181954
                                                                                                                                                       |                                                                                                                                                  
gAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGcac                                                                                                                                                    <  2:241199/149‑1 (MQ=31)
     tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAAGAGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATg                                                                                                                                                 >  1:294450/1‑148 (MQ=12)
         gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTTTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGt                                                                                                                                            <  5:296988/148‑1 (MQ=17)
            gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATg                                                                                                                                                 <  3:69029/140‑1 (MQ=32)
            gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATg                                                                                                                                                 >  4:69029/1‑140 (MQ=32)
            gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAgg‑‑‑g                                                                                                                                     <  8:167752/149‑1 (MQ=32)
              ccATAGAGCATGGTAACCACCTGATCCCATCACGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTTTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGt                                                                                                                                   <  1:157070/149‑1 (MQ=255)
                  agagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTg                                                                                                                               <  2:251867/149‑1 (MQ=32)
                     gCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCCGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTc                                                                                                                             <  1:284462/149‑1 (MQ=17)
                     gCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCt                                                                                                                            <  2:110548/149‑1 (MQ=32)
                       aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCtt                                                                                                                           >  6:56452/1‑148 (MQ=32)
                        ttGGAACCACCTGATCCCATCCCGAACGCAGTAGTGAAA‑CGATGCATAGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCttt                                                                                                                          >  1:207257/1‑148 (MQ=17)
                          ggAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTaa                                                                                                                        <  6:91099/148‑1 (MQ=32)
                           gAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAg                                                                                                                       >  6:145563/1‑148 (MQ=32)
                            aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTa                                                                                                                     <  7:261551/149‑1 (MQ=32)
                               caccTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAgg                                                                                                                   >  5:163727/1‑148 (MQ=32)
                               caccTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAgg                                                                                                                   >  4:214587/1‑148 (MQ=32)
                               caccTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGa                                                                                                                  <  7:84101/149‑1 (MQ=33)
                                accTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGaa                                                                                                                 <  7:2884/149‑1 (MQ=34)
                                 ccTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAAc                                                                                                                >  1:195717/1‑149 (MQ=35)
                                 ccTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAAc                                                                                                                >  5:6411/1‑149 (MQ=35)
                                        ccatcccGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACagag                                                                                                         <  3:285928/149‑1 (MQ=40)
                                        ccatcccGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACa                                                                                                            >  2:368911/1‑146 (MQ=38)
                                        ccatcccGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTACC‑TGCATGTGCAGATAgg‑‑‑g                                                                                                                                     <  7:92421/121‑1 (MQ=255)
                                        ccatcccGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTACC‑TGCATGTGCAGATAgg‑‑‑g                                                                                                                                     >  8:92421/1‑121 (MQ=255)
                                         catcccGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAg                                                                                                                    >  4:334/1‑137 (MQ=32)
                                         catcccGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAg                                                                                                                    <  3:334/137‑1 (MQ=32)
                                              cGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCg                                                                                                    >  1:4969/1‑148 (MQ=255)
                                               gAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGTTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGTAACTACAGAGATTCGCa                                                                                                  <  3:214587/149‑1 (MQ=255)
                                                   tCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCa                                                                                              >  6:110062/1‑149 (MQ=255)
                                                       tagtGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTc                                                                                          <  1:18424/149‑1 (MQ=255)
                                                        agtGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTcc                                                                                         >  5:255513/1‑149 (MQ=255)
                                                            aaa‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCa                                                                                              <  7:283958/140‑1 (MQ=255)
                                                            aaa‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCa                                                                                              >  8:283958/1‑140 (MQ=255)
                                                                  aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                <  2:4969/149‑1 (MQ=255)
                                                                  aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                                <  5:110062/149‑1 (MQ=255)
                                                                     cATCGCCGATGGTAGTGTGGGGTATCACCATGTGAGAGTACGTCAT‑CGTCAAGATTCATTTCGCAAAACCCATATC‑TGCACATGCAGGAAGG‑‑‑GGTTTTCTGTTTAAGTAg                                                                                                                    >  2:35657/1‑110 (MQ=11)
                                                                      aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCCCGTCCCTAGGACGTGcct                                                                            <  6:6411/149‑2 (MQ=17)
                                                                       tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg                                                                           >  2:19361/1‑149 (MQ=255)
                                                                             aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacaca                                                                     <  5:45076/149‑1 (MQ=255)
                                                                                gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGaa                                                                  <  2:195717/149‑1 (MQ=255)
                                                                                   gtgtGGGGTCTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAtt                                                                <  5:230260/148‑1 (MQ=255)
                                                                                      tgGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt                                                             <  5:145563/148‑1 (MQ=255)
                                                                                        gggTTTGCCCATGTGCCTTTAGGGCAT‑CGTCAACATTCATTTCCCAAAACCCCTTTCTTCCACATGCAGGTAGGCTCCCTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTc                                                              <  1:116707/149‑1 (MQ=255)
                                                                                         ggTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgac                                                         <  6:163727/149‑1 (MQ=255)
                                                                                                ccATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                    >  4:113365/1‑147 (MQ=255)
                                                                                                ccATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATa                                                  >  7:159040/1‑149 (MQ=255)
                                                                                                       agagTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                           <  2:207257/149‑1 (MQ=255)
                                                                                                       agagTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                           <  1:368911/149‑1 (MQ=255)
                                                                                                         agTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                    >  8:304841/1‑138 (MQ=255)
                                                                                                         agTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                    <  7:304841/138‑1 (MQ=255)
                                                                                                          gTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTgg                                         <  6:255513/148‑1 (MQ=255)
                                                                                                          gTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGa                                        <  5:56452/149‑1 (MQ=255)
                                                                                                                     gTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                              >  6:257030/1‑149 (MQ=255)
                                                                                                                      tCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                             >  6:57026/1‑149 (MQ=255)
                                                                                                                                tttCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                   >  5:175730/1‑149 (MQ=255)
                                                                                                                                tttCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                   >  3:120857/1‑149 (MQ=255)
                                                                                                                                  tCGCAAAACCCCAATC‑TGCACATGCAGGTAGT‑‑‑GGATTTGCCGTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTATTGACGACCATAGAGCATTGGAACCTCCGGATCCCATCCCGAACt                 >  3:39365/1‑149 (MQ=37)
                                                                                                                                  agacAAAACCCCTACC‑TGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                 <  7:219585/146‑1 (MQ=37)
                                                                                                                                  agacAAAACCCCTACC‑TGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                 >  8:219585/4‑149 (MQ=37)
                                                                                                                                    gCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt                 >  6:191077/1‑147 (MQ=255)
                                                                                                                                      aaaaCCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag              >  3:97950/1‑148 (MQ=255)
                                                                                                                                       aaaCCCCAAAA‑TCCACATGCGGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAAGTACAGAGATACCCCGGCACGACCTAGGACCGGCCCAGCACACAGAATTTCTTGACGACCTTAACGCATTGTAACCACCTCATCCAAGCCCGaaagaccac            >  2:107912/1‑142 (MQ=255)
                                                                                                                                           cccTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  2:121374/1‑148 (MQ=255)
                                                                                                                                           cccTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  7:260100/1‑148 (MQ=255)
                                                                                                                                           cccTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  6:188622/1‑148 (MQ=255)
                                                                                                                                           cccTACC‑TGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  3:211416/1‑148 (MQ=255)
                                                                                                                                           cccTACC‑TGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg         >  8:21918/1‑148 (MQ=255)
                                                                                                                                             caaTC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa      >  5:90275/3‑149 (MQ=255)
                                                                                                                                             cTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa      >  7:271746/1‑149 (MQ=255)
                                                                                                                                             cTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa      >  1:105929/1‑149 (MQ=255)
                                                                                                                                                     cacaTGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt  <  6:175730/146‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
GAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181663‑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: ‑