Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,781,020 +G intergenic (‑144/+65) 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,0191.G99.5% 618.5 / ‑4.3 203intergenic (‑143/+66)ypjA/b2648putative ATP‑binding component of a transport system/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (88/114);  total (89/114)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 221 total aligned reads displayed.

GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106
                                                                                             |                                                                                        
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                           >  1:4402656/1‑93 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                           >  2:1888296/1‑93 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                           >  1:2886047/1‑93 (MQ=255)
agCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                           <  1:3233629/92‑1 (MQ=255)
 gCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAg                                                                                         <  1:531222/93‑1 (MQ=255)
 gCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAg                                                                                         <  2:1986644/93‑1 (MQ=255)
   aTTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAggc                                                                                       <  2:70769/93‑3 (MQ=255)
   aTTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAggc                                                                                       >  2:2347290/1‑91 (MQ=255)
       tttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGtt                                                                                   <  1:3598119/93‑1 (MQ=255)
        ttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTg                                                                                  <  2:3897028/93‑1 (MQ=255)
        ttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTg                                                                                  >  1:364942/1‑93 (MQ=255)
            aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaa                                                                              >  2:3032104/1‑93 (MQ=255)
            aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaa                                                                              >  1:2310004/1‑93 (MQ=255)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            <  1:4326093/93‑1 (MQ=255)
               cAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaat                                                                           >  2:2205892/1‑93 (MQ=255)
                aGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaata                                                                          <  2:807415/93‑1 (MQ=25)
                    cgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaat                                                                      <  1:1790589/93‑1 (MQ=25)
                    cgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaat                                                                      <  2:1385406/93‑1 (MQ=25)
                     gcTGATATGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAAAACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaacaa                                                                     <  2:4223352/93‑4 (MQ=25)
                     gcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATg                                                                     >  1:2116848/1‑93 (MQ=255)
                      cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCTACAATGAGGCCGTTGGCAAAATAAAATGa                                                                    <  2:2050832/93‑1 (MQ=25)
                      cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                    <  1:1179828/93‑1 (MQ=255)
                      cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                    <  1:1234513/93‑1 (MQ=255)
                       tGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                   >  1:1372146/1‑93 (MQ=255)
                          aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTa                                                                <  1:2778744/93‑1 (MQ=37)
                           gtgtGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTAc                                                               >  1:3840345/1‑93 (MQ=255)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  2:2962105/93‑1 (MQ=37)
                             gtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACAt                                                             <  1:1597768/93‑1 (MQ=37)
                              tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                            <  2:3450081/93‑1 (MQ=37)
                               ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATaa                                                           >  2:4250868/1‑93 (MQ=255)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          >  2:1248130/1‑93 (MQ=255)
                                 ttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                         <  2:4205787/93‑1 (MQ=255)
                                  tGGATAAAAAAGTGAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCg                                                        >  2:3640951/1‑93 (MQ=37)
                                   ggATAAAAAGGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                       >  2:3143336/1‑93 (MQ=37)
                                   ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                       <  1:1588184/93‑1 (MQ=37)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      <  2:4713098/93‑1 (MQ=37)
                                      tAAATAAGGCAATTCATTCAGGAAATATGGGCGTAGTCTTTTCTTTCGACAATGAGGCGGTGGGCAAAATAAAATGATTTACATAATCGTTTc                                                    <  1:806626/93‑1 (MQ=21)
                                      tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                    <  2:2079932/93‑1 (MQ=37)
                                       aaaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCt                                                   <  1:3106092/93‑1 (MQ=38)
                                        aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtg                                                  <  2:1843220/93‑1 (MQ=255)
                                         aaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCACAATAAAATGATTTACATAATCGTTTCtgt                                                 >  1:1083262/1‑92 (MQ=37)
                                         aaaaGACAATCCAATCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGTCCGCTGGCAAAATAAAATGATTTACATAATCGTTTCtga                                                 >  2:2717336/1‑93 (MQ=25)
                                          aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                <  1:570188/93‑1 (MQ=37)
                                             gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                             >  2:1148275/1‑93 (MQ=255)
                                                aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                          <  1:2013191/93‑1 (MQ=25)
                                                   ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                       <  1:1191951/93‑1 (MQ=255)
                                                   ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                       <  2:1652639/93‑1 (MQ=255)
                                                   ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                       >  2:1218564/1‑93 (MQ=255)
                                                   ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtctt                                       >  2:2290607/1‑93 (MQ=255)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:614793/93‑1 (MQ=255)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:3932286/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  2:4515246/93‑1 (MQ=25)
                                                      ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                    >  2:1541946/1‑93 (MQ=255)
                                                      ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                    <  1:989853/93‑1 (MQ=255)
                                                      ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                    >  1:222485/1‑93 (MQ=255)
                                                      ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                    >  1:1224167/1‑93 (MQ=255)
                                                       tCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGc                                   >  1:3454462/1‑93 (MQ=255)
                                                       tCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGATTCTGATGAATATCTTCTGc                                   >  2:3886958/1‑93 (MQ=21)
                                                        cAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCt                                  >  2:2321651/1‑93 (MQ=255)
                                                        cAGGAAAGACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCt                                  <  1:1265075/93‑1 (MQ=21)
                                                          ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTca                                <  1:2522512/93‑1 (MQ=21)
                                                           gAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcac                               >  2:3281744/1‑93 (MQ=255)
                                                            aaaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcaca                              <  1:2261603/93‑1 (MQ=21)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            >  1:2343350/1‑93 (MQ=18)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            >  1:2694153/1‑93 (MQ=18)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            >  2:1063024/1‑93 (MQ=255)
                                                               tACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaa                           >  1:843906/1‑93 (MQ=18)
                                                                  gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                        <  2:2145000/93‑1 (MQ=18)
                                                                  gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                        >  1:477634/1‑93 (MQ=255)
                                                                   ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                       <  1:3331867/93‑1 (MQ=255)
                                                                   ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                       <  2:856964/93‑1 (MQ=255)
                                                                    gCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc                      <  2:629110/93‑1 (MQ=255)
                                                                    gCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc                      <  1:1912103/93‑1 (MQ=255)
                                                                    gCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc                      <  1:4496875/93‑1 (MQ=255)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  1:2591196/1‑93 (MQ=18)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  1:4691508/1‑93 (MQ=255)
                                                                        aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                  <  1:3319376/93‑1 (MQ=255)
                                                                        aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                  >  2:2508123/1‑93 (MQ=18)
                                                                        aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                  <  2:2841564/93‑1 (MQ=255)
                                                                           cttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAAt               >  1:3205896/1‑93 (MQ=17)
                                                                           cttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAAt               <  1:2617099/93‑1 (MQ=255)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              >  2:1007287/1‑93 (MQ=17)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              <  2:3699169/93‑1 (MQ=255)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              >  1:2492128/1‑93 (MQ=255)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              >  2:2618467/1‑93 (MQ=17)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              >  1:3067313/1‑93 (MQ=17)
                                                                             tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa             <  2:1936396/93‑1 (MQ=255)
                                                                               tctttcGACAATGAGGCCGGTGGCAAAATAAAATGATTTACATAATCGTTTGTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           >  2:787237/1‑93 (MQ=14)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  1:1897116/93‑1 (MQ=255)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  2:4338097/93‑1 (MQ=255)
                                                                                 tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt         >  1:4045636/1‑93 (MQ=17)
                                                                                 tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt         >  1:1660320/1‑93 (MQ=17)
                                                                                   tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga       <  2:3983542/93‑1 (MQ=255)
                                                                                   tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga       <  1:1516781/93‑1 (MQ=255)
                                                                                     gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga     >  1:257878/1‑93 (MQ=18)
                                                                                     gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga     <  2:4569150/93‑1 (MQ=255)
                                                                                      aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt    >  2:3846325/1‑93 (MQ=17)
                                                                                      aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt    <  2:3275087/93‑1 (MQ=255)
                                                                                       gaaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATc   <  2:3221660/92‑1 (MQ=18)
                                                                                       cAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATc   >  2:4436364/1‑93 (MQ=18)
                                                                                             |                                                                                        
GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106

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