Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,588 G→C 15.6% intergenic (‑79/+28) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5880GC15.6% 517.0 / 84.3 192intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (91/71);  new base C (22/8);  total (113/79)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 210 total aligned reads displayed.

CCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCC  >  NC_000913/3439489‑3439687
                                                                                                   |                                                                                                   
 cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                    <  2:2000762/100‑1 (MQ=255)
 cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                    >  2:522702/1‑100 (MQ=255)
 cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                    >  1:3215855/1‑100 (MQ=255)
 cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                    <  1:2013513/100‑1 (MQ=255)
 cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                    <  1:2001913/100‑1 (MQ=255)
   cTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGCATAAAAAAACCCGCCgggg                                                                                                  <  1:3928284/100‑1 (MQ=255)
     ggTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                                >  2:1241264/1‑100 (MQ=255)
       tCAAGTAACCCGAGAAATACTGGCGGGGATGAACAAAAAGAGTACAGCTATTGTATATAAGGGGGGATACTACGGGTATAAAGAAACCCGCCGGGGCgg                                                                                               <  2:1906569/99‑1 (MQ=255)
        cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                              >  2:3395812/1‑99 (MQ=255)
          aggaaCCACATAAATACTCCAGGGGATGAACAAAAAGAGCCCAGCTATTGTAGATCAGTGTGGATACTACGGGTATAAAAAAACCCGCCGGGGCgggggt                                                                                           >  2:3714806/4‑97 (MQ=255)
          aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtt                                                                                            >  1:3807798/1‑99 (MQ=255)
             aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                        <  1:2838982/100‑1 (MQ=255)
             aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtttttt                                                                                        <  1:536671/100‑1 (MQ=255)
              aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttttt                                                                                       <  1:3338793/100‑1 (MQ=255)
                cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                     >  2:2022579/1‑100 (MQ=255)
                cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                     <  1:2409676/100‑1 (MQ=255)
                   aTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGtt                                                                                  <  2:3231492/100‑1 (MQ=255)
                     aaaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTACATAAGTGGGGATACTACGGGTGTAAAAAAACCCGCCCCGGCGGGTTTTTTTatacccg                                                                                <  2:2199034/100‑7 (MQ=25)
                       aTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                              >  1:2763240/1‑100 (MQ=255)
                       aTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                              <  1:3711200/100‑1 (MQ=255)
                         aCTCCAGAGGATGAACACAAAGAGAACAGCTATAGTAGCTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTc                                                                             <  2:3320683/99‑1 (MQ=255)
                          tcttccgatctGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg                                                                           <  1:3480802/90‑1 (MQ=255)
                          cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCa                                                                            >  1:178260/1‑99 (MQ=255)
                          cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg                                                                           >  1:1659649/1‑100 (MQ=255)
                          cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAg                                                                           >  1:165031/1‑100 (MQ=255)
                           tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGATTTTTTACGTTGCTTCAGa                                                                          >  2:3741391/1‑100 (MQ=255)
                            ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAt                                                                         >  1:339479/1‑100 (MQ=255)
                              aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                       >  2:1570020/1‑100 (MQ=255)
                              aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                       >  2:266080/1‑100 (MQ=255)
                              aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                       >  1:107671/1‑100 (MQ=255)
                               ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAc                                                                      >  2:2270325/1‑100 (MQ=255)
                                gggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                     >  2:260530/1‑100 (MQ=255)
                                 ggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActc                                                                    >  2:2295646/1‑100 (MQ=255)
                                  gATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActct                                                                   >  1:2531048/1‑100 (MQ=255)
                                    tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCTGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                                 <  2:1110562/100‑1 (MQ=255)
                                    tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                                 >  2:3065390/1‑100 (MQ=255)
                                    tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                                 >  1:80879/1‑100 (MQ=255)
                                    tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTCTTTTTACGTTGATTCAAAttccact                                                                   >  2:1845121/1‑93 (MQ=255)
                                     gAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGTTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGagatcggaaga                                                                >  2:3480776/1‑90 (MQ=255)
                                      aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                               <  2:4000261/100‑1 (MQ=255)
                                      aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                               <  2:2816278/100‑1 (MQ=255)
                                      aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                               <  1:2980195/100‑1 (MQ=255)
                                       aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagc                                                              >  1:1056332/1‑100 (MQ=255)
                                       aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGAGACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagc                                                              <  2:2930126/100‑1 (MQ=255)
                                        cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagca                                                             <  1:1293395/100‑1 (MQ=255)
                                        cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagc                                                              >  1:1708399/1‑99 (MQ=255)
                                        cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTACAGATTACTCTgcagca                                                             <  2:3282085/100‑1 (MQ=255)
                                         aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGcagcag                                                            >  2:2511625/1‑100 (MQ=255)
                                           aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCt                                                          <  2:1607865/100‑1 (MQ=255)
                                             agagTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                        <  1:3066468/100‑1 (MQ=255)
                                             agagTACAGCTATTGTAGACAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTATGTTGCTGCAGAGTAATCTgaagcaacg                                                          <  2:1937974/98‑4 (MQ=38)
                                                gTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGAGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGc                                                     >  2:4105923/1‑100 (MQ=255)
                                                  aCAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCGGATTACTCTGCAGCAGCTTCTGCtt                                                   >  1:2189111/1‑100 (MQ=255)
                                                   cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                  >  1:230685/1‑100 (MQ=255)
                                                   cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                  >  2:795321/1‑100 (MQ=255)
                                                   cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                  >  1:685566/1‑100 (MQ=255)
                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtc                                                 <  1:2636419/100‑1 (MQ=255)
                                                     gCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGGGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCt                                                    >  2:2940151/1‑96 (MQ=255)
                                                      cTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                               <  2:1415544/100‑1 (MQ=255)
                                                      cTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtct                                                >  2:562934/1‑99 (MQ=255)
                                                        aTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                             <  1:2349527/100‑1 (MQ=255)
                                                           gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGTGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          <  1:1387053/100‑1 (MQ=255)
                                                           gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                          >  1:1105904/1‑100 (MQ=255)
                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCTGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACg                                         <  2:2772842/100‑1 (MQ=255)
                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACg                                         >  1:1220355/1‑100 (MQ=255)
                                                              gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGAt                                       >  1:487043/1‑100 (MQ=255)
                                                              gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGAt                                       >  1:2953973/1‑100 (MQ=255)
                                                               aTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                      <  1:3954170/100‑1 (MQ=255)
                                                                tAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCa                                     <  2:2149703/100‑1 (MQ=255)
                                                                 aaGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCaa                                    <  1:3036831/100‑1 (MQ=255)
                                                                   gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                  >  1:401096/1‑100 (MQ=255)
                                                                     ggggATACTACGGGTATAAAGAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAg                                >  2:2544434/1‑100 (MQ=255)
                                                                      gggATACTACGGGTATAAAAAAGCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGc                               >  2:2411779/1‑100 (MQ=255)
                                                                       ggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCt                              >  1:1806597/1‑100 (MQ=255)
                                                                       ggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCt                              >  1:3159376/1‑100 (MQ=255)
                                                                                gagTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                     <  2:3383775/98‑1 (MQ=255)
                                                                                  acgtaAAAAAACCCGCCCCGGCGGGCTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                   >  1:2248759/4‑100 (MQ=38)
                                                                                   cgtaAAAAAACCCGCCCTGGCGGGATTTTCTAAGTTGCTTCAAATTCCTCTGCAGAAGCTTCTGCTTTCTCTGAACGATGAACCAGCTCGATGTAAGcc                   >  2:3254203/3‑99 (MQ=255)
                                                                                   cgtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCa                  >  2:2973498/3‑100 (MQ=255)
                                                                                    gtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAt                 >  1:416783/2‑100 (MQ=255)
                                                                                    gtaAAAAAACCCGCCCCGGCGGGTTTGTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAGGATCAACCAGCTCGATGTAAGCCAt                 >  1:2497708/2‑100 (MQ=38)
                                                                                      aaaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg               >  1:3461515/1‑100 (MQ=255)
                                                                                      aaaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg               >  1:2100844/1‑100 (MQ=255)
                                                                                      aaaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg               >  1:3590142/1‑100 (MQ=255)
                                                                                       aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg              >  1:738060/1‑100 (MQ=255)
                                                                                       aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg              >  2:2137673/1‑100 (MQ=255)
                                                                                       aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg              >  2:3599410/1‑100 (MQ=255)
                                                                                       aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg              >  1:969604/1‑100 (MQ=255)
                                                                                        aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc             >  2:2051503/1‑100 (MQ=255)
                                                                                        aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc             >  1:1900951/1‑100 (MQ=255)
                                                                                        aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc             >  2:496721/1‑100 (MQ=255)
                                                                                         aaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcg            >  2:620550/1‑100 (MQ=255)
                                                                                         aaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcg            <  2:3435396/100‑1 (MQ=255)
                                                                                           aaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg          <  2:1548507/100‑1 (MQ=255)
                                                                                            aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGt         <  1:3831501/100‑1 (MQ=255)
                                                                                            aCCCGCCGGGGCGGGTTTTTTTACGTTGCTCCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGCACAATCAGCTCGATGTACGCCATCGGCGCGt         >  1:1848966/1‑100 (MQ=255)
                                                                                              tcgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTATGCAGCAGCTTCTGCTTTCTCTGAACGATCAACGAGCTCGATGTAAGTCATCGGCGCGTTg       >  2:1715057/2‑100 (MQ=255)
                                                                                                 ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCg    <  2:3845074/100‑1 (MQ=255)
                                                                                                 ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCg    >  1:2867036/1‑100 (MQ=255)
                                                                                                   ggggCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc  <  1:2042290/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCC  >  NC_000913/3439489‑3439687

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