Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,781,027 C→G intergenic (‑151/+58) ypjA ← / ← b2648 putative ATP‑binding component of a transport system/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,781,0270CG99.6% 735.5 / ‑4.8 239intergenic (‑151/+58)ypjA/b2648putative ATP‑binding component of a transport system/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (96/142);  minor base A (0/1);  total (96/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

Only 100 of 240 total aligned reads displayed.

AACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTT  >  NC_000913/2780938‑2781116
                                                                                          |                                                                                         
aaCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCa                                                                                         >  2:3132607/1‑93 (MQ=255)
  cGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaa                                                                                       <  2:6229456/93‑1 (MQ=255)
  cGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaa                                                                                       >  2:6435891/1‑93 (MQ=255)
     aGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaata                                                                                    <  1:2916585/93‑1 (MQ=25)
       tGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaa                                                                                  >  1:4379800/1‑93 (MQ=255)
        gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaa                                                                                 >  2:1248437/1‑93 (MQ=255)
        gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaa                                                                                 <  1:3299718/93‑1 (MQ=25)
        gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaa                                                                                 >  2:4865766/1‑93 (MQ=255)
        gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaa                                                                                 >  1:1614085/1‑93 (MQ=255)
           cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                              >  1:4113680/1‑93 (MQ=255)
            tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                             <  1:1907504/93‑1 (MQ=255)
            tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                             <  1:589578/93‑1 (MQ=25)
            tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                             >  1:3767043/1‑93 (MQ=255)
             gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                            <  1:1004084/93‑1 (MQ=25)
             gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                            <  1:973391/93‑1 (MQ=25)
             gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                            <  2:6100125/93‑1 (MQ=255)
             gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                            <  1:3932933/93‑1 (MQ=255)
              aaGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAttt                                                                           <  1:1149593/93‑1 (MQ=255)
                gtgtGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTAc                                                                         >  1:4029242/1‑93 (MQ=255)
                 tgtgGTTGGATAAAGAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                                        <  2:1123109/93‑1 (MQ=37)
                 tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                                        <  2:3834022/93‑1 (MQ=255)
                 tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                                        <  2:1542950/93‑1 (MQ=37)
                    ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATaa                                                                     >  2:5718584/1‑93 (MQ=255)
                     gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                                    <  2:3405542/93‑1 (MQ=255)
                     gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                                    <  1:3357377/93‑1 (MQ=37)
                     gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                                    >  1:4320287/1‑93 (MQ=255)
                        ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                                 >  1:5065871/1‑93 (MQ=255)
                         gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                                <  1:3902974/93‑1 (MQ=255)
                          aTAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                                               <  2:3863824/93‑1 (MQ=255)
                           tAAAAAAGTCAATCCATTCAGGAAGTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                              <  2:2687987/93‑1 (MQ=37)
                           tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                              <  1:895707/93‑1 (MQ=38)
                           tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                              <  1:4641947/93‑1 (MQ=37)
                             aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtg                                                            <  1:2463845/93‑1 (MQ=38)
                             aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtg                                                            <  1:1773831/93‑1 (MQ=255)
                              aaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtga                                                           >  1:4422503/1‑93 (MQ=255)
                              aaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtga                                                           <  1:5461983/93‑1 (MQ=255)
                              aaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtga                                                           <  2:5812548/93‑1 (MQ=38)
                               aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                          <  1:4593481/93‑1 (MQ=37)
                                 aGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatga                                                        >  1:5642977/1‑93 (MQ=255)
                                  gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                                       <  1:521319/93‑1 (MQ=37)
                                   tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                                      >  1:3693319/1‑93 (MQ=255)
                                   tCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                                      <  1:2116936/93‑1 (MQ=255)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatatt                                                   >  1:4674758/1‑92 (MQ=25)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   >  1:1486682/1‑93 (MQ=255)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   <  2:6218105/93‑1 (MQ=255)
                                      aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtc                                                   <  2:813989/93‑1 (MQ=255)
                                       tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtct                                                  <  1:652847/93‑1 (MQ=255)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                >  2:2998899/1‑93 (MQ=255)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  2:1143336/93‑1 (MQ=25)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  2:683271/93‑1 (MQ=255)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  2:1416484/93‑1 (MQ=25)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  2:6498246/93‑1 (MQ=25)
                                         cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                                <  2:4272698/93‑1 (MQ=255)
                                          aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                               >  2:5485751/1‑93 (MQ=255)
                                           ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                              <  1:1471357/93‑1 (MQ=255)
                                           ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                              >  2:1393215/1‑93 (MQ=255)
                                              aGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTc                                           >  2:3994541/1‑93 (MQ=255)
                                               ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTca                                          >  1:1141237/1‑93 (MQ=255)
                                               ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTca                                          >  1:358124/1‑93 (MQ=255)
                                                  aaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACAt                                       <  1:6661685/93‑1 (MQ=255)
                                                     aCGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaa                                    <  2:362051/93‑1 (MQ=255)
                                                      cGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaa                                   <  1:4757694/93‑1 (MQ=18)
                                                       gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                                  <  2:893720/93‑1 (MQ=255)
                                                       gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                                  <  2:1746665/93‑1 (MQ=255)
                                                       gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                                  >  1:765027/1‑93 (MQ=255)
                                                       gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                                  <  2:3471096/93‑1 (MQ=18)
                                                        ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 >  1:1923060/1‑93 (MQ=255)
                                                        ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 >  1:553403/1‑93 (MQ=17)
                                                        ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                                 >  2:3245326/1‑93 (MQ=18)
                                                         gCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc                                <  2:3087158/93‑1 (MQ=18)
                                                          cGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATca                               >  1:4949529/1‑93 (MQ=255)
                                                             aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                            >  1:2211724/1‑93 (MQ=17)
                                                             aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                            <  2:2570874/93‑1 (MQ=255)
                                                             aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                            <  1:2733274/93‑1 (MQ=17)
                                                              ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                           >  2:2364947/1‑93 (MQ=255)
                                                               tcttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACaa                          <  2:4493720/93‑1 (MQ=17)
                                                                 ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa                        <  2:466658/93‑1 (MQ=255)
                                                                  tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa                       <  1:6369040/93‑1 (MQ=255)
                                                                  tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa                       >  2:5245444/1‑93 (MQ=17)
                                                                   ttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAc                      <  1:5314831/93‑1 (MQ=255)
                                                                   ttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAc                      <  1:408041/93‑1 (MQ=255)
                                                                      tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt                   >  1:2643771/1‑93 (MQ=17)
                                                                      tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt                   >  2:5565891/1‑93 (MQ=17)
                                                                        tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga                 >  1:6085532/1‑93 (MQ=17)
                                                                        tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga                 <  2:1682473/93‑1 (MQ=255)
                                                                        tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga                 <  1:5801519/93‑1 (MQ=255)
                                                                          gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga               <  2:4982923/93‑1 (MQ=255)
                                                                          gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga               <  2:2629972/93‑1 (MQ=255)
                                                                           aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACACTAACTTTGAGAt              <  1:5345890/93‑1 (MQ=18)
                                                                           aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt              <  1:452418/93‑1 (MQ=255)
                                                                           aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt              >  1:244785/1‑93 (MQ=255)
                                                                           aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt              <  2:2149722/93‑1 (MQ=255)
                                                                           aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt              >  2:3172686/1‑93 (MQ=255)
                                                                            cAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATc             <  2:4209724/93‑1 (MQ=255)
                                                                             aaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTCGAGATCg            >  1:2001698/1‑93 (MQ=18)
                                                                                gaggCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAg         >  2:5172964/4‑93 (MQ=255)
                                                                                  ggCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGAt       >  1:170196/2‑93 (MQ=255)
                                                                                     cGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGt    <  2:4411104/93‑1 (MQ=255)
                                                                                     cGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGt    >  2:1623492/1‑93 (MQ=255)
                                                                                       ttGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGttt  >  1:4194679/1‑93 (MQ=255)
                                                                                          |                                                                                         
AACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTT  >  NC_000913/2780938‑2781116

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