Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,805 T→C 8.4% 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 8.4% 187.5 / 7.3 72intergenic (+19/‑87)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (34/31);  new base C (4/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAA  >  NC_002947/181662‑181950
                                                                                                                                                 |                                                                                                                                                     
aGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGca                                                                                                                                                 >  1:241199/1‑149 (MQ=11)
 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                                                                                                            >  5:6411/1‑149 (MQ=35)
                                  ccTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAAc                                                                                                            >  1:195717/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                                                                                                                                 >  8:92421/1‑121 (MQ=255)
                                         ccatcccGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTACC‑TGCATGTGCAGATAgg‑‑‑g                                                                                                                                 <  7:92421/121‑1 (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                                                                            <  5:110062/149‑1 (MQ=255)
                                                                   aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACgg                                                                            <  2:4969/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                                       <  1:368911/149‑1 (MQ=255)
                                                                                                        agagTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                       <  2:207257/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     >  7:260100/1‑148 (MQ=255)
                                                                                                                                            cccTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg     >  2:121374/1‑148 (MQ=255)
                                                                                                                                            cccTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg     >  6:188622/1‑148 (MQ=255)
                                                                                                                                            cccTACC‑TGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg     >  8:21918/1‑148 (MQ=255)
                                                                                                                                            cccTACC‑TGCATGTGCAGATAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg     >  3:211416/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)
                                                                                                                                                 |                                                                                                                                                     
AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAA‑CGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCAT‑CGTCAAGATTCATTTCGCAAAACCCCTATC‑TGCACATGCAGGTAGG‑‑‑GGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAA  >  NC_002947/181662‑181950

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: ‑