Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,471,402 G→T 8.8% R478R (CGG→AGG)  alg8 ← glycosyltransferase Alg8

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,471,4020GT 8.8% 249.4 / 8.8 79R478R (CGG→AGG) alg8glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base T (3/4);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GGCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAT  >  NC_002947/1471254‑1471546
                                                                                                                                                    |                                                                                                                                                
ggCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCcgc                                                                                                                                                 >  1:956865/1‑150 (MQ=255)
   gcTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAg                                                                                                                                              >  1:103810/1‑150 (MQ=255)
       ggCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGAcc                                                                                                                                          <  2:956865/150‑1 (MQ=255)
        gCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCa                                                                                                                                         <  1:640218/150‑1 (MQ=255)
           tCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatga                                                                                                                                      <  1:104500/150‑4 (MQ=255)
           tCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgt                                                                                                                                      <  2:103810/150‑1 (MQ=255)
            cGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatgac                                                                                                                                     <  1:1273544/150‑5 (MQ=255)
                    tGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCa                                                                                                                             >  1:856860/1‑150 (MQ=255)
                      cACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATGCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                           >  2:1415537/1‑150 (MQ=255)
                      cACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatgacc                                                                                                                                    >  2:302709/1‑136 (MQ=255)
                      cACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatgacc                                                                                                                                    <  1:302709/141‑6 (MQ=255)
                         gACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                           >  2:1004777/1‑147 (MQ=255)
                         gACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                           <  1:1004777/147‑1 (MQ=255)
                            aTTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGtt                                                                                                                                   <  1:1141185/136‑1 (MQ=255)
                            aTTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGtt                                                                                                                                   >  2:1141185/1‑136 (MQ=255)
                                 cGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTg                                                                                                                <  1:1058359/150‑1 (MQ=255)
                                      aCGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAg                                                                                                           >  2:948174/1‑150 (MQ=255)
                                            gTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAGTGTTGAAAGCTGGCGAGGTCAcc                                                                                                     >  2:1242389/1‑149 (MQ=255)
                                              aTTCATGGTGGTTTGTTCATGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGTGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTAACGCt                                                                                                   >  1:336783/1‑150 (MQ=255)
                                               ttCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                              >  1:468222/1‑138 (MQ=255)
                                               ttCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                              <  2:468222/138‑1 (MQ=255)
                                                   tggtggTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgtgt                                                                                                                                    <  2:910346/112‑1 (MQ=255)
                                                   tggtggTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgtgt                                                                                                                                    >  1:910346/1‑112 (MQ=255)
                                                    ggtggtCTGTCCCAGTTCGACCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAggg                                                                                             <  2:443256/150‑1 (MQ=255)
                                                     gtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGc                                                                                            <  2:565179/150‑1 (MQ=255)
                                                     gtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGc                                                                                            >  2:686666/1‑150 (MQ=255)
                                                            gTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                           <  2:1212208/112‑1 (MQ=255)
                                                            gTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                           >  1:1212208/1‑112 (MQ=255)
                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCatga                                                                                                                                      <  2:950838/100‑4 (MQ=255)
                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCatga                                                                                                                                      >  1:950838/1‑97 (MQ=255)
                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAgg                                                                                                                                       >  2:497884/1‑99 (MQ=255)
                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAgg                                                                                                                                       <  1:497884/99‑1 (MQ=255)
                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCt                                                                                    <  2:856860/150‑1 (MQ=255)
                                                                 tGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCCGCCGAGAAGGTCATGGTCCGCGAGGACCAGGGGTTGAACCATTGTTGAACGCTGGCGAGGCTACGCTGGCGGGCGGTGGTCTGAcg                                                                                >  2:661950/1‑150 (MQ=255)
                                                                  gTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGATGGCAGTGGTCTGAcgc                                                                               <  2:1002345/150‑1 (MQ=255)
                                                                       atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAggg                                                                                             <  1:857970/131‑1 (MQ=255)
                                                                       atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAggg                                                                                             >  2:857970/1‑131 (MQ=255)
                                                                       atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCa                                                                          <  2:142986/150‑1 (MQ=255)
                                                                        tccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGtt                                                                                                                                   >  1:653820/1‑92 (MQ=255)
                                                                        tccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGtt                                                                                                                                   <  2:653820/92‑1 (MQ=255)
                                                                                       gCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGa                                                                                                            <  2:752345/100‑1 (MQ=255)
                                                                                       gCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGa                                                                                                            >  1:752345/1‑100 (MQ=255)
                                                                                       gCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAgg                                                          >  1:238997/1‑150 (MQ=255)
                                                                                        cTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGc                                                         <  1:676137/150‑1 (MQ=255)
                                                                                         tCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGt                                                                                         <  1:144682/117‑1 (MQ=255)
                                                                                         tCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGt                                                                                         >  2:144682/1‑117 (MQ=255)
                                                                                         tCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCg                                                        >  1:330713/1‑150 (MQ=255)
                                                                                          cacaCGACCATGAACAGCCCAGCAACGAAGATGCTGGCAGCCGAGGAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCgg                                                       <  1:661950/150‑1 (MQ=255)
                                                                                           acacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                           >  2:9127/1‑81 (MQ=255)
                                                                                           acacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                           <  1:9127/81‑1 (MQ=255)
                                                                                           acacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAgg                                                                                              >  1:966067/1‑110 (MQ=255)
                                                                                           acacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAgg                                                                                              <  2:966067/110‑1 (MQ=255)
                                                                                             acGACCATGAACAGCACAGCAACGCAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaag                                                    >  1:1272969/1‑150 (MQ=255)
                                                                                             acGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaag                                                    >  2:845648/1‑150 (MQ=255)
                                                                                             acGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaag                                                    >  2:1381426/1‑150 (MQ=255)
                                                                                               gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaagaa                                                  <  1:845648/150‑1 (MQ=255)
                                                                                                 ccATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGAGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAaca                                                >  2:634340/1‑150 (MQ=255)
                                                                                                  cATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                               <  1:1019945/102‑1 (MQ=255)
                                                                                                  cATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                               >  2:1019945/1‑102 (MQ=255)
                                                                                                   aTGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACg                                              >  2:1295348/1‑150 (MQ=255)
                                                                                                     gAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAc                                                                                 <  1:79851/113‑1 (MQ=255)
                                                                                                     gAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAc                                                                                 >  2:79851/1‑113 (MQ=255)
                                                                                                      aaCAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                              >  1:1116773/1‑83 (MQ=255)
                                                                                                      aaCAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                              <  2:1116773/83‑1 (MQ=255)
                                                                                                            acagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                     <  1:948174/150‑1 (MQ=255)
                                                                                                                    gAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAgcgc                             >  2:175068/1‑150 (MQ=255)
                                                                                                                         tGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCga                        <  1:686666/150‑1 (MQ=255)
                                                                                                                            tGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCgatga                     >  2:1032051/1‑150 (MQ=255)
                                                                                                                            tGGCAGCCGACAAGGTAATGGTCCGCGAGGACCAGGTGTGGAACGATTGTTGAAAGATGGCGAGGTCAAGTTTGAGGGCAGTAAGGTGACGCGTGCACGACTGCTTGTAGAGGCGGAAGAATACGTAGATCTTCATCAGCGCGCCgatga                     >  2:90473/1‑150 (MQ=255)
                                                                                                                              gCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATc                   <  2:1272969/150‑1 (MQ=255)
                                                                                                                                  ccGAGAAGGTCATGGTACGCGAGGACCAGGTGTTGAACCAGTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCCTGTCGAGGCGGAAGTACACGTAGATCTTCATCAGCGCGCCGATGATCTGGt               >  2:87730/1‑150 (MQ=255)
                                                                                                                                    gagaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTg             >  2:1172238/1‑150 (MQ=255)
                                                                                                                                        aGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAAt         <  2:238997/150‑1 (MQ=255)
                                                                                                                                         acacctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATa        <  1:537953/145‑1 (MQ=255)
                                                                                                                                          cacctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATag       >  2:1163819/5‑150 (MQ=255)
                                                                                                                                           acctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATaga      >  2:371468/4‑150 (MQ=255)
                                                                                                                                            cATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGAt                                         >  2:1144164/1‑114 (MQ=255)
                                                                                                                                            cATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGAt                                         <  1:1144164/114‑1 (MQ=255)
                                                                                                                                              tGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGaa   <  1:1415537/150‑1 (MQ=255)
                                                                                                                                               ggTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt  <  1:1295348/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAT  >  NC_002947/1471254‑1471546

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