Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 431,111 A→T 15.8% intergenic (+159/+18) yajD → / ← tsx HNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913431,1110AT15.8% 245.8 / 35.6 95intergenic (+159/+18)yajD/tsxHNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base A (43/37);  new base T (4/11);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 8.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

ATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCT  >  NC_000913/431013‑431210
                                                                                                  |                                                                                                   
atTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAAt                                                                                                    >  1:44483/1‑100 (MQ=255)
   gAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATccc                                                                                                 >  2:1013659/1‑100 (MQ=255)
   gAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATccc                                                                                                 <  2:976463/100‑1 (MQ=255)
   gAATTTACGTGCTTCTATTTGCCGGACAAAGCGTTTACGCCGCATCCGGCATTAACAAAGCACACGTTGTTAACAAGCAGAAATGCCGGGAATAAATccc                                                                                                 <  1:2052821/100‑1 (MQ=37)
     aTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCgg                                                                                               <  2:2340463/100‑1 (MQ=255)
     aTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCgg                                                                                               <  2:307988/100‑1 (MQ=255)
       ttACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCa                                                                                             <  1:1304283/100‑1 (MQ=255)
        tACGTGCTTTTGTTGGCCGGACAAAGCGTTTACTCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAt                                                                                            <  1:1356435/100‑1 (MQ=37)
        tACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAt                                                                                            >  2:977484/1‑100 (MQ=255)
        tACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAt                                                                                            <  1:2488544/100‑1 (MQ=255)
        tACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAt                                                                                            >  1:1138197/1‑100 (MQ=255)
         aCGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtt                                                                                           <  1:1845866/100‑1 (MQ=37)
         aCGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtt                                                                                           <  1:2090157/100‑1 (MQ=37)
         aCGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAtt                                                                                           >  2:2252443/1‑100 (MQ=255)
          cGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAttt                                                                                          >  1:653793/1‑100 (MQ=255)
           gTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttc                                                                                         <  2:1544694/100‑2 (MQ=37)
            tGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAttttt                                                                                        <  2:1752381/100‑2 (MQ=37)
            tGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttct                                                                                        <  2:668870/100‑3 (MQ=37)
            tGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttct                                                                                        <  1:328919/100‑3 (MQ=25)
            tGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttct                                                                                        <  2:697339/100‑3 (MQ=37)
            tGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttct                                                                                        <  1:185216/100‑3 (MQ=37)
            tGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTc                                                                                        <  1:1999995/100‑1 (MQ=255)
             gCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttctg                                                                                       <  1:1066479/100‑4 (MQ=37)
             gCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCa                                                                                       <  1:2218182/100‑1 (MQ=255)
              cTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCAt                                                                                      >  2:1609751/1‑100 (MQ=255)
              cTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCAt                                                                                      <  2:1930826/100‑1 (MQ=255)
               ttttGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATa                                                                                     <  2:137625/100‑1 (MQ=255)
               ttttGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATa                                                                                     <  2:1888068/100‑1 (MQ=255)
               ttttGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATa                                                                                     >  1:697811/1‑100 (MQ=255)
                 ttgttgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttctgattg                                                                                   <  2:2193937/100‑8 (MQ=25)
                  tgttgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATc                                                                                  <  2:1727370/100‑1 (MQ=255)
                     tgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGa                                                                               >  2:127356/1‑100 (MQ=255)
                      ggCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGaa                                                                              >  1:2096660/1‑100 (MQ=255)
                       gCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAg                                                                             >  2:1714418/1‑100 (MQ=255)
                       gCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAg                                                                             >  2:2073838/1‑100 (MQ=255)
                       gCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAg                                                                             >  2:2433004/1‑100 (MQ=255)
                       gCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAg                                                                             >  2:1803253/1‑100 (MQ=255)
                       gCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAg                                                                             >  2:108055/1‑100 (MQ=255)
                        ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGt                                                                            <  2:795484/100‑1 (MQ=255)
                        ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGt                                                                            <  2:1387414/100‑1 (MQ=255)
                           gACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGt                                                                         >  2:813923/1‑100 (MQ=255)
                               aaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAAcc                                                                     <  1:1155870/100‑1 (MQ=255)
                                aGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATACATCCCGGCATTTTCATAATCAGAAGTTGTAACct                                                                    >  2:1069404/1‑100 (MQ=255)
                                      tACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTAcc                                                              >  2:140116/1‑100 (MQ=255)
                                       aCGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCa                                                             >  1:225940/1‑100 (MQ=255)
                                         gccgcATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAgg                                                           <  2:423665/100‑1 (MQ=255)
                                          ccgcATCCGGCATGAACAAAGCACACGTTGGTAGCAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCAACTACCAGGt                                                          >  1:1700263/1‑100 (MQ=255)
                                           cgcATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAGTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTa                                                         <  1:43822/100‑1 (MQ=255)
                                             cATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAac                                                       <  2:1825367/100‑1 (MQ=255)
                                              aTCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAacc                                                      <  2:469094/100‑1 (MQ=255)
                                              aTCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAacc                                                      >  2:2209531/1‑100 (MQ=255)
                                               tCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAacca                                                     >  2:1470757/1‑100 (MQ=255)
                                               tCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAacca                                                     <  2:1499301/100‑1 (MQ=255)
                                                ccGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccac                                                    >  1:1810092/1‑100 (MQ=255)
                                                  ggCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                  <  2:1064678/100‑1 (MQ=255)
                                                   gCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAcccc                                                 <  1:152404/100‑1 (MQ=255)
                                                    cATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCa                                                >  2:2159474/1‑100 (MQ=255)
                                                      tGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGc                                              >  2:2236369/1‑100 (MQ=255)
                                                       gAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGcc                                             <  2:776172/100‑1 (MQ=255)
                                                       gAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGcc                                             >  2:464394/1‑100 (MQ=255)
                                                       gAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGcc                                             <  2:589063/100‑1 (MQ=255)
                                                       gAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAAACAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGcc                                             >  2:1921370/1‑100 (MQ=255)
                                                          cAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGt                                          >  2:2252562/1‑100 (MQ=255)
                                                           aaaGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTa                                         <  1:2309141/100‑1 (MQ=255)
                                                             aGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTaga                                       >  2:1008072/1‑100 (MQ=255)
                                                             aGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTaga                                       >  2:1471361/1‑100 (MQ=255)
                                                               cacaCGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGc                                     <  1:1226004/100‑1 (MQ=255)
                                                                   cGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAAc                                 <  1:2258575/100‑1 (MQ=255)
                                                                    gttgttAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACg                                >  2:1350856/1‑100 (MQ=255)
                                                                      tgttAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGtt                              <  2:1602978/100‑1 (MQ=255)
                                                                       gttAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTg                             <  2:1110019/100‑1 (MQ=255)
                                                                        ttAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTGATAATCAGAAGTTGTAACCTACTACCAGGTAACTACCCCAGGCGGTAGAGCGAACgg                               <  1:1822715/97‑2 (MQ=255)
                                                                          aaCAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAg                          >  2:1100710/1‑100 (MQ=255)
                                                                            cAATCAGAAATGACGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgtt                        >  2:1704524/1‑100 (MQ=255)
                                                                             aaTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttg                       >  2:69565/1‑100 (MQ=255)
                                                                              aTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTCACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgc                      >  1:2154869/1‑100 (MQ=255)
                                                                                cAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccg                    <  1:2089898/100‑1 (MQ=255)
                                                                                cAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccg                    <  2:2524987/100‑1 (MQ=255)
                                                                                   aaaTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGACgttgccgttg                 >  2:2542102/1‑100 (MQ=38)
                                                                                   aaaTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttg                 <  1:648276/100‑1 (MQ=255)
                                                                                    aaTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgc                >  2:1278697/1‑100 (MQ=255)
                                                                                    aaTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgc                >  2:1651686/1‑100 (MQ=39)
                                                                                    aaTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgc                >  1:740521/1‑100 (MQ=255)
                                                                                    aaTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgc                <  2:2126011/100‑1 (MQ=255)
                                                                                     aTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgcc               >  2:6089/1‑100 (MQ=255)
                                                                                     aTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgcc               >  2:2235769/1‑100 (MQ=255)
                                                                                       gCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGa             >  2:805252/1‑100 (MQ=255)
                                                                                        ccGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGaa            >  1:343965/1‑100 (MQ=255)
                                                                                          gggAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagt          >  1:1129148/1‑100 (MQ=255)
                                                                                            gAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttc        <  2:1528784/100‑1 (MQ=255)
                                                                                             aaTAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttca       <  1:603112/100‑1 (MQ=255)
                                                                                              aTAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcag      <  2:653796/100‑1 (MQ=255)
                                                                                              aTAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcag      >  2:2336157/1‑100 (MQ=255)
                                                                                              aTAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCAGTAAAGCGAACGTTGAAGTTGCCGTTACCGAagttca       >  1:996782/1‑99 (MQ=255)
                                                                                               tAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagt     >  2:1568940/1‑100 (MQ=255)
                                                                                                 aaTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagttc   >  1:2438939/1‑100 (MQ=255)
                                                                                                  aTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCt  >  2:314303/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCT  >  NC_000913/431013‑431210

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