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.6% 781.1 / ‑3.9 256intergenic (‑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 (115/140);  total (116/140)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 282 total aligned reads displayed.

GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106
                                                                                             |                                                                                       
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  1:6521418/93‑1 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  1:7394702/93‑1 (MQ=255)
 gCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAg                                                                                        >  1:6722570/1‑93 (MQ=255)
 gCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAg                                                                                        <  2:5085217/93‑1 (MQ=255)
 gCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAg                                                                                        >  2:575444/1‑93 (MQ=255)
   aTTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAggc                                                                                      <  2:6822824/93‑3 (MQ=255)
        ttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTg                                                                                  >  2:4959850/1‑93 (MQ=255)
        ttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTg                                                                                  <  1:7073989/93‑1 (MQ=255)
         tAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTgg                                                                                 <  1:6962021/93‑2 (MQ=255)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCATGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            >  1:4717845/1‑93 (MQ=25)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            >  1:1317271/1‑93 (MQ=255)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            <  2:4375483/93‑1 (MQ=37)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            >  1:7235857/1‑93 (MQ=255)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            <  1:389030/93‑1 (MQ=25)
                  tGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaa                                                                        <  1:2279439/93‑1 (MQ=255)
                   gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaa                                                                       >  1:2795456/1‑93 (MQ=255)
                   gcgcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataaaa                                                                       >  1:2686668/1‑93 (MQ=255)
                      cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                    >  1:6063818/1‑93 (MQ=255)
                      cTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                                    <  1:2540144/93‑1 (MQ=25)
                       tGAAGTGTGGGTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                                   <  1:5526631/93‑1 (MQ=25)
                        gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                  >  2:72345/1‑93 (MQ=255)
                         aaGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAttt                                                                 >  2:6456645/1‑93 (MQ=255)
                          aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTa                                                                <  2:2986372/93‑1 (MQ=37)
                          aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTa                                                                <  2:2412404/93‑1 (MQ=37)
                          aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTa                                                                <  1:953890/93‑1 (MQ=37)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  1:3943126/93‑1 (MQ=37)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              >  2:5856580/1‑93 (MQ=255)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  1:4888797/93‑1 (MQ=37)
                             gtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACAt                                                             <  2:1148768/93‑1 (MQ=37)
                              tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                            >  2:4118500/1‑93 (MQ=255)
                              tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                            <  1:2625773/93‑1 (MQ=37)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  2:144805/93‑1 (MQ=37)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  1:6140167/93‑1 (MQ=37)
                                  tGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCg                                                        >  2:1964527/1‑93 (MQ=255)
                                   ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                       <  1:1833650/93‑1 (MQ=38)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      <  2:3255453/93‑1 (MQ=38)
                                     aTAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                                     <  2:1822658/93‑1 (MQ=38)
                                      tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                    <  2:3481514/93‑1 (MQ=38)
                                      tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTc                                                    <  1:3705622/93‑1 (MQ=38)
                                           aaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatg                                               >  1:1074025/1‑93 (MQ=255)
                                            aGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatga                                              <  1:3394792/93‑1 (MQ=37)
                                             gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                             >  1:6280173/1‑93 (MQ=255)
                                             gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                             <  1:2980143/93‑1 (MQ=37)
                                                aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                          <  2:161979/93‑1 (MQ=37)
                                                  tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtct                                        <  1:5803134/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:5930810/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:4169650/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:2438403/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  1:2584909/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      >  2:2769393/1‑93 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     >  2:3171383/1‑93 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     >  1:3452702/1‑93 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     >  2:3843282/1‑93 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     >  2:6539760/1‑93 (MQ=255)
                                                         aGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTc                                 >  2:5379182/1‑93 (MQ=255)
                                                         aGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTc                                 >  1:5592129/1‑93 (MQ=255)
                                                          ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTca                                <  2:2304747/93‑1 (MQ=21)
                                                            aaaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcaca                              >  2:2179800/1‑93 (MQ=21)
                                                             aaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATCTCTTCTGCTCACAt                             <  2:5190939/93‑1 (MQ=21)
                                                             aaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACAt                             >  2:1216520/1‑93 (MQ=18)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            <  2:3354421/93‑1 (MQ=21)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            <  2:6995390/93‑1 (MQ=21)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            <  1:217909/93‑1 (MQ=21)
                                                                aCGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaa                          <  2:447987/93‑1 (MQ=18)
                                                                 cGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaa                         >  1:4848075/1‑93 (MQ=255)
                                                                   ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                       <  1:3338456/93‑1 (MQ=18)
                                                                    gCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc                      >  1:7179431/1‑93 (MQ=17)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  1:1468588/1‑93 (MQ=18)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  2:5493746/1‑93 (MQ=18)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  2:7050942/1‑93 (MQ=18)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  2:7361811/1‑93 (MQ=18)
                                                                        aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                  >  1:425155/1‑93 (MQ=17)
                                                                        aTTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacac                  <  2:1264884/93‑1 (MQ=17)
                                                                         ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                 >  2:6155391/1‑93 (MQ=17)
                                                                         ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                 <  1:5373711/93‑1 (MQ=17)
                                                                         ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                 <  1:7526083/93‑1 (MQ=17)
                                                                          tcttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACaa                <  2:5429850/93‑1 (MQ=17)
                                                                          tcttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACaa                >  1:4052465/1‑93 (MQ=17)
                                                                           cttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAAt               >  1:1132893/1‑93 (MQ=17)
                                                                           cttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAAt               >  1:4991610/1‑93 (MQ=14)
                                                                              ttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAc            >  2:3586749/1‑93 (MQ=255)
                                                                              ttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAc            <  2:2596576/93‑1 (MQ=255)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           >  2:4873352/1‑93 (MQ=17)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  2:7490709/93‑1 (MQ=17)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  2:3753267/93‑1 (MQ=17)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  2:2877816/93‑1 (MQ=255)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  2:1683969/93‑1 (MQ=255)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          >  1:2040991/1‑93 (MQ=17)
                                                                                 tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt         <  2:4534698/93‑1 (MQ=17)
                                                                                 tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt         >  2:1579652/1‑93 (MQ=17)
                                                                                   tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga       <  1:166868/93‑1 (MQ=255)
                                                                                   tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga       <  1:6879974/93‑1 (MQ=255)
                                                                                   tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTga       <  1:6092420/93‑1 (MQ=255)
                                                                                    cGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgag      <  1:4577329/93‑1 (MQ=255)
                                                                                     gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga     >  2:798739/1‑93 (MQ=18)
                                                                                     gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga     >  2:6209044/1‑93 (MQ=18)
                                                                                     gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga     >  2:6403739/1‑93 (MQ=18)
                                                                                     gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgaga     <  2:3972283/93‑1 (MQ=255)
                                                                                       cAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATc   >  1:5838953/1‑93 (MQ=18)
                                                                                        aaTGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCg  >  2:4700834/1‑93 (MQ=18)
                                                                                             |                                                                                       
GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106

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