Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 335,166 T→A 17.0% intergenic (‑144/‑114) yahC ← / → yahD putative inner membrane protein/ankyrin repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913335,1660TA17.0% 493.3 / 79.0 200intergenic (‑144/‑114)yahC/yahDputative inner membrane protein/ankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base T (90/75);  new base A (14/20);  total (105/95)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 212 total aligned reads displayed.

TGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAA  >  NC_000913/335068‑335265
                                                                                                  |                                                                                                   
tgCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtatt                                                                                                    <  1:2828672/100‑1 (MQ=255)
tgCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtatt                                                                                                    >  1:2774941/1‑100 (MQ=255)
tgCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtatt                                                                                                    >  2:3507997/1‑100 (MQ=255)
  cAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtattat                                                                                                  <  1:1986399/100‑1 (MQ=255)
    gACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTc                                                                                                >  2:2471209/1‑100 (MQ=255)
    gACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTc                                                                                                <  2:895566/100‑1 (MQ=255)
    gACTATAGAGAAATAAAACCTGATTGAGGGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTc                                                                                                >  2:2511741/1‑100 (MQ=255)
      cTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaa                                                                                              <  2:203995/100‑1 (MQ=255)
        ataGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaaca                                                                                            <  1:3419876/100‑1 (MQ=255)
          agagAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaat                                                                                          <  1:2053196/100‑2 (MQ=255)
          agagAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaat                                                                                          <  2:2953119/100‑2 (MQ=255)
           gagaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaa                                                                                         >  2:3983963/1‑100 (MQ=255)
           gagaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaa                                                                                         >  1:124031/1‑100 (MQ=255)
           gagaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatc                                                                                         <  2:4216588/100‑3 (MQ=255)
            agaAATAAAACCTGATTTAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatct                                                                                        <  1:4288684/100‑4 (MQ=255)
              aaaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttt                                                                                      <  2:1715948/100‑6 (MQ=255)
               aaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAg                                                                                       <  1:3493244/98‑1 (MQ=255)
                aTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCg                                                                                    >  2:2099464/1‑100 (MQ=255)
                aTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCg                                                                                    >  2:2842738/1‑100 (MQ=255)
                 tAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCga                                                                                   <  2:1088968/100‑1 (MQ=255)
                 tAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttttcg                                                                                   <  2:2171056/100‑9 (MQ=38)
                 tAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttttc                                                                                    <  2:811136/99‑8 (MQ=39)
                 tAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAACAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCga                                                                                   >  1:731144/1‑100 (MQ=255)
                  aaaaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAGCAAAAAGGCgag                                                                                  <  2:2889873/100‑1 (MQ=255)
                  aaaaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCgag                                                                                  <  2:4075373/100‑1 (MQ=255)
                  aaaaCCTGATTGAGTCAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttttcgc                                                                                  <  2:4342771/100‑10 (MQ=38)
                    aaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGAt                                                                                >  1:931361/1‑100 (MQ=255)
                        tGATTGAGTGAAATGGTATCGGGTCAAAGAAAGAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGTCGAGATACgg                                                                            <  1:2868443/100‑1 (MQ=255)
                             gAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCGACAAAAAGGCGAGATACGGCGAAt                                                                       >  2:2741814/1‑100 (MQ=255)
                             gAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAt                                                                       <  2:1176542/100‑1 (MQ=255)
                             gAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAt                                                                       >  1:3625025/1‑100 (MQ=255)
                             gAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAt                                                                       >  2:3409229/1‑100 (MQ=255)
                             gAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCACACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAt                                                                       <  1:3847527/100‑1 (MQ=255)
                              aGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAtt                                                                      <  2:2700975/100‑1 (MQ=255)
                              aGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAtt                                                                      <  2:1477782/100‑1 (MQ=255)
                                  aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAt                                                                  <  1:3517548/100‑1 (MQ=255)
                                  aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAt                                                                  <  1:3919267/100‑1 (MQ=255)
                                  aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAt                                                                  <  2:1291787/100‑1 (MQ=255)
                                  aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAt                                                                  <  2:1229507/100‑1 (MQ=255)
                                   aaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAtt                                                                 >  2:2343807/1‑100 (MQ=255)
                                     gttgaTCGTGTCAAAGATACAATTCTTGCACTATCTCACCCTGCTAAAAGATTGTTGACTATTATTCCACCAAAAGGCGAGATACGGCGAATTaaan                                                                  <  1:1975139/93‑2 (MQ=255)
                                     gggTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTa                                                               >  2:3046331/1‑100 (MQ=255)
                                      ggTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTaa                                                              >  2:2544276/1‑100 (MQ=255)
                                       gTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAg                                                             >  1:1136755/1‑100 (MQ=255)
                                         aTCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCg                                                           <  1:3073790/100‑1 (MQ=255)
                                         aTCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCg                                                           >  2:942129/1‑100 (MQ=255)
                                          tCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCgg                                                          <  2:949353/100‑1 (MQ=255)
                                          tCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCgg                                                          <  1:415897/100‑1 (MQ=255)
                                          tCGGGTCAAAGAAACAAATAGTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCgg                                                          <  1:1485318/100‑1 (MQ=255)
                                           cGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCggg                                                         <  1:4190927/100‑1 (MQ=255)
                                            gggTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCggg                                                         >  1:2011543/1‑99 (MQ=255)
                                            gggTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGt                                                        >  1:866485/1‑100 (MQ=255)
                                             ggTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGtt                                                       >  1:2263080/1‑100 (MQ=255)
                                             ggTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGtt                                                       >  2:3748377/1‑100 (MQ=255)
                                             ggTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGt                                                        <  2:1758449/99‑1 (MQ=255)
                                                 aaaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCttt                                                   >  2:4187736/1‑100 (MQ=255)
                                                 aaaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCttt                                                   <  2:1984698/100‑1 (MQ=255)
                                                  aaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtttt                                                  >  1:1503532/1‑100 (MQ=255)
                                                   aGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTa                                                 <  1:765491/100‑1 (MQ=255)
                                                   aGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTa                                                 <  1:3534367/100‑1 (MQ=255)
                                                      aaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCa                                              <  1:2544263/100‑1 (MQ=255)
                                                      aaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCa                                              >  1:2954848/1‑100 (MQ=255)
                                                      aaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCa                                              <  1:1003594/100‑1 (MQ=255)
                                                      aaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCa                                              >  1:3304303/1‑100 (MQ=255)
                                                      aaCAAATATTGCACTACCGCACACTGCCAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCa                                              <  1:4033408/100‑1 (MQ=255)
                                                        cAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGa                                            >  1:1398773/1‑100 (MQ=255)
                                                           atatTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTa                                         >  2:3177199/1‑100 (MQ=255)
                                                             atTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATc                                       >  1:603926/1‑100 (MQ=255)
                                                                gCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGtt                                    >  1:1633200/1‑100 (MQ=255)
                                                                  aCTACCGAACACTGCGAAAAGATTGTTGAAAATTATTCAACAAAAAGGAGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  <  2:648236/100‑1 (MQ=255)
                                                                     aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTAt                               >  1:3865069/1‑100 (MQ=255)
                                                                     aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTAt                               >  1:3353625/1‑100 (MQ=255)
                                                                       cGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATtc                             <  1:1054982/100‑1 (MQ=255)
                                                                         cacaCTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATtctc                           >  1:2535539/1‑100 (MQ=255)
                                                                         cacaCTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATtctc                           <  2:382549/100‑1 (MQ=255)
                                                                            acTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAAt                        <  2:1423666/100‑1 (MQ=255)
                                                                                cGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCAt                    >  1:792052/1‑100 (MQ=255)
                                                                                 gAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATc                   <  1:802040/100‑1 (MQ=255)
                                                                                 aaaaaaaTTGTTTAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATc                   <  1:2848296/99‑1 (MQ=255)
                                                                                  gcctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATc                   >  1:178474/7‑99 (MQ=255)
                                                                                  gcctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  >  2:3324571/7‑100 (MQ=255)
                                                                                  gcctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  >  2:520151/7‑100 (MQ=255)
                                                                                  aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  <  1:2629082/100‑1 (MQ=255)
                                                                                  aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  <  2:348026/100‑1 (MQ=255)
                                                                                  aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  <  1:826623/100‑1 (MQ=255)
                                                                                   cctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGa                 >  2:371630/6‑100 (MQ=255)
                                                                                   cctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGa                 >  1:2335775/6‑100 (MQ=255)
                                                                                    ctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATTAGATTATCGTTGCTATTCTCAATTCATCGa                 >  2:1812876/5‑99 (MQ=39)
                                                                                     aGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGaaa               >  2:460239/1‑100 (MQ=255)
                                                                                     aGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATAAAATTTAAGCGGGATCTTTTATCAGATTATCGTTGCTATTCTTAACTCATCGaaa               >  1:4121486/1‑100 (MQ=255)
                                                                                       atggttgAATATTATTTAACAGAAAGGCGTGATACGGCGAATTAAATTTTACTGGGTTCTTTTATGAGATTATCGTGGCTATTCTCAAtt                       <  1:2718044/87‑1 (MQ=255)
                                                                                         tgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAg           >  1:1613713/1‑100 (MQ=255)
                                                                                          gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  2:2733846/1‑100 (MQ=255)
                                                                                          gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  1:5477/1‑100 (MQ=255)
                                                                                            tgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCca        >  1:998018/1‑100 (MQ=255)
                                                                                             gAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCcac       >  2:1645198/1‑100 (MQ=255)
                                                                                             gAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCcac       <  1:2143924/100‑1 (MQ=255)
                                                                                                tattatTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaaa    <  2:2997727/100‑1 (MQ=255)
                                                                                                 attattCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaaaa   <  1:953437/100‑1 (MQ=255)
                                                                                                  ttattCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaaaaa  <  2:499798/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAA  >  NC_000913/335068‑335265

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: