Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,806,042 (C)6→5 100% coding (42/759 nt) ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,0370C.77.1% 195.4 / 79.3 83coding (47/759 nt)ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base . (35/29);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01

CGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACT  >  NC_000913/1805939‑1806134
                                                                                                  |                                                                                                 
cGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcc                                                                                                  >  1:1338494/1‑100 (MQ=255)
 gTTCAGATTACCTTCGAAAGGTTTTTTTTTTTTTTTAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGTCGCAAACATGccc                                                                                                 >  1:1883955/1‑100 (MQ=255)
 gTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAATAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccc                                                                                                 >  1:1869794/1‑100 (MQ=255)
  ttCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcccc                                                                                                <  1:2478216/100‑1 (MQ=255)
      gATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcccccgcc                                                                                            <  1:1781760/100‑4 (MQ=255)
        ttACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCaa                                                                                         <  2:804104/100‑1 (MQ=255)
        ttACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCaa                                                                                         <  2:243496/100‑1 (MQ=255)
           ccTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACAt                                                                                      >  2:538344/1‑100 (MQ=255)
           ccTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACAt                                                                                      >  2:1275929/1‑100 (MQ=255)
           ccTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACAt                                                                                      >  1:655084/1‑100 (MQ=255)
               cGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACg                                                                                  >  2:2092073/1‑100 (MQ=255)
                gAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGa                                                                                 >  2:1662493/1‑100 (MQ=255)
                gAAAGGCTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGa                                                                                 >  1:1129851/1‑100 (MQ=255)
                 aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGAt                                                                                <  2:1207700/100‑1 (MQ=255)
                   aGGTTTTTTAGCGCTGGCAAGTTCATCCATGACAGTAAAACCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGAtn                                                                               <  2:1093710/99‑2 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGc                                                                             <  1:1857838/100‑1 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGc                                                                             >  2:628677/1‑100 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGc                                                                             >  1:1418590/1‑100 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGc                                                                             >  1:2036999/1‑100 (MQ=255)
                     gTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCa                                                                            <  1:2225085/100‑1 (MQ=255)
                     gTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCa                                                                            >  1:2380251/1‑100 (MQ=255)
                     gTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCa                                                                            <  1:81558/100‑1 (MQ=255)
                      ttttttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAg                                                                           <  2:1324700/100‑1 (MQ=255)
                      gatcttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAg                                                                           <  1:211345/96‑1 (MQ=255)
                        ttttAGCGCTGCCAGGATCATCCGTTACAGTAAAAACGGAACCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCa                                                                            <  1:125667/97‑1 (MQ=255)
                          ttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAAc                                                                       <  2:2052459/100‑1 (MQ=255)
                           tAGCGCTGGCAGGATCCACCATGACAGTCAAAACGGAATTAGCGACAGCCGCACTCAGTGACGCAAACATG‑CCCTCGCCAACATAACGAGTGCAg                                                                           >  1:2424647/1‑95 (MQ=255)
                           tAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAgagatc                                                                      >  2:211347/1‑95 (MQ=255)
                           tAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACt                                                                      <  1:1345460/100‑1 (MQ=255)
                             gcgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTg                                                                     <  2:1433002/100‑1 (MQ=255)
                              cgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTc                                                                   >  1:135412/1‑100 (MQ=255)
                              cgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTc                                                                   <  1:1275924/100‑1 (MQ=255)
                               gcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCt                                                                  <  2:1914597/100‑1 (MQ=255)
                                cTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCt                                                                  <  2:2582590/100‑1 (MQ=255)
                                cTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCtt                                                                 <  1:2411119/100‑1 (MQ=255)
                                 tGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTc                                                                <  2:301978/100‑1 (MQ=255)
                                  ggCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCCACATAACGATGGCAGGAACTGTCTTc                                                                >  1:2245505/1‑100 (MQ=255)
                                  ggCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTc                                                                <  2:2337089/100‑1 (MQ=255)
                                   gCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCa                                                               <  1:1242755/100‑1 (MQ=255)
                                     aGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGTCAACATAACGATGGCAGGAACTGTCTTCaaa                                                             <  1:555521/100‑1 (MQ=255)
                                     aGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCaaaa                                                            <  2:1686423/100‑1 (MQ=255)
                                       gATCATCCATTACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                          <  1:1631973/100‑1 (MQ=255)
                                        atcatcCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                          >  1:283200/1‑100 (MQ=255)
                                          catcCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTc                                                       <  2:1349098/100‑1 (MQ=255)
                                            tcCCTGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGTCCCCAGCCAACATAACGATGGCAGGAACAGTCGTCAAAAGCTTCa                                                      >  1:17707/1‑100 (MQ=255)
                                            tcCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAt                                                     <  1:672020/100‑1 (MQ=255)
                                             ccATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAt                                                     <  2:1238538/100‑1 (MQ=255)
                                              cATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAttt                                                   >  1:240234/1‑100 (MQ=255)
                                              cATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAttt                                                   <  1:1870873/100‑1 (MQ=255)
                                                tGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTa                                                  >  2:959367/1‑100 (MQ=255)
                                                tGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTAt                                                 >  1:1889351/1‑100 (MQ=255)
                                                tGACAGTAAAAACGGAACCATTGGCAGCTGCATTCACTGACGCAAACAAG‑CCCACGCCAACATAACGACGGCAGGAACTGTCTTCAAAAGCTTCATTTa                                                  <  1:2374477/99‑1 (MQ=255)
                                                  aCAGTAAATACGGAATTATCGGCAGCAGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGGACTGTCTTCAAAAGCTTCATTTATCa                                               >  1:561202/1‑100 (MQ=255)
                                                   cAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCaa                                              >  1:399988/1‑100 (MQ=255)
                                                    aGTAAAAACGGAATCATCGGCAGCTGTATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAACAGCTTCATTTATCAag                                             <  1:389745/100‑1 (MQ=255)
                                                     gTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAag                                             >  1:2576317/1‑100 (MQ=255)
                                                     gTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAaga                                            >  1:37123/1‑100 (MQ=255)
                                                       aaaaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAagag                                           >  2:844722/1‑99 (MQ=255)
                                                        aaaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTc                                         <  1:64180/100‑1 (MQ=255)
                                                         aaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTc                                         >  2:1621758/1‑100 (MQ=255)
                                                         aaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTcc                                        >  1:1616429/1‑100 (MQ=255)
                                                          aaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCg                                       >  2:674332/1‑100 (MQ=255)
                                                           aCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGt                                      >  2:936631/1‑100 (MQ=255)
                                                              gAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                    <  2:825915/100‑1 (MQ=255)
                                                              gAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTaca                                   >  2:2245018/1‑100 (MQ=255)
                                                                 tcatcGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTacaaca                                >  1:22512/1‑100 (MQ=255)
                                                                  catcGGCAGCGGCAGTCAGTGACGCAAACATG‑ACCCCGCCAACAAAACGAAGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                    <  1:6341/95‑1 (MQ=255)
                                                                   atcGGCAGCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACaaa                              >  1:684816/1‑100 (MQ=255)
                                                                          gCTGCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAagag                       >  2:2389588/1‑100 (MQ=255)
                                                                             gCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGAcc                    >  2:1930879/1‑100 (MQ=255)
                                                                             gCATTCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGAcc                    <  2:2364186/100‑1 (MQ=255)
                                                                                 tCAGTGACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATcg                >  2:150962/1‑100 (MQ=255)
                                                                                   aGTGACGCAAACATGCCCCCCGCCAACATAACGAGGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATcgc               <  2:2162889/100‑1 (MQ=255)
                                                                                    ttGAAGCAAACATG‑CCCCCGTCAACATATCGATGGCAGGAACTATCTTCAAAAGCTTCATTTATCAAGAGCCCGTACAACAAAAAAAGAGACCATCGCgg             <  1:2156823/99‑1 (MQ=255)
                                                                                      gACGCAAACATG‑CCCCCGCCAACCTAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGACCGTACAACAAAAAAAGAGACCATAGCGGTc           >  2:2058338/1‑100 (MQ=255)
                                                                                      gACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTc           >  1:2359333/1‑100 (MQ=255)
                                                                                      gACGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGACACCATCGCGGTc           >  1:958239/1‑100 (MQ=255)
                                                                                       aCGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTcc          <  2:70968/100‑1 (MQ=255)
                                                                                        cGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTcc          >  2:473183/1‑100 (MQ=255)
                                                                                        cGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTcc          >  2:1331197/1‑100 (MQ=255)
                                                                                        cGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTccc         >  1:1128960/1‑100 (MQ=255)
                                                                                        cGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTccc         >  1:630792/1‑100 (MQ=255)
                                                                                        cGCAAACATG‑CCCCCGCCAACATAACGATGGCAGGAAATGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTccc         >  1:220529/1‑100 (MQ=255)
                                                                                          cAAACATG‑CCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCgg       <  1:94623/100‑1 (MQ=255)
                                                                                           aaaCATG‑CCCCCGCCAACATAACGATGGCAGGAAATGTCTTCAAAAGCTTGATTTATCAGGAGTCCGTACATCAAAAAAAGAGACCATGGCGGTCCCGGa      >  1:1155902/1‑100 (MQ=255)
                                                                                              cATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaaa    >  1:2099634/1‑100 (MQ=255)
                                                                                              cATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaaa    >  2:1406243/1‑100 (MQ=255)
                                                                                                tGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACt  >  2:1610918/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
CGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACT  >  NC_000913/1805939‑1806134

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