Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 696,413 A→G 6.7% intergenic (+137/+17) 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,4130AG 6.7% 328.2 / 10.1 105intergenic (+137/+17)ubiF/glnX2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base A (41/56);  new base G (3/4);  total (44/61)
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.06e-01

Only 100 of 106 total aligned reads displayed.

TTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTAC  >  NC_000913/696265‑696553
                                                                                                                                                    |                                                                                                                                            
ttAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGaa                                                                                                                                              <  3:139151/149‑1 (MQ=255)
      ttGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGc                                                                                                                                        >  7:214095/1‑149 (MQ=255)
         tAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCtt                                                                                                                                      <  7:269534/148‑1 (MQ=255)
            ccTTACAAGATTGCCGGATGCGTGCCAACCGTAGGTCTGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCtttttt                                                                                                                                  <  2:394106/149‑1 (MQ=255)
            ccTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCTTCGGCGAGCttttt                                                                                                                                   <  6:123941/148‑1 (MQ=255)
             cTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTa                                                                                                                                 <  8:238858/149‑1 (MQ=255)
                acaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAtg                                                                                                                              >  2:5925/1‑149 (MQ=255)
                acaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAt                                                                                                                               <  3:204882/148‑1 (MQ=255)
                 caacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAtgt                                                                                                                             <  3:306580/149‑1 (MQ=255)
                 caacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGAGAGCTTTTTTAAtgt                                                                                                                             >  3:9524/1‑149 (MQ=255)
                  aacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTCAATGACCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAtgtg                                                                                                                            >  5:66249/1‑149 (MQ=255)
                      ttGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTg                                                                                                                        <  8:201309/149‑1 (MQ=255)
                         ccGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTggg                                                                                                                      <  8:316010/148‑1 (MQ=255)
                          cGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTgggg                                                                                                                     >  2:389506/1‑148 (MQ=255)
                          cGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGt                                                                                                                    <  7:155563/149‑1 (MQ=255)
                           ggATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTa                                                                                                                   >  1:391598/1‑149 (MQ=255)
                              tGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGa                                                                                                                <  2:391598/149‑1 (MQ=255)
                                 gTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGa                                                                                                             >  6:116489/1‑149 (MQ=255)
                                   gCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGAtt                                                                                                           >  7:152895/1‑149 (MQ=255)
                                    ccAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTc                                                                                                          >  3:190963/1‑149 (MQ=255)
                                      aaCCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCg                                                                                                         >  5:377554/1‑148 (MQ=255)
                                      aaCCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCg                                                                                                         >  5:255009/1‑148 (MQ=255)
                                       aCCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCAATCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGaa                                                                                                       <  8:178170/149‑1 (MQ=255)
                                        ccGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACAGTGAAGGATAAGACGTGTCAACAGCGCATTCGACATTGAATGCACGCAGAAAAGCAAAAAGCTCGCCGCAGCGAGCGTGTTTAAGGTGGCTGGGGTACGAGGAGTCGacc                                                                                                      >  7:55372/1‑147 (MQ=38)
                                             ggTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCg                                                                                                  >  4:292224/1‑148 (MQ=255)
                                               tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGaa                                                                                               >  7:108686/1‑149 (MQ=255)
                                                cGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                                                                                              <  1:323752/149‑1 (MQ=255)
                                                  gATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGc                                                                                            >  4:210113/1‑149 (MQ=255)
                                                  gATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGc                                                                                            >  2:47114/1‑149 (MQ=255)
                                                  gATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGc                                                                                            >  2:320413/1‑149 (MQ=255)
                                                  gATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGc                                                                                            <  1:253578/149‑1 (MQ=255)
                                                       gACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGaa                                                                                       >  5:75057/1‑149 (MQ=255)
                                                       gACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGaa                                                                                       >  4:130490/1‑149 (MQ=255)
                                                         cgcgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGGCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATc                                                                                     <  2:365878/149‑1 (MQ=255)
                                                          gcgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGTAATGCCGGAATCa                                                                                    <  2:252276/149‑1 (MQ=255)
                                                          gcgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCa                                                                                    <  2:195943/149‑1 (MQ=255)
                                                          gcgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCa                                                                                    <  1:5925/149‑1 (MQ=255)
                                                          gcgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCa                                                                                    <  8:214095/149‑1 (MQ=255)
                                                           cgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAg                                                                                   <  3:130490/149‑1 (MQ=255)
                                                                  gTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                            <  3:292224/149‑1 (MQ=255)
                                                                   tCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                                           <  8:359291/149‑1 (MQ=255)
                                                                   tCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                                           <  8:108686/149‑1 (MQ=255)
                                                                   tCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                                           <  5:321301/149‑1 (MQ=255)
                                                                       aTCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCACGGAATCAGCGAATCAGAATCCGGTGCCt                                                                       >  7:185402/1‑149 (MQ=39)
                                                                        tCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCt                                                                       >  5:380546/1‑148 (MQ=255)
                                                                           taCATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcg                                                                   >  1:351785/2‑148 (MQ=255)
                                                                               ttGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                            >  6:334662/1‑136 (MQ=255)
                                                                               ttGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                            <  5:334662/136‑1 (MQ=255)
                                                                                  aaGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTgg                                                             <  6:82494/148‑1 (MQ=255)
                                                                                  aaGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                            <  1:47114/149‑1 (MQ=255)
                                                                                  aaGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                            >  7:14553/1‑149 (MQ=255)
                                                                                    ggATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGa                                                          <  6:377554/149‑1 (MQ=255)
                                                                                    ggATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGa                                                          >  3:96274/1‑149 (MQ=255)
                                                                                         aGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATAccc                                                     <  1:320413/149‑1 (MQ=255)
                                                                                         aGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATAcc                                                      <  1:389506/148‑1 (MQ=255)
                                                                                          gACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATAccc                                                     >  1:337939/1‑148 (MQ=255)
                                                                                             gtgtCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAc                                                 >  1:240679/1‑149 (MQ=255)
                                                                                              tgtCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATGAGAAGCCGGTGCCTTACCGGTTGGc                                                            >  4:344129/1‑137 (MQ=255)
                                                                                              tgtCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                            <  3:344129/137‑1 (MQ=255)
                                                                                               gtCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTg                                               >  4:171806/1‑149 (MQ=255)
                                                                                                tCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTgg                                              <  5:116489/149‑1 (MQ=255)
                                                                                                  aaCATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTg                                           >  7:258819/1‑149 (MQ=255)
                                                                                                  aaCATCGCATGCGACAGTGCAGGCACGCAGTACAGCAAAAATCGCGCCGACGCGAGCTTTTTTAATGTGGCTGTGGTACGAGCATTCGAACCGCGGAATGCCGGAATCAGAATCCGTTTCCTTACCGCTTGGCGAGACCCCAACTGGGt                                            <  2:319205/149‑1 (MQ=18)
                                                                                                      tCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCAc                                        <  6:75057/149‑1 (MQ=255)
                                                                                                        gCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACt                                       <  2:351785/148‑1 (MQ=255)
                                                                                                          aTTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTAc                                    <  2:240679/149‑1 (MQ=255)
                                                                                                           ttCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACa                                   <  4:96274/149‑1 (MQ=255)
                                                                                                              gACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAg                                 <  4:190963/148‑1 (MQ=255)
                                                                                                                cATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTa                              <  3:171806/149‑1 (MQ=255)
                                                                                                                 aTTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGt                               >  4:311359/1‑147 (MQ=255)
                                                                                                                    gaatgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGc                           <  8:14553/148‑1 (MQ=255)
                                                                                                                     aatgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGt                         >  3:191554/1‑149 (MQ=255)
                                                                                                                      atgaaCGCAGAAAATCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCAGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTc                        <  6:380546/149‑1 (MQ=255)
                                                                                                                      atgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTa                              >  5:299175/1‑143 (MQ=255)
                                                                                                                      atgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTa                              <  6:299175/143‑1 (MQ=255)
                                                                                                                         aaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCtt                      <  8:152895/148‑1 (MQ=255)
                                                                                                                           cGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGaa                   <  7:91955/149‑1 (MQ=255)
                                                                                                                              aGAAAAGCATAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                            <  4:205490/105‑1 (MQ=255)
                                                                                                                              aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAgg                                <  8:158343/133‑1 (MQ=255)
                                                                                                                              aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAgg                                >  7:158343/1‑133 (MQ=255)
                                                                                                                                     cAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATACGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTa                                                                     <  6:34415/89‑1 (MQ=255)
                                                                                                                                     cAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATACGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTa                                                                     >  5:34415/1‑89 (MQ=255)
                                                                                                                                       aaaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAg                                 >  4:380966/1‑123 (MQ=255)
                                                                                                                                       aaaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAg                                 <  3:380966/123‑1 (MQ=255)
                                                                                                                                       aaaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                            <  5:328543/80‑1 (MQ=255)
                                                                                                                                       aaaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                            >  6:328543/1‑80 (MQ=255)
                                                                                                                                       aaaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATa                                                        <  4:243627/100‑1 (MQ=255)
                                                                                                                                       aaaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATa                                                        >  3:243627/1‑100 (MQ=255)
                                                                                                                                        aaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACa                                   <  3:110668/120‑1 (MQ=18)
                                                                                                                                        aaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACa                                   <  3:116002/120‑1 (MQ=18)
                                                                                                                                        aaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACa                                   >  4:116002/1‑120 (MQ=18)
                                                                                                                                        aaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACa                                   >  4:110668/1‑120 (MQ=18)
                                                                                                                                        aaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGt                                            >  7:225175/1‑111 (MQ=255)
                                                                                                                                        aaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGt                                            <  8:225175/111‑1 (MQ=255)
                                                                                                                                        aaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTggg      <  2:337939/149‑1 (MQ=255)
                                                                                                                                        aaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTggg      <  8:25140/149‑1 (MQ=255)
                                                                                                                                          aGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAATGTAAGCGTCTTGAATAAATTGGCTGGGGt    <  4:9524/149‑1 (MQ=255)
                                                                                                                                            cTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGTGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGTCTGGGGTAc  <  4:191554/149‑1 (MQ=255)
                                                                                                                                              cGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGc          >  3:336986/1‑139 (MQ=255)
                                                                                                                                               gCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCg                                                                                          <  6:92495/58‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTAC  >  NC_000913/696265‑696553

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