Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,818 G→A 13.2% 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,8180GA13.2% 172.8 / 14.6 68intergenic (+32/‑74)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (5/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATG  >  NC_002947/181672‑181966
                                                                                                                                                  |                                                                                                                                                       
gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt                                                                                                                                                        <  2:200884/148‑1 (MQ=32)
 acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                                      >  6:316967/1‑149 (MQ=32)
   gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCGTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAg                                                                                                                                                      >  8:157861/1‑147 (MQ=11)
   gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCGTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAg                                                                                                                                                      <  7:157861/147‑1 (MQ=11)
   gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                                    >  2:319644/1‑149 (MQ=32)
   gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                                    >  1:25317/1‑149 (MQ=32)
   gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                                    >  6:333403/1‑149 (MQ=11)
     ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAGGGGt                                                                                                                                                  <  7:165346/149‑1 (MQ=255)
       aTAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGttt                                                                                                                                                <  2:103929/149‑1 (MQ=32)
          gagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGt                                                                                                                                             <  5:316967/149‑1 (MQ=32)
               ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCttt                                                                                                                                         >  2:250024/1‑148 (MQ=32)
               ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTa                                                                                                                                        <  1:319644/149‑1 (MQ=32)
                 ggAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAg                                                                                                                                      >  6:332006/1‑149 (MQ=32)
                   aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCAGCTTTGATATTAGGTAATCGTAAAAATTAATTTAGCAAAACCCCAATATGCACATGAATGTAGGGGTTTTGTCTTTAAGt                                                                                                                                     >  1:262069/1‑148 (MQ=11)
                        ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAc                                                                                                                               <  3:49209/149‑1 (MQ=35)
                           gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑a                                                                                                                           >  7:289395/1‑148 (MQ=37)
                             tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑a                                                                                                                           >  1:12482/1‑146 (MQ=37)
                                  tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGCGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGAtt                                                                                                                  >  2:100321/1‑149 (MQ=255)
                                      gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCg                                                                                                                >  8:293501/1‑147 (MQ=255)
                                       aaCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAg                                                                                                             <  5:332006/149‑1 (MQ=255)
                                             gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAGGGGt                                                                                                                                                  >  1:196568/1‑109 (MQ=11)
                                             gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAGGGGt                                                                                                                                                  <  2:196568/109‑1 (MQ=11)
                                               agtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTc                                                                                                      >  5:21827/1‑148 (MQ=255)
                                                     aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCAc                                                                                                         >  1:109739/1‑139 (MQ=255)
                                                     aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCAc                                                                                                         <  2:109739/139‑1 (MQ=255)
                                                           cATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGc                                                                                          >  8:214016/1‑148 (MQ=255)
                                                            aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGc                                                                                          >  5:98285/1‑147 (MQ=255)
                                                            aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGCAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCa                                                                                        <  7:293501/149‑1 (MQ=255)
                                                                ccGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGca                                                                                     >  2:259404/1‑148 (MQ=255)
                                                                   aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacaca                                                                                 <  2:142823/149‑1 (MQ=255)
                                                                      gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGaa                                                                              <  2:12482/149‑1 (MQ=255)
                                                                      gTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGa                                                                               <  2:25317/148‑1 (MQ=255)
                                                                        aGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAtt                                                                            >  8:195253/1‑149 (MQ=255)
                                                                            tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGc                                                                                                               <  4:187908/110‑1 (MQ=255)
                                                                            tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGc                                                                                                               >  3:187908/1‑110 (MQ=255)
                                                                              tggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTATTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                                      <  4:278061/148‑1 (MQ=255)
                                                                                   tCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacgac                                                                  <  6:98285/148‑1 (MQ=255)
                                                                                     cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ATCAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                                <  1:250024/149‑1 (MQ=255)
                                                                                      ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                                >  6:169483/1‑147 (MQ=255)
                                                                                       cATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATag                                                             <  7:214016/149‑1 (MQ=255)
                                                                                       cATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGACAGCACACAGAATTTCTTGACGACCATag                                                             <  6:294912/149‑1 (MQ=255)
                                                                                        aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATaga                                                            <  1:259404/149‑1 (MQ=255)
                                                                                          gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATagag                                                           >  2:221990/1‑148 (MQ=255)
                                                                                                gTAGGTCATCGTCAATATTCATTTCGCAAAACCCCTATCTGCACATTCATGTAGGGGTTTTGTCTTTAAGTAGTAACTTCCCAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                       <  1:100321/148‑1 (MQ=25)
                                                                                                       aTCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccacc                                              <  2:262069/148‑1 (MQ=255)
                                                                                                       aTCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACAGGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCt                                             >  4:152450/1‑149 (MQ=255)
                                                                                                          gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGc                                                                                      <  6:47060/105‑1 (MQ=255)
                                                                                                          gTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCAAGGGACGGGCCAGc                                                                                      >  5:47060/1‑105 (MQ=255)
                                                                                                           tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTc                                                                          >  6:4914/1‑116 (MQ=255)
                                                                                                           tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTc                                                                          <  5:4914/116‑1 (MQ=255)
                                                                                                            cAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                                        >  7:55552/1‑149 (MQ=255)
                                                                                                              aGATTCATATCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGTAACT‑‑ACAGAGATTCGCAGGCACTTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatccca                                      <  8:289395/149‑1 (MQ=255)
                                                                                                                aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                                         <  3:309118/144‑1 (MQ=255)
                                                                                                                aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatc                                         >  4:309118/1‑144 (MQ=255)
                                                                                                                 ttCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatc                                    <  7:195253/148‑1 (MQ=255)
                                                                                                                   ttaaagacAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACCTCCTAGGGCCGGGCCAGCACACAGAATTTCATGACGACCATAGACCATTGGAACCACCTGaacacatccc                                  >  6:320723/8‑148 (MQ=37)
                                                                                                                    aTTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTg                                                      >  5:241302/1‑127 (MQ=255)
                                                                                                                    aTTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCg                                 <  1:221990/148‑1 (MQ=255)
                                                                                                                    aTTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACACAATTTCTTGACGACCATAGAGCATTg                                                      <  6:241302/127‑1 (MQ=255)
                                                                                                                     tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                               >  6:48642/1‑149 (MQ=255)
                                                                                                                           aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag                          >  2:219359/1‑148 (MQ=255)
                                                                                                                            aaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcagca                        <  6:21827/149‑1 (MQ=255)
                                                                                                                               ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGa                                           <  6:53638/127‑1 (MQ=255)
                                                                                                                               ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGa                                           >  5:53638/1‑127 (MQ=255)
                                                                                                                                cccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATACCATCCCGAACTCAGCAGTg                     >  5:313340/1‑148 (MQ=255)
                                                                                                                                 ccTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGa                    >  6:64597/1‑148 (MQ=255)
                                                                                                                                            aTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                                          <  3:180587/115‑1 (MQ=255)
                                                                                                                                            aTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                                          >  4:180587/1‑115 (MQ=255)
                                                                                                                                                 ggTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAt   <  5:48642/149‑1 (MQ=255)
                                                                                                                                                  gTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATg  >  2:133093/1‑149 (MQ=255)
                                                                                                                                                  gTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATg  >  4:253334/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATG  >  NC_002947/181672‑181966

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