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,7490CA99.0% 359.6 / ‑3.6 103D578Y (GAT→TAT) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (51/51);  total (51/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 103 total aligned reads displayed.

GAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGC  >  NC_000913/843653‑843845
                                                                                                |                                                                                                
cgatcttgtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATcc                                                                                               <  1:841802/94‑1 (MQ=255)
  atgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAg                                                                                             <  2:2240379/100‑1 (MQ=255)
  atgatgGGCTACGCCCCGGTCTCCTGTCTTACACCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAACATCCAg                                                                                             <  2:1018359/100‑1 (MQ=255)
   tgatgTGATACGCCCCGGTATCCTGTCGTTCGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgt                                                                                            <  2:1384575/100‑1 (MQ=255)
   tgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgt                                                                                            >  2:2060327/1‑100 (MQ=255)
    gatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtg                                                                                           <  1:1010805/100‑1 (MQ=255)
    gatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtg                                                                                           <  1:1914353/100‑1 (MQ=255)
    gatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtg                                                                                           <  2:359612/100‑1 (MQ=255)
    gatgTGACACGACCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtg                                                                                           >  2:1209455/1‑100 (MQ=255)
     cgatcttACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgt                                                                                          <  1:601229/94‑1 (MQ=255)
     atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgt                                                                                          >  1:503267/1‑100 (MQ=255)
     atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgt                                                                                          >  2:1122737/1‑100 (MQ=255)
     atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgt                                                                                            <  2:1943982/98‑1 (MQ=255)
      tgtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgatcg                                                                                         >  2:841800/1‑96 (MQ=255)
      tgtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGtt                                                                                         >  2:1069486/1‑100 (MQ=255)
        tgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCa                                                                                       >  1:325809/1‑100 (MQ=255)
         gATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAt                                                                                      >  2:1098222/1‑100 (MQ=255)
         gATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAt                                                                                      >  2:2112963/1‑100 (MQ=255)
          aTACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTAAACGGCCCGATAGTCACCAAATATCCAGTGTTCa                                                                                       <  2:1869239/98‑1 (MQ=255)
           tACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgtagatcg                                                                                    >  2:601230/1‑94 (MQ=255)
           tACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGc                                                                                    <  1:2112971/100‑1 (MQ=255)
            aCGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGcc                                                                                   >  1:520462/1‑100 (MQ=255)
             caccccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCg                                                                                  >  1:1552738/3‑100 (MQ=255)
             cGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCg                                                                                  >  2:1894537/1‑100 (MQ=255)
               ccccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtt                                                                                >  2:1004223/1‑100 (MQ=255)
                  cGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTc                                                                             <  2:473785/100‑1 (MQ=255)
                   ggTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCa                                                                            >  2:565964/1‑100 (MQ=255)
                   ggTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCa                                                                            >  1:1781325/1‑100 (MQ=255)
                      aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAAc                                                                         <  1:1122739/100‑1 (MQ=255)
                       atcTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACt                                                                        <  1:204328/98‑1 (MQ=255)
                        ccTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTg                                                                       <  1:1354090/100‑1 (MQ=255)
                        ccTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTg                                                                       <  1:657041/100‑1 (MQ=255)
                          tGTCGTACGCCCACGGAGCGAATCTTCATCCCTTCCACTGTGCCGGTCAACGGCCCGTTAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAAc                                                                         <  2:505797/96‑1 (MQ=255)
                          tGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaa                                                                     <  1:565963/100‑1 (MQ=255)
                           gTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACtagat                                                                    >  2:204329/1‑96 (MQ=255)
                           gTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaat                                                                    <  1:358105/100‑1 (MQ=255)
                           gTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaat                                                                    <  1:8019/100‑1 (MQ=255)
                            tCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaata                                                                   <  1:916565/100‑1 (MQ=255)
                            tCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaata                                                                   <  1:1041720/100‑1 (MQ=255)
                            tCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaata                                                                   <  2:427470/100‑1 (MQ=255)
                             cGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataa                                                                  >  2:1887607/1‑100 (MQ=255)
                             cGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataa                                                                  >  2:348046/1‑100 (MQ=255)
                               tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                                <  1:1069488/100‑1 (MQ=255)
                               tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                                >  1:58543/1‑100 (MQ=255)
                                aCGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaata                                                               >  1:127017/1‑100 (MQ=255)
                                aCGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaata                                                               <  1:2110380/100‑1 (MQ=255)
                                 cGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaata                                                               >  2:1432804/1‑99 (MQ=255)
                                  gCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAcc                                                             <  2:127018/100‑1 (MQ=255)
                                   cttccgatctgAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACCGAATAAATAccc                                                            <  1:1509812/90‑1 (MQ=255)
                                    ccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAcccc                                                           >  2:841662/1‑100 (MQ=255)
                                     cACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCg                                                          >  1:1558044/1‑100 (MQ=255)
                                     cACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCg                                                          >  2:1566280/1‑100 (MQ=255)
                                         gAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGa                                                      <  2:1911458/100‑1 (MQ=255)
                                           gCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGata                                                    >  1:354824/1‑100 (MQ=255)
                                            cGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataa                                                   >  1:1526337/1‑100 (MQ=255)
                                               aTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAatat                                                <  1:986934/100‑1 (MQ=255)
                                                 cGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCt                                              <  2:325811/100‑1 (MQ=255)
                                                 cGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCt                                              <  1:1894544/100‑1 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCACACTGAATAAATACCCCGGTGATAATATCtt                                             >  1:293673/1‑100 (MQ=255)
                                                  ctcATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                             <  2:1480004/98‑1 (MQ=255)
                                                  ccgatctgTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATGATATCtt                                             <  1:1212426/93‑1 (MQ=255)
                                                    cATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCttag                                           >  1:1480008/1‑98 (MQ=255)
                                                    cATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                           <  2:1224083/100‑1 (MQ=255)
                                                    cATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                           <  2:911070/100‑1 (MQ=255)
                                                      tCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAc                                         <  1:1832314/100‑1 (MQ=255)
                                                       ccGTTCTACTGTGCCGGTCAACGGCTCGTTAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                           <  2:1959340/97‑1 (MQ=255)
                                                       ccGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAcc                                        >  2:1189003/1‑100 (MQ=255)
                                                       ccGTTCCAATGTGCCGGTCAACGACCCGATAGACACCAAATATCCAGTGTTCATGCCGTTATCAACCATAATAAATACCCCGTTGATAATATCTTTCAcc                                        >  2:188267/1‑100 (MQ=255)
                                                         gTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAg                                      >  1:637117/1‑100 (MQ=255)
                                                         gTTCCACTGTGACGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACGGAATAAATACCCCGGTGATAATATCttagatcgg                                      >  2:1212423/1‑93 (MQ=255)
                                                             cACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGtt                                  >  1:609269/1‑100 (MQ=255)
                                                              aCTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGttt                                 >  2:1073880/1‑100 (MQ=255)
                                                                tgtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGc                               <  1:841664/100‑1 (MQ=255)
                                                                 ntgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                              >  2:168734/2‑100 (MQ=255)
                                                                 gtgCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                              <  1:2077029/100‑1 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                            >  1:1445439/1‑100 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                            >  1:62809/1‑100 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                            >  1:599082/1‑100 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                            >  1:188156/1‑100 (MQ=255)
                                                                     cGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAcc                          >  2:2037166/1‑100 (MQ=255)
                                                                         cAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAc                      <  1:628790/100‑1 (MQ=255)
                                                                          aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                     <  2:637116/100‑1 (MQ=255)
                                                                          aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                     >  2:124897/1‑100 (MQ=255)
                                                                          aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                     >  2:497545/1‑100 (MQ=255)
                                                                           aCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACga                    <  2:62809/100‑1 (MQ=255)
                                                                              gCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                 >  2:1222383/1‑100 (MQ=255)
                                                                                ccGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAg               <  2:520463/100‑1 (MQ=255)
                                                                                 cGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGc              <  1:1222386/100‑1 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCTAACGAGATAGcc             <  1:1509704/100‑1 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc             >  1:1224738/1‑100 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCATTATCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc             >  2:1509700/1‑100 (MQ=255)
                                                                                   aTAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCa            <  1:168733/100‑1 (MQ=255)
                                                                                   aTAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCa            <  2:1781320/100‑1 (MQ=255)
                                                                                     aGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGa          <  2:293674/100‑1 (MQ=255)
                                                                                       tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc        <  1:1432808/100‑1 (MQ=255)
                                                                                       tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc        <  1:1566284/100‑1 (MQ=255)
                                                                                        cACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCt       >  1:1594003/1‑100 (MQ=255)
                                                                                         aCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTa      >  2:883312/1‑100 (MQ=255)
                                                                                           cAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAt    >  1:1168468/1‑100 (MQ=255)
                                                                                             aaTATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGc  >  2:2000280/1‑100 (MQ=255)
                                                                                                |                                                                                                
GAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGC  >  NC_000913/843653‑843845

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: