Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 843,749 C→A D578Y (GAT→TAT)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,7490CA96.7% 322.0 / 2.4 92D578Y (GAT→TAT) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (56/33);  total (58/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATG  >  NC_000913/843655‑843844
                                                                                              |                                                                                               
atgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAg                                                                                            >  2:493727/1‑100 (MQ=255)
   atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgt                                                                                         >  1:633296/1‑100 (MQ=255)
     gtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTc                                                                                       <  1:1727456/100‑1 (MQ=255)
     gtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTc                                                                                       >  2:1631481/1‑100 (MQ=255)
       gATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAt                                                                                     <  1:576627/100‑1 (MQ=255)
       gATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAt                                                                                     >  1:1784486/1‑100 (MQ=255)
       gATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAt                                                                                     >  2:1065649/1‑100 (MQ=255)
        aTACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATg                                                                                    <  1:1057585/100‑1 (MQ=255)
         tACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGc                                                                                   >  2:1392344/1‑100 (MQ=255)
          aCGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGcc                                                                                  >  1:1170914/1‑100 (MQ=255)
          aCGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGcc                                                                                  <  2:1385901/100‑1 (MQ=255)
            gCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCg                                                                                 >  1:1688845/1‑99 (MQ=255)
               ncGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtttt                                                                             >  1:2048565/2‑100 (MQ=255)
                 ggTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCa                                                                           >  2:1567806/1‑100 (MQ=255)
                 ggTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCa                                                                           <  1:1065646/100‑1 (MQ=255)
                    aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAAc                                                                        >  2:2286727/1‑100 (MQ=255)
                    aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAAc                                                                        <  1:763403/100‑1 (MQ=255)
                      ccTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTg                                                                      <  2:741257/100‑1 (MQ=255)
                       cTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGa                                                                     >  1:254707/1‑100 (MQ=255)
                        tGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaa                                                                    >  1:1453102/1‑100 (MQ=255)
                         gTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaat                                                                   <  1:1556639/100‑1 (MQ=255)
                           cGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataa                                                                 >  2:1010371/1‑100 (MQ=255)
                            gTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATCGTCGCCAAATATCCAGTGGTCATGCCGTTTTCAAACTGaataaa                                                                >  1:1550664/1‑100 (MQ=255)
                             tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                               >  1:844071/1‑100 (MQ=255)
                             tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                               <  2:65023/100‑1 (MQ=255)
                               cGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGCCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAc                                                             >  2:1091475/1‑100 (MQ=255)
                                 cccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAccc                                                           >  1:2334760/1‑100 (MQ=255)
                                   cACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCg                                                         >  1:1828226/1‑100 (MQ=255)
                                   cACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCg                                                         >  1:2359389/1‑100 (MQ=255)
                                      ggAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTg                                                      >  1:997035/1‑100 (MQ=255)
                                      ggAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTg                                                      <  1:1840927/100‑1 (MQ=255)
                                      ggAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTAAATAAATACCCCGGTg                                                      >  1:108870/1‑100 (MQ=255)
                                       gAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGa                                                     <  1:899256/100‑1 (MQ=255)
                                         gCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGata                                                   <  1:469173/100‑1 (MQ=255)
                                         gCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGata                                                   <  1:1291485/100‑1 (MQ=255)
                                           gAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataat                                                 >  2:884644/1‑100 (MQ=255)
                                            aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                                >  1:567754/1‑100 (MQ=255)
                                            aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                                >  2:2202584/1‑100 (MQ=255)
                                            aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                                >  1:1345989/1‑100 (MQ=255)
                                              tCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATc                                              >  2:441497/1‑100 (MQ=255)
                                               cGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAatatn                                              <  2:2048563/99‑2 (MQ=255)
                                               cGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTGTTCAAACTGAATAAATACCCCGGTGATAATATCt                                             >  2:347893/1‑100 (MQ=255)
                                                 aCATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCttt                                           <  1:1356156/100‑1 (MQ=255)
                                                    tCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAc                                        <  2:1449631/100‑1 (MQ=255)
                                                    tCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAc                                        <  2:1170916/100‑1 (MQ=255)
                                                     ccGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAcc                                       >  2:424554/1‑100 (MQ=255)
                                                     ccGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAcc                                       <  1:493727/100‑1 (MQ=255)
                                                      cGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCa                                      >  1:1706777/1‑100 (MQ=255)
                                                      cGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCa                                      >  1:714593/1‑100 (MQ=255)
                                                       gTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAg                                     >  2:638054/1‑100 (MQ=255)
                                                       gTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAg                                     <  1:193653/100‑1 (MQ=255)
                                                             cTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTg                               <  2:1784487/100‑1 (MQ=255)
                                                              tgtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGc                              <  2:254707/100‑1 (MQ=255)
                                                               gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                             <  1:778204/100‑1 (MQ=255)
                                                               gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                             >  2:610729/1‑100 (MQ=255)
                                                               gatctggTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                             <  2:1095874/95‑1 (MQ=255)
                                                                tgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGa                            >  2:1433072/1‑100 (MQ=255)
                                                                  ccGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAc                          <  1:1010369/100‑1 (MQ=255)
                                                                  ccGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAAGAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                           >  2:2014620/1‑99 (MQ=255)
                                                                   cGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAcc                         >  1:455101/1‑100 (MQ=255)
                                                                    ggTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGagatc                        >  1:1095871/1‑96 (MQ=255)
                                                                    ggTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCa                        >  2:2169652/1‑100 (MQ=255)
                                                                    ggTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCa                        >  1:1140987/1‑100 (MQ=255)
                                                                     gTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCaa                       <  2:633297/100‑1 (MQ=255)
                                                                       cAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAc                     <  2:1640411/100‑1 (MQ=255)
                                                                        aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                    >  1:1213014/1‑100 (MQ=255)
                                                                        aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                    <  1:1739554/100‑1 (MQ=255)
                                                                         aCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACga                   >  1:738662/1‑100 (MQ=255)
                                                                         aCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACga                   >  1:1390107/1‑100 (MQ=255)
                                                                         aCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACga                   >  1:1592815/1‑100 (MQ=255)
                                                                         aCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                    >  1:738699/1‑99 (MQ=255)
                                                                            gCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                >  1:1721782/1‑100 (MQ=255)
                                                                            gCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                >  1:1950049/1‑100 (MQ=255)
                                                                            gCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                >  1:803529/1‑100 (MQ=255)
                                                                              ccGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAg              >  2:2048479/1‑100 (MQ=255)
                                                                              ccGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAg              <  2:2359386/100‑1 (MQ=255)
                                                                               cGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGc             <  2:398154/100‑1 (MQ=255)
                                                                                   aGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGa         <  1:1437956/100‑1 (MQ=255)
                                                                                    gTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAc        >  2:50933/1‑100 (MQ=255)
                                                                                    gTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAc        >  1:2304231/1‑100 (MQ=255)
                                                                                    gTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAc        >  1:2324074/1‑100 (MQ=255)
                                                                                    gTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAc        <  2:1261366/100‑1 (MQ=255)
                                                                                     tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc       >  1:144056/1‑100 (MQ=255)
                                                                                     tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc       >  2:1604678/1‑100 (MQ=255)
                                                                                      cACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCt      <  1:82120/100‑1 (MQ=255)
                                                                                       aCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTa     >  2:2355312/1‑100 (MQ=255)
                                                                                        ccAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTaa    >  1:840223/1‑100 (MQ=255)
                                                                                        ccAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTaa    >  2:1798660/1‑100 (MQ=255)
                                                                                         cAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAt   <  2:2162316/100‑1 (MQ=255)
                                                                                         cAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAt   >  1:1876956/1‑100 (MQ=255)
                                                                                          aaaTCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATg  >  1:2392201/1‑100 (MQ=255)
                                                                                          aaaTATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATg  <  2:1592816/100‑1 (MQ=255)
                                                                                              |                                                                                               
ATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATG  >  NC_000913/843655‑843844

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