Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,811 C→T 13.4% 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,8110CT13.4% 168.0 / 15.6 67intergenic (+25/‑81)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base T (5/4);  total (34/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 = 6.29e-01

TCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGA  >  NC_002947/181668‑181948
                                                                                                                                               |                                                                                                                                            
tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGc                                                                                                                                          <  4:166224/148‑1 (MQ=32)
    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)
                                   ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGcac                                                                                                                                              >  6:298186/1‑109 (MQ=31)
                                   ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGcac                                                                                                                                              <  5:298186/109‑1 (MQ=31)
                                      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                                                                                                                                <  2:196568/109‑1 (MQ=11)
                                                 gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAGGGGt                                                                                                                                >  1:196568/1‑109 (MQ=11)
                                                   agtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTc                                                                                    >  5:21827/1‑148 (MQ=255)
                                                         aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCAc                                                                                       <  2:109739/139‑1 (MQ=255)
                                                         aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCAc                                                                                       >  1:109739/1‑139 (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                                                        <  5:4914/116‑1 (MQ=255)
                                                                                                               tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTc                                                        >  6:4914/1‑116 (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)
                                                                                                                                               |                                                                                                                                            
TCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACT‑‑ACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGA  >  NC_002947/181668‑181948

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: