Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,553 T→C 7.8% L171P (CTG→CCG) ‡ cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,5530TC 7.8% 411.7 / 18.2 128L171P (CTG→CCG) ‡cydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base T (61/57);  new base C (7/3);  total (68/60)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 137 total aligned reads displayed.

GGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGC  >  NC_000913/773405‑773698
                                                                                                                                                    |                                                                                                                                                 
 gCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgct                                                                                                                                                   >  2:716/1‑148 (MQ=255)
          ttcggaGCTTCGTTCCGCAGGGGGAAAGTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                                                                                                                                         >  4:561800/6‑149 (MQ=255)
              tAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggtg                                                                                                                                      >  2:635975/1‑148 (MQ=255)
                gCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGc                                                                                                                                   >  3:779327/1‑149 (MQ=255)
                  ttCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGCCGTGGTGAGCGt                                                                                                                                 >  3:110496/1‑149 (MQ=255)
                   tCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                                                                                                                                 <  1:617452/148‑1 (MQ=255)
                    cGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTa                                                                                                                                <  2:370013/148‑1 (MQ=255)
                     gTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAgg                                                                                                                              <  4:393969/149‑1 (MQ=255)
                      ttCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAggg                                                                                                                             <  3:100339/149‑1 (MQ=255)
                        ccgccgCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                                                                                                                           <  4:582246/149‑1 (MQ=255)
                          gccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatg                                                                                                                          >  1:533249/1‑148 (MQ=255)
                           ccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatga                                                                                                                         <  2:13773/148‑1 (MQ=255)
                           ccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatga                                                                                                                         >  1:711187/1‑148 (MQ=255)
                               tGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCAtt                                                                                                                    <  3:414450/149‑1 (MQ=255)
                               tGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCAt                                                                                                                     >  1:567088/1‑148 (MQ=255)
                                   aaTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTc                                                                                                                <  3:724957/149‑1 (MQ=255)
                                    aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTc                                                                                                                <  1:580764/148‑1 (MQ=255)
                                        gtgtAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGgc                                                                                                           <  3:879665/149‑1 (MQ=255)
                                            aGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCa                                                                                                        <  2:854070/148‑1 (MQ=255)
                                              cGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAAc                                                                                                      >  1:50710/1‑148 (MQ=255)
                                                ttCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTa                                                                                                   >  4:397638/1‑149 (MQ=255)
                                                  cGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTAt                                                                                                  <  1:861983/148‑1 (MQ=255)
                                                  cGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATc                                                                                                 <  4:864376/149‑1 (MQ=255)
                                                      gaCCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGAGCATTACTCAGGGCGCAACCTACCAGCa                                                                                             >  4:583382/2‑149 (MQ=255)
                                                        ccTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCATCCTATCTGCaa                                                                                            >  1:38004/1‑148 (MQ=255)
                                                        ccTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCaaa                                                                                           <  3:691153/149‑1 (MQ=255)
                                                         cTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAAt                                                                                          >  4:220853/1‑149 (MQ=255)
                                                         cTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCCAGCAGGccgaa                                                                                                                                        <  4:117923/103‑5 (MQ=255)
                                                           gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATg                                                                                         <  2:567088/148‑1 (MQ=255)
                                                           gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATg                                                                                         >  1:136654/1‑148 (MQ=255)
                                                           gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATg                                                                                         >  1:565528/1‑148 (MQ=255)
                                                                    cGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTgg                                                                               <  4:288108/149‑1 (MQ=255)
                                                                       aCCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGc                                                                             >  1:114668/1‑148 (MQ=255)
                                                                        ccGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCg                                                                            >  1:636444/1‑148 (MQ=255)
                                                                          gTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGaa                                                                          <  2:98767/148‑1 (MQ=255)
                                                                          gTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAAc                                                                         <  3:33939/149‑1 (MQ=255)
                                                                           ttCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACt                                                                        >  3:420682/1‑149 (MQ=255)
                                                                           ttCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACt                                                                        >  4:702818/1‑149 (MQ=255)
                                                                            tCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACt                                                                        >  1:728871/1‑148 (MQ=255)
                                                                               aCGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCa                                                                     >  2:740020/1‑148 (MQ=255)
                                                                                cGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCAc                                                                    <  1:831734/148‑1 (MQ=255)
                                                                                cGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCAc                                                                    <  2:636444/148‑1 (MQ=255)
                                                                                    gatgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCg                                                               >  3:299836/1‑149 (MQ=255)
                                                                                     atgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGt                                                              >  3:537356/1‑149 (MQ=255)
                                                                                     atgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGt                                                              >  4:846720/1‑149 (MQ=255)
                                                                                       gaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTa                                                             >  2:692241/1‑148 (MQ=255)
                                                                                            tCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGt                                                        <  1:716/148‑1 (MQ=255)
                                                                                            tCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTg                                                       >  3:213878/1‑149 (MQ=255)
                                                                                            tCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTg                                                       >  3:366887/1‑149 (MQ=255)
                                                                                               gCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCa                                                     >  2:101883/1‑148 (MQ=255)
                                                                                                  tCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACg                                                  <  1:871433/148‑1 (MQ=255)
                                                                                                   cTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGt                                                                                 >  2:888370/1‑116 (MQ=255)
                                                                                                   cTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGt                                                                                 <  1:888370/116‑1 (MQ=255)
                                                                                                    tGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGc                                                <  2:50710/148‑1 (MQ=255)
                                                                                                    tGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCt                                               >  3:422368/1‑149 (MQ=255)
                                                                                                         tacACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGt                                          <  4:761270/149‑1 (MQ=255)
                                                                                                         tacACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGt                                          <  3:702818/149‑1 (MQ=255)
                                                                                                         tacACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTAGCTGCAAATGCGTACCGTGGGCGAACGGCACCTGCGTACCCGTGCAACGGCTCAGGt                                          <  4:421590/149‑1 (MQ=255)
                                                                                                                    cttcttCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAgg                                           >  1:389553/1‑137 (MQ=255)
                                                                                                                    cttcttCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAgg                                           <  2:389553/137‑1 (MQ=255)
                                                                                                                    cttcttCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTg                                <  1:101883/148‑1 (MQ=255)
                                                                                                                     ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTgg                               >  1:316115/1‑148 (MQ=255)
                                                                                                                     ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTgg                               >  2:658369/1‑148 (MQ=255)
                                                                                                                     ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                              >  3:144782/1‑149 (MQ=255)
                                                                                                                     ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                              <  4:276537/149‑1 (MQ=255)
                                                                                                                      tcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                              >  1:302207/1‑148 (MQ=255)
                                                                                                                      tcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                              >  2:697763/1‑148 (MQ=255)
                                                                                                                       cttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTg                             <  2:741802/148‑1 (MQ=255)
                                                                                                                       cttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTg                             <  1:830103/148‑1 (MQ=255)
                                                                                                                       cttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGa                            >  3:73343/1‑149 (MQ=255)
                                                                                                                         tcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGac                           <  2:302207/148‑1 (MQ=255)
                                                                                                                             gTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                      >  4:894938/1‑149 (MQ=255)
                                                                                                                             gTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                      >  4:318692/1‑149 (MQ=255)
                                                                                                                             gTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTg                       >  1:620119/1‑148 (MQ=255)
                                                                                                                              ttGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                      >  2:357588/1‑148 (MQ=255)
                                                                                                                               tGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                      <  2:603160/147‑1 (MQ=255)
                                                                                                                               tGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                      >  1:603160/1‑147 (MQ=255)
                                                                                                                               tGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTc                    <  4:422368/149‑1 (MQ=255)
                                                                                                                                gCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCt                   >  4:185846/1‑149 (MQ=255)
                                                                                                                                  ttAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGt                 >  3:772490/1‑149 (MQ=255)
                                                                                                                                  ttAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGt                 >  3:669361/1‑149 (MQ=255)
                                                                                                                                   tAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGtt                <  4:420682/149‑1 (MQ=255)
                                                                                                                                   tAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGt                 >  2:271687/1‑148 (MQ=255)
                                                                                                                                     aCCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGttt               <  2:889277/148‑1 (MQ=255)
                                                                                                                                      cccGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTc              >  1:422488/1‑148 (MQ=255)
                                                                                                                                      cccGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTc              >  2:401523/1‑148 (MQ=255)
                                                                                                                                       ccGTTCGGCCTGCTGTCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGc            >  4:872017/1‑149 (MQ=255)
                                                                                                                                          ccacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAct         >  4:719744/5‑149 (MQ=255)
                                                                                                                                          ccacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAc          <  4:627487/144‑1 (MQ=255)
                                                                                                                                          ccacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAc          >  3:627487/5‑148 (MQ=255)
                                                                                                                                           cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGc                                                                <  3:661423/90‑1 (MQ=255)
                                                                                                                                           cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGc                                                                >  4:661423/4‑93 (MQ=255)
                                                                                                                                           cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActg        >  3:586139/4‑149 (MQ=255)
                                                                                                                                           cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActg        >  4:853753/4‑149 (MQ=255)
                                                                                                                                           cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAct         >  1:807346/4‑148 (MQ=255)
                                                                                                                                              gCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActggc      <  1:282262/148‑1 (MQ=255)
                                                                                                                                                 tgctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActggctgg   >  2:290185/1‑148 (MQ=255)
                                                                                                                                                 tgctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActggctgg   >  2:11460/1‑148 (MQ=255)
                                                                                                                                                 tgctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACtggctggc  <  3:699230/149‑1 (MQ=255)
                                                                                                                                                  gctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACtggctggc  <  1:635975/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGC  >  NC_000913/773405‑773698

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