Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 696,328 A→C 12.5% intergenic (+52/+102) ubiF → / ← glnX 2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913696,3280AC12.5% 219.7 / 23.7 72intergenic (+52/+102)ubiF/glnX2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base A (28/35);  new base C (3/6);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

ACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAA  >  NC_000913/696183‑696427
                                                                                                                                                 |                                                                                                   
accacTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                                     <  4:8933/146‑1 (MQ=255)
accacTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                                     >  3:8933/1‑146 (MQ=255)
  cacTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTc                                                                                                <  1:645414/149‑1 (MQ=255)
   acTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCg                                                                                               <  2:149817/149‑1 (MQ=255)
   acTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCg                                                                                               <  2:1386025/149‑1 (MQ=255)
   acTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCg                                                                                               >  3:824827/1‑149 (MQ=255)
    cTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGTCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTctt                                                                                              <  1:146162/149‑3 (MQ=255)
    cTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGTCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTctt                                                                                              <  1:409895/149‑3 (MQ=255)
     tGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTctta                                                                                             <  3:709195/149‑4 (MQ=255)
       cGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                                           <  3:439017/149‑1 (MQ=255)
       cGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATc                                                                                           <  3:1749668/149‑1 (MQ=255)
                 gTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTc                                                                                                >  2:193120/1‑134 (MQ=255)
                 gTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTc                                                                                                <  1:193120/134‑1 (MQ=255)
                 gTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                 <  2:1473471/149‑1 (MQ=255)
                 gTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                 <  3:128941/149‑1 (MQ=255)
                 gTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                 <  3:1679029/149‑1 (MQ=255)
                   aaTCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAgg                                                                               <  4:824827/149‑1 (MQ=255)
                          ggTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                 <  1:1323083/140‑1 (MQ=255)
                          ggTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                 >  2:1323083/1‑140 (MQ=255)
                            ttAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACgt                                                                      >  1:306489/1‑149 (MQ=255)
                                tGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                                            >  2:1327250/1‑119 (MQ=255)
                                tGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                                            <  1:1327250/123‑5 (MQ=255)
                                           cGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                                            <  4:1155357/112‑5 (MQ=255)
                                           cGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATGCGACGCTGCCGCGTcttat                                                                                            >  3:1155357/1‑108 (MQ=255)
                                                    gtgtTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaa                                              <  2:306489/149‑1 (MQ=255)
                                                    gtgaTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaa                                              >  1:1618166/1‑149 (MQ=255)
                                                     tgtTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaat                                             >  4:1697189/1‑149 (MQ=255)
                                                       ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCg                                                      <  3:1237810/138‑1 (MQ=255)
                                                       ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCg                                                      >  2:977181/1‑138 (MQ=255)
                                                       ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCg                                                      <  1:977181/138‑1 (MQ=255)
                                                       ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTAGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCg                                                      >  4:1237810/1‑138 (MQ=255)
                                                        tGAAACGTGAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaatgaa                                          <  2:1618166/149‑1 (MQ=255)
                                                         gAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACAt                                                                                     <  1:1278372/105‑1 (MQ=255)
                                                         gAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACAt                                                                                     >  2:1278372/1‑105 (MQ=255)
                                                         gAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAAc                                         >  2:134058/1‑149 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCTGCCGCGTcttat                                                                                            >  4:869070/1‑93 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCTGCCGCGTcttat                                                                                            <  3:869070/97‑5 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGc                                                                                              <  2:1881300/95‑1 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGc                                                                                              >  1:1881300/1‑95 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATGAGACGTGTc                                                                   <  2:408529/122‑1 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTc                                                                   >  1:408529/1‑122 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTc                                                                   >  4:799247/1‑122 (MQ=255)
                                                          aaaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTc                                                                   <  3:799247/122‑1 (MQ=255)
                                                           aaCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGc                                       >  4:563564/1‑149 (MQ=255)
                                                                 aGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                                 >  4:843589/1‑85 (MQ=255)
                                                                 aGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGt                                                                                                 <  3:843589/85‑1 (MQ=255)
                                                                 aGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGAc                                                    <  2:1429473/130‑1 (MQ=255)
                                                                 aGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGAc                                                    >  1:1429473/1‑130 (MQ=255)
                                                                 aGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaaa                                 <  3:563564/149‑1 (MQ=255)
                                                                 aGGAGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAAACGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaaa                                 <  1:134058/149‑1 (MQ=255)
                                                                       tGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACa                                                                                      >  2:1547973/1‑90 (MQ=255)
                                                                       tGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACa                                                                                      <  1:1547973/90‑1 (MQ=255)
                                                                       tGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaa                                              <  2:612721/130‑1 (MQ=255)
                                                                       tGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATtgaa                                              >  1:612721/1‑130 (MQ=255)
                                                                            tatGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTc                      <  1:1165681/149‑1 (MQ=255)
                                                                            tatGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTc                      >  2:1165681/1‑149 (MQ=255)
                                                                            tatGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTc                      <  3:1697189/149‑1 (MQ=255)
                                                                             atGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGCACGCAGAAAAGCAAAAAGCTc                      >  4:286593/1‑148 (MQ=255)
                                                                             atGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTc                      <  3:286593/148‑1 (MQ=255)
                                                                                 gTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCaaa                            >  3:756657/1‑138 (MQ=255)
                                                                                 gTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCaaa                            <  4:756657/138‑1 (MQ=255)
                                                                                        ttGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaaa                                 <  3:1241470/126‑1 (MQ=255)
                                                                                        ttGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGaaaa                                 >  4:1241470/1‑126 (MQ=255)
                                                                                          gTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                                     <  2:495097/56‑1 (MQ=255)
                                                                                          gTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCa                                                                                                     >  1:495097/1‑56 (MQ=255)
                                                                                           tAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCa                              >  2:835911/1‑126 (MQ=255)
                                                                                           tAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCa                              <  1:835911/126‑1 (MQ=255)
                                                                                                ttACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTaa  >  1:1133770/1‑149 (MQ=255)
                                                                                                  acaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                 <  4:1579460/68‑1 (MQ=255)
                                                                                                  acaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGaa                                                                                 >  3:1579460/1‑68 (MQ=255)
                                                                                                     acaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCGGaaaa                                 <  4:1705920/113‑1 (MQ=255)
                                                                                                     acaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCGGaaaa                                 >  3:1705920/1‑113 (MQ=255)
                                                                                                                          ccGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAtt                                                 <  3:1440564/76‑1 (MQ=255)
                                                                                                                          ccGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACAtt                                                 >  4:1440564/1‑76 (MQ=255)
                                                                                                                                         gACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGaa                >  4:783116/1‑94 (MQ=255)
                                                                                                                                         gACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGaa                <  3:783116/94‑1 (MQ=255)
                                                                                                                                                 |                                                                                                   
ACCACTGCGTTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAA  >  NC_000913/696183‑696427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: