Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,956,973:1 +A 15.4% coding (848/894 nt) ilvY ← transcriptional activator of ilvC; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,956,9731.A15.4% 543.5 / 110.4 169coding (848/894 nt)ilvYtranscriptional activator of ilvC; autorepressor
Reads supporting (aligned to +/- strand):  ref base . (73/70);  new base A (8/18);  total (81/88)
Fisher's exact test for biased strand distribution p-value = 8.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

Only 100 of 178 total aligned reads displayed.

AGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGC  >  NC_000913/3956874‑3957072
                                                                                                    |                                                                                                   
 gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑c                                                                                                    >  1:2069955/1‑100 (MQ=255)
 gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑c                                                                                                    >  2:1204250/1‑100 (MQ=255)
 gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGGTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCt                                                                                                      >  1:3681560/1‑99 (MQ=255)
  ggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑Ca                                                                                                   >  1:1889445/1‑100 (MQ=255)
   gAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CAt                                                                                                  >  1:2917887/1‑100 (MQ=255)
    aaaaaTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATg                                                                                                 <  2:2173036/100‑1 (MQ=255)
     aaaaTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGt                                                                                                >  1:1402510/1‑100 (MQ=255)
      aaaTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGATCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCG‑CATGTa                                                                                               >  2:691537/1‑100 (MQ=255)
        aTGGCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAg                                                                                              <  2:1554378/99‑1 (MQ=255)
           ccTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGc                                                                                          >  2:1460604/1‑100 (MQ=255)
            cTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGCt                                                                                         <  1:1814714/100‑1 (MQ=255)
             tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGcca                                                                                        <  1:2819683/100‑3 (MQ=255)
             tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCtt                                                                                        >  1:302123/1‑100 (MQ=255)
             tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCCGCT‑CATGTAGCCGCtt                                                                                        >  2:2449807/1‑100 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGCtaat                                                                                      <  1:1246756/100‑4 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGCt                                                                                         <  2:210233/97‑1 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCtttt                                                                                      <  2:378218/100‑1 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCtttt                                                                                      >  1:232104/1‑100 (MQ=255)
               aTAGCGCTGCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGCtaan                                                                                      <  2:2538777/100‑4 (MQ=255)
                  gcgcTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCt                                                                                   >  2:3313392/1‑100 (MQ=255)
                  gcgcTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCt                                                                                   >  1:1277766/1‑100 (MQ=255)
                   cgcttcgcttATCAGGCCTACCCGCTCGACATCGTCATTTGTGGTTCGGCAAAAGCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGCtaattgag                                                                                  <  2:1970096/100‑8 (MQ=37)
                   cgcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATAAGCGGCTACATG‑AGCCGCtaattgag                                                                                  <  2:640360/100‑8 (MQ=37)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGCtaattgagg                                                                                 <  1:1241314/100‑9 (MQ=38)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATG‑AGCCGCtaattgagg                                                                                 <  1:3269484/100‑9 (MQ=38)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTgc                                                                                 <  1:1873422/100‑1 (MQ=255)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTgc                                                                                 <  2:528108/100‑1 (MQ=255)
                     cttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATTCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTgcg                                                                                >  1:3991140/1‑100 (MQ=255)
                      ttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTgcgc                                                                               <  1:530621/100‑1 (MQ=255)
                       tcgcttATCATGCCTAAGCGCTCCAAACCGTCATTTGTGGTTCGGCAAAATCTTCCATAATGCCTCAATTAGCGGCTACATG‑AGCCGCtaattgaggc                                                                                <  2:958850/98‑10 (MQ=16)
                       tcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGca                                                                              >  2:3024900/1‑100 (MQ=255)
                       tcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGca                                                                              <  1:645078/100‑1 (MQ=255)
                        cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGcac                                                                             >  1:965917/1‑100 (MQ=255)
                        cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGcac                                                                             <  2:2333376/100‑1 (MQ=255)
                        cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGcac                                                                             >  2:603841/1‑100 (MQ=255)
                         gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGcaca                                                                            >  1:3688071/1‑100 (MQ=255)
                         gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGcaca                                                                            <  2:1613349/100‑1 (MQ=255)
                         gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGcaca                                                                            >  1:2685689/1‑100 (MQ=255)
                           ttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGcacaca                                                                          <  1:1588074/100‑1 (MQ=255)
                            tATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCacacac                                                                         >  2:57719/1‑100 (MQ=255)
                            tATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCacacac                                                                         >  2:1085155/1‑100 (MQ=255)
                              tCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGc                                                                       >  1:2473974/1‑100 (MQ=255)
                              tCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCACAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGc                                                                       >  2:3958895/1‑100 (MQ=255)
                               cAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGcc                                                                      >  1:966093/1‑100 (MQ=255)
                                 ggCCTACCCGCGCGACAACGTCATTTGTGCTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGccc                                                                     <  2:3137973/99‑1 (MQ=255)
                                   ccTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGc                                                                  <  2:402213/100‑1 (MQ=255)
                                      actcgcgcgACAACGTCATTTGGGGTTCGGCGAAATATTGCCGAATGACTCAGTTAGCGGCT‑CATGTAGCCGCTTTGTCTGCGCACACACGCCCAg                                                                   >  2:2563669/4‑96 (MQ=255)
                                      aCCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCa                                                               >  2:2292715/1‑100 (MQ=255)
                                       cccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCa                                                               >  2:3669643/1‑99 (MQ=255)
                                        ccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCaaa                                                             >  1:3031880/1‑100 (MQ=255)
                                        ccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTACAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACAAACGCCCAGCTCaaa                                                             >  2:1040801/1‑100 (MQ=255)
                                            gcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGc                                                         <  1:224063/100‑1 (MQ=255)
                                             cgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCg                                                        <  1:1086338/100‑1 (MQ=255)
                                             cgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCg                                                        <  1:3159100/100‑1 (MQ=255)
                                                cAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGttt                                                     <  1:3775935/100‑1 (MQ=255)
                                                cAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGttt                                                     <  1:797526/100‑1 (MQ=255)
                                                 aaCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGtttt                                                    >  2:421178/1‑100 (MQ=255)
                                                  aCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtc                                                   >  1:1692436/1‑100 (MQ=255)
                                                  aCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtc                                                   <  1:675974/100‑1 (MQ=255)
                                                    gTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtctc                                                 <  1:2431550/100‑1 (MQ=255)
                                                    gTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCCCTTTTTCTGCGCACACACGCCCAGCTCACACGGCGTTTtcac                                                 <  1:2837867/100‑3 (MQ=255)
                                                     tCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                >  2:2257937/1‑100 (MQ=255)
                                                     tCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                >  1:2372477/1‑100 (MQ=255)
                                                       aTTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATc                                              <  1:1318360/100‑1 (MQ=255)
                                                        tttGTGGTTCGGCAAAATCTTCCAGAACGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCg                                             >  2:1118044/1‑100 (MQ=255)
                                                         ttGTGGTTCGGCAAAATCTTCCAGACTGCCTCAATTAGCCGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGc                                            <  2:3929745/100‑1 (MQ=255)
                                                         ttGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGc                                            <  2:2590676/100‑1 (MQ=255)
                                                           gtgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTg                                          <  1:3091720/100‑1 (MQ=255)
                                                            tgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgc                                         >  2:531025/1‑100 (MQ=255)
                                                                  ggCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTa                                   >  1:3350854/1‑100 (MQ=255)
                                                                   gCAACATCTTCCCGCTTGCCTCAGCTAGCGGCT‑CAGGTAGCCGCTTTTTCTGCGCACGCGCGCCCAGCTTAAACGGCGTTTTCTCATCGCTGCGCTCTa                                   <  2:3810289/99‑1 (MQ=255)
                                                                   gCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTaa                                  >  2:110857/1‑100 (MQ=255)
                                                                     aaaaTCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTaaaa                                <  2:2372472/100‑1 (MQ=255)
                                                                     aaaaTCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTaaaa                                <  1:3879750/100‑1 (MQ=255)
                                                                       aaTCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatc                              <  2:4154778/100‑1 (MQ=255)
                                                                        aTCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatca                             <  2:1374187/100‑1 (MQ=255)
                                                                          cTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatcat                            >  1:1503895/1‑99 (MQ=255)
                                                                            tCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAc                         >  1:3565212/1‑100 (MQ=255)
                                                                              ccgAATGTCTCAATTACCGGCT‑CATGTAGCCGCTTTTTCTCCGCACACACGCCCAGCTCAAACCGCGTTTTCTCATCGCTGCGCTCTAAAATAATCAggc                       >  2:3152178/3‑100 (MQ=255)
                                                                              cAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAcgc                       <  2:3582984/100‑1 (MQ=255)
                                                                               aGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACGCGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAAGCAGCAcgcg                      <  2:820981/100‑1 (MQ=255)
                                                                                  aTGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGtt                   <  1:648763/100‑1 (MQ=255)
                                                                                   gcgcagaaaaaGCGGCTACATG‑AGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTa                  >  2:1736021/11‑100 (MQ=37)
                                                                                     gcagaaaaaGCGGCTACATG‑AGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACg                >  1:2349949/9‑100 (MQ=38)
                                                                                     ccTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACg                <  2:2477914/100‑1 (MQ=255)
                                                                                     ccTCAATTAGCGGCC‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCCGCTCAACCGGCGCTTTCTCATCCCTGCGCTCTAAAATCATCACGCGGTTCCg                >  1:116529/1‑100 (MQ=255)
                                                                                      cTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACGCGCCCAGCTCAAACGGCGTCTTCTCACCGCTGCGCTCTAAAATCATCACGCGGTTACGc               <  2:281237/100‑1 (MQ=255)
                                                                                      cTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGc               <  1:456509/100‑1 (MQ=255)
                                                                                      cTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGc               >  1:110806/1‑100 (MQ=255)
                                                                                        gaaaaaGCGGCTACATG‑AGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAc             >  1:1976100/6‑100 (MQ=255)
                                                                                        gaaaaaGCGGCTACATG‑AGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAc             >  1:1089156/6‑100 (MQ=255)
                                                                                        cAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCa              <  2:362195/99‑1 (MQ=255)
                                                                                         aaTTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAcc            <  2:2802031/100‑1 (MQ=255)
                                                                                          aaaaGCGGCTACATG‑AGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCg           >  1:4134863/4‑100 (MQ=255)
                                                                                          aaaaGCGGCTACATG‑AGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCg           >  2:779704/4‑100 (MQ=255)
                                                                                             aGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGtt        >  1:841352/1‑100 (MQ=255)
                                                                                               cGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCg      <  2:858419/100‑1 (MQ=255)
                                                                                                ggCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgg     >  1:2536570/1‑100 (MQ=255)
                                                                                                 gCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCggg    >  1:3973943/1‑100 (MQ=255)
                                                                                                  nt‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgggg   >  1:229834/2‑100 (MQ=255)
                                                                                                    |                                                                                                   
AGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCT‑CATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGC  >  NC_000913/3956874‑3957072

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: ‑