Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,769,943 G→A 6.7% T290T (ACG→ACA)  ydiK → UPF0118 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,9430GA 6.7% 277.7 / 6.0 89T290T (ACG→ACA) ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (43/39);  new base A (4/2);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

ACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCG  >  NC_000913/1769797‑1770089
                                                                                                                                                  |                                                                                                                                                  
aCGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCt                                                                                                                                                  <  1:911670/149‑1 (MQ=255)
       tAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTTGTTGGCACGCTGGATAAc                                                                                                                                           >  2:1006243/1‑149 (MQ=255)
       tAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAAc                                                                                                                                           >  2:1047269/1‑149 (MQ=255)
        aaCGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACg                                                                                                                                          <  2:1196181/149‑1 (MQ=255)
         aCGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGt                                                                                                                                         >  1:611857/1‑149 (MQ=255)
           ggTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAAc                                                                                                                                           <  1:1466897/145‑1 (MQ=255)
           ggTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAAc                                                                                                                                           >  2:1466897/1‑145 (MQ=255)
            gtgtTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCAt                                                                                                                                      <  1:436982/149‑1 (MQ=255)
            gtgtTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCAt                                                                                                                                      <  2:920472/149‑1 (MQ=255)
              gtTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATcc                                                                                                                                    >  1:446681/1‑149 (MQ=255)
                 aaTGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGcc                                                                                                                                 >  2:1552128/1‑149 (MQ=255)
                 aaTGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGcc                                                                                                                                 >  2:1017610/1‑149 (MQ=255)
                        ctctCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGtt                                                                                                                          >  2:342366/1‑149 (MQ=255)
                        ctctCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACACTGGATAACGTCATCCACCCAATGtt                                                                                                                          >  2:345202/1‑149 (MQ=255)
                        ctctCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACACTGGATAACGTCATCCACCCAATGtt                                                                                                                          >  2:339470/1‑149 (MQ=255)
                        ctctCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATAACACCTGGGGAACGGTATTGGTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCGTCCGGCCAATGtt                                                                                                                          >  1:84412/1‑149 (MQ=255)
                         tctcCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTa                                                                                                                         >  1:1237291/1‑149 (MQ=255)
                          ctcCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTaa                                                                                                                        >  2:1545737/1‑149 (MQ=255)
                           accTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGGGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAAt                                                                                                                       <  1:1138889/148‑1 (MQ=255)
                               gCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGc                                                                                                                   >  2:67311/1‑149 (MQ=255)
                                ccTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCa                                                                                                                  >  1:1592485/1‑149 (MQ=255)
                                ccTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCa                                                                                                                  >  1:1593323/1‑149 (MQ=255)
                                ccTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCa                                                                                                                  >  1:1586407/1‑149 (MQ=255)
                                 ctcctCCACCTTGGCCCGTTGCCGGTACTGATTCCCGCGAATAGCTGGCGCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGGTAAGTCGCAt                                                                                                                 <  2:236463/145‑1 (MQ=255)
                                    ggccAGCTTGGCCCGTTGCCGGTACTTATTCCCGCGATGATCTGGCTCTACTTGTCTGTCGTTACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATggg                                                                                                              <  1:993302/147‑1 (MQ=255)
                                    gTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATggg                                                                                                              <  1:380200/149‑1 (MQ=255)
                                       cAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGc                                                                                                           >  2:947467/1‑149 (MQ=255)
                                        aGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGcc                                                                                                          <  1:707685/149‑1 (MQ=255)
                                         gCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCg                                                                                                         <  2:486715/149‑1 (MQ=255)
                                         gCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCg                                                                                                         >  1:595708/1‑149 (MQ=255)
                                           ttGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGAt                                                                                                       >  1:975138/1‑149 (MQ=255)
                                            tGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGAtt                                                                                                      <  1:1017610/149‑1 (MQ=255)
                                              gCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTa                                                                                                    >  1:626389/1‑149 (MQ=255)
                                              gCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTa                                                                                                    >  1:620764/1‑149 (MQ=255)
                                               cccGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTAc                                                                                                   >  2:1509595/1‑149 (MQ=255)
                                               cccGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTAc                                                                                                   <  2:1326734/149‑1 (MQ=255)
                                                ccGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGAGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTAcc                                                                                                  >  2:1049230/1‑149 (MQ=255)
                                                 cGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCg                                                                                                 <  1:1552128/149‑1 (MQ=255)
                                                      ccGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGAt                                                                                            >  1:592042/1‑149 (MQ=255)
                                                      ccGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGAGTGCTTGGCACTCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTTCCGCTgct                                                                                            >  2:290015/1‑147 (MQ=255)
                                                       cGGAACTGATACCGGCGATTACCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTCGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATc                                                                                           <  2:1089217/149‑1 (MQ=255)
                                                         gTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGACGATTTACCGCTGATCCt                                                                                         >  2:601460/1‑149 (MQ=255)
                                                         gTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTAGATCTCTGTGGTAGCCAAAAAATAAAAAAAAAAAAGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCt                                                                                         >  2:188500/1‑149 (MQ=255)
                                                           aCTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGa                                                                                       <  2:595708/149‑1 (MQ=255)
                                                                  ccGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATtctctc                                                                                <  2:1338057/149‑1 (MQ=255)
                                                                    ggCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAGCGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTg                                                                              >  2:1528408/1‑149 (MQ=255)
                                                                     gCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTgg                                                                             >  2:48080/1‑149 (MQ=255)
                                                                      cGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGc                                                                            <  2:626389/149‑1 (MQ=255)
                                                                      cGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGc                                                                            <  2:620764/149‑1 (MQ=255)
                                                                        aTTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTTGATAAAGTCATCCGCCCAAGGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGt                                                                          >  1:292656/1‑149 (MQ=255)
                                                                           aTCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCTCATGGGTGCCGATTTACCGCTGAGCCTGATTCTCTCTGGCGTTAt                                                                       >  2:792391/1‑149 (MQ=255)
                                                                            tCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTAtt                                                                      >  2:1293748/1‑149 (MQ=255)
                                                                             cTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTTGATCGCTGTGGTAGGCCAACAGGAAAAATTAATCAGCCCAATGTTAATTCGCATGTGTGCCGATTTACCGATGATCGTGATTCTCTATGGCGTTATtg                                                                     >  2:938824/1‑149 (MQ=255)
                                                                             cTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATtg                                                                     >  1:857136/1‑149 (MQ=255)
                                                                               ggCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTTGCACGCTGGATAACTTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATtggt                                                                   <  1:509904/149‑1 (MQ=255)
                                                                                  tctACTGGACTGGCGATACCAGGTGGGGAAAGGTTTTGTTTGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTggtggt                                                                <  2:1387386/149‑1 (MQ=255)
                                                                                  tctACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTggtggt                                                                >  2:1021549/1‑149 (MQ=255)
                                                                                  cctACTGGACTGGCGATACCACCTGGGGAACCGTATTGGTAGTGTGGATCGGGGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTggtggt                                                                <  2:874276/148‑1 (MQ=255)
                                                                                    tACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGttt                                                              >  1:626973/1‑149 (MQ=255)
                                                                                      ctggactggCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACACTGGATAACGTCATCCACCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGa                                                            <  1:339470/149‑1 (MQ=255)
                                                                                      ctggactggCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACACTGGATAACGTCATCCACCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGa                                                            <  1:345202/149‑1 (MQ=255)
                                                                                       tggactggCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGAt                                                           <  1:48080/149‑1 (MQ=255)
                                                                                            tggCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCtt                                                      >  1:1199694/1‑149 (MQ=255)
                                                                                            tggCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCtt                                                      <  1:1200821/149‑1 (MQ=255)
                                                                                                      accTGTGGAACGGTATTGTTAGAGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGgatgat                                            <  1:1006243/149‑1 (MQ=255)
                                                                                                        cTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCg                                          <  1:861937/149‑1 (MQ=255)
                                                                                                           gggAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCACTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTc                                       <  2:1237291/149‑1 (MQ=255)
                                                                                                              aaCGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGt                                    <  2:611857/149‑1 (MQ=255)
                                                                                                                cGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGttt                                  >  2:1260050/1‑149 (MQ=255)
                                                                                                                cGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTAGCATGGGTGCCGATTTACCGCTGATACTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTATGttt                                  >  1:1531020/1‑149 (MQ=255)
                                                                                                                 ggTAATGTTAGGGTGGAGCGGTGGGGTTGGCACGCTGGATAACGGCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTCCCGCTGAACCTGTATCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGGCTGTTTa                                 <  2:291230/149‑1 (MQ=255)
                                                                                                                  gTATTGTTAGTGTGGAGCTGTGGGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTTCCGGTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTAt                                <  2:1256948/149‑1 (MQ=255)
                                                                                                                   tATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTAtt                               >  1:15688/1‑149 (MQ=255)
                                                                                                                    aTTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTg                              <  1:1545737/149‑1 (MQ=255)
                                                                                                                     ttGTTAGTGTGGAGCGGTGGGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTgg                             <  1:67311/149‑1 (MQ=255)
                                                                                                                        ttAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTcc                          <  1:1047269/149‑1 (MQ=255)
                                                                                                                            tgtgGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGtt                      <  1:1509595/149‑1 (MQ=255)
                                                                                                                               ggAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTg                   <  2:446681/149‑1 (MQ=255)
                                                                                                                                     gtgtGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTTATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGcc             <  1:1021549/149‑1 (MQ=255)
                                                                                                                                     gtgtGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGcc             <  2:1593323/149‑1 (MQ=255)
                                                                                                                                     gtgtGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGcc             >  2:1214908/1‑149 (MQ=255)
                                                                                                                                          gTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCGGATTTACCGCTGCTCCTGATTCTCTCTGGCGTTTTTGGTGGTTTGATTGCTTTCGGGCTGCTCGGTCTGTGTATTGGTCGGGTTCTGTTAGCCGTTTc        >  1:301128/1‑149 (MQ=255)
                                                                                                                                          gTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTc        >  1:115968/1‑149 (MQ=255)
                                                                                                                                          gTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTc        >  2:406526/1‑149 (MQ=255)
                                                                                                                                             ggCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTg     <  1:1262992/149‑1 (MQ=255)
                                                                                                                                               cACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGc   <  1:947467/149‑1 (MQ=255)
                                                                                                                                                ccGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTGATTGCTGGTTTGATTGCTTTCGGGATGATAGGTCTGTTTATTGGTACGGTTATGTTAGCCGTTTCATGGCg  >  2:82386/2‑149 (MQ=255)
                                                                                                                                                aCGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCg  >  1:1057193/1‑149 (MQ=255)
                                                                                                                                                aCGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCg  >  2:968029/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCG  >  NC_000913/1769797‑1770089

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