Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 286,949 Δ1 bp coding (161/1611 nt) yagH → CP4‑6 prophage; putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913286,9460C.98.7% 317.6 / ‑2.9 76coding (158/1611 nt)yagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (47/28);  minor base G (0/1);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACTCGC  >  NC_000913/286808‑287077
                                                                                                                                          |                                                                                                                                   
tACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCAcc                                                                                                                                    >  1:420618/1‑140 (MQ=255)
tACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCAcc                                                                                                                                    >  2:453214/1‑140 (MQ=255)
     aCCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCAcccgttg                                                                                                                               >  2:465012/1‑136 (MQ=255)
     aCCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCAcccgttg                                                                                                                               >  1:232434/1‑136 (MQ=255)
      ccGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTgg                                                                                                                             >  2:649401/1‑140 (MQ=255)
      ccGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTgg                                                                                                                             >  2:365913/1‑140 (MQ=255)
          cTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACcg                                                                                                                         >  1:462866/1‑140 (MQ=255)
          cTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACcg                                                                                                                         >  1:319353/1‑140 (MQ=255)
          cTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACcg                                                                                                                         >  2:1169964/1‑140 (MQ=255)
          cTTCAACCCGGACCCGTCCCTGCGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACcg                                                                                                                         >  2:247167/1‑140 (MQ=255)
               aCCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCgtgt                                                                                                                    <  2:232434/140‑1 (MQ=255)
                   ggACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGAt                                                                                                                >  1:769246/1‑140 (MQ=255)
                    gACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCgtgt                                                                                                                    <  1:581344/135‑1 (MQ=255)
                    gACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCgtgt                                                                                                                    >  2:581344/1‑135 (MQ=255)
                                tgCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTg                                                                                                            <  2:1182854/131‑1 (MQ=255)
                                tgCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTg                                                                                                            >  1:1182854/1‑131 (MQ=255)
                                     ccAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCaa                                                                                              >  1:1254124/1‑140 (MQ=255)
                                     ccAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCaa                                                                                              >  2:1133149/1‑140 (MQ=255)
                                             aGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACt                                                                                      >  1:969286/1‑140 (MQ=255)
                                               gACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTcc                                                                                    <  2:211740/140‑1 (MQ=255)
                                               gACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTcc                                                                                    <  2:909252/140‑1 (MQ=255)
                                                    ctacATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCcggcgg                                                                               <  1:649401/140‑1 (MQ=255)
                                                    ctacATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCcggcgg                                                                               <  2:1082378/140‑1 (MQ=255)
                                                      acATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCa                                                                             <  2:420618/140‑1 (MQ=255)
                                                       cATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCcgg                                                                                  <  1:964980/134‑1 (MQ=255)
                                                       cATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCAt                                                                            >  1:1369631/1‑140 (MQ=255)
                                                       cATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCAt                                                                            >  1:885254/1‑140 (MQ=255)
                                                       cATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATACTGGACATGAAGGGCAACCCGGACTCcgg                                                                                  >  2:964980/1‑134 (MQ=255)
                                                        aTCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATc                                                                           >  2:999617/1‑140 (MQ=255)
                                                          cGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTg                                                                         >  2:88114/1‑140 (MQ=255)
                                                          cGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTg                                                                         >  2:1304002/1‑140 (MQ=255)
                                                          cGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTg                                                                         >  2:1226510/1‑140 (MQ=255)
                                                          cGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTg                                                                         >  2:482685/1‑140 (MQ=255)
                                                            ccACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTggg                                                                       <  1:573442/140‑1 (MQ=255)
                                                            ccACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTggg                                                                       >  2:9885/1‑140 (MQ=255)
                                                            ccACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCATTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTggg                                                                       >  2:422493/1‑140 (MQ=255)
                                                              aCCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgcg                                                                     >  1:147121/1‑140 (MQ=255)
                                                              aCCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgcg                                                                     >  1:811391/1‑140 (MQ=255)
                                                              aCCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgcg                                                                     <  1:453214/140‑1 (MQ=255)
                                                               ccTCGACCTTCGAGTGGTTTCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgcgc                                                                    <  2:1389449/140‑1 (MQ=255)
                                                               ccTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgcgc                                                                    >  1:432045/1‑140 (MQ=255)
                                                               ccTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgcgc                                                                    >  2:658054/1‑140 (MQ=255)
                                                                 tCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCg                                                                  >  2:1265330/1‑140 (MQ=255)
                                                                    aCCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGc                                                               <  1:365913/140‑1 (MQ=255)
                                                                    aCCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGc                                                               >  1:764900/1‑140 (MQ=255)
                                                                    aCCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGc                                                               <  2:1369631/140‑1 (MQ=255)
                                                                         cGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAg                                                          <  1:247167/140‑1 (MQ=255)
                                                                          gAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGc                                                               >  1:459075/1‑134 (MQ=255)
                                                                          gAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGc                                                               <  2:459075/134‑1 (MQ=255)
                                                                               gTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGc                                                    >  1:502486/1‑140 (MQ=255)
                                                                               gTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGCGCCTGAGCTACGc                                                    >  2:764865/1‑140 (MQ=255)
                                                                                      ggCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGt                                             <  1:1289789/140‑1 (MQ=255)
                                                                                          tGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAAt                                         >  1:1090823/1‑140 (MQ=255)
                                                                                          tGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAAt                                         >  2:1302819/1‑140 (MQ=255)
                                                                                          tGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAAt                                         >  1:514719/1‑140 (MQ=255)
                                                                                               aTCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTgg                                    <  1:1304002/140‑1 (MQ=255)
                                                                                               aTCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTgg                                    <  2:514719/140‑1 (MQ=255)
                                                                                               aTCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTgg                                    <  1:482685/140‑1 (MQ=255)
                                                                                                tCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGgc                                   <  1:1169964/140‑1 (MQ=255)
                                                                                                 cTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGgct                                  >  2:832519/1‑140 (MQ=255)
                                                                                                 cTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGgct                                  >  2:460683/1‑140 (MQ=255)
                                                                                                  tACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGgctg                                 <  1:1133149/140‑1 (MQ=255)
                                                                                                    ccACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGgctgct                               <  1:676124/140‑1 (MQ=255)
                                                                                                           cGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGgctgct                               >  1:1342404/1‑133 (MQ=255)
                                                                                                           cGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGgctgct                               <  2:1342404/133‑1 (MQ=255)
                                                                                                           cGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACAcc                        >  1:1238639/1‑140 (MQ=255)
                                                                                                           cGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACAcc                        <  2:462866/140‑1 (MQ=255)
                                                                                                           cGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACAcc                        >  1:120744/1‑140 (MQ=255)
                                                                                                           cGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCCGGCTGCTCTACAcc                        <  1:1265330/140‑1 (MQ=255)
                                                                                                            gTGACCTGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCg                       >  2:749168/1‑140 (MQ=255)
                                                                                                                  tGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGa                 >  1:578593/1‑140 (MQ=255)
                                                                                                                  tGAAAAACTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGa                 <  1:1226510/140‑1 (MQ=255)
                                                                                                                    aaaaaCTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAg               >  2:108263/1‑140 (MQ=255)
                                                                                                                    aaaaaCTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAg               >  2:401397/1‑140 (MQ=255)
                                                                                                                    aaaaaCTGGTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGCATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGa              >  1:612561/1‑140 (MQ=255)
                                                                                                                       aaCTGGTCGCTGGTCAGCGCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGAtt            <  1:764865/140‑1 (MQ=255)
                                                                                                                            gTCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGa       >  1:1360846/1‑140 (MQ=255)
                                                                                                                             tCGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGCATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACt     <  2:612561/140‑1 (MQ=255)
                                                                                                                              cGCTGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACt     >  1:1287365/1‑140 (MQ=255)
                                                                                                                                 tGGTCAGCACCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACTCGc  >  2:486876/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                   
TACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCTGAGCTACGCCGACGGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACTCGC  >  NC_000913/286808‑287077

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