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.

AGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGC  >  NC_000913/3956874‑3957072
                                                                                                    |                                                                                                   
 gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTc                                                                                                    >  1:2069955/1‑100 (MQ=255)
 gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTc                                                                                                    >  2:1204250/1‑100 (MQ=255)
 gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGGTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCt                                                                                                      >  1:3681560/1‑99 (MQ=255)
  ggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCa                                                                                                   >  1:1889445/1‑100 (MQ=255)
   gAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCAt                                                                                                  >  1:2917887/1‑100 (MQ=255)
    aaaaaTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATg                                                                                                 <  2:2173036/100‑1 (MQ=255)
     aaaaTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGt                                                                                                >  1:1402510/1‑100 (MQ=255)
      aaaTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGATCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCGCATGTa                                                                                               >  2:691537/1‑100 (MQ=255)
        aTGGCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAg                                                                                              <  2:1554378/99‑1 (MQ=255)
           ccTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGc                                                                                          >  2:1460604/1‑100 (MQ=255)
            cTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCt                                                                                         <  1:1814714/100‑1 (MQ=255)
             tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGcca                                                                                        <  1:2819683/100‑3 (MQ=255)
             tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCtt                                                                                        >  1:302123/1‑100 (MQ=255)
             tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCCGCTCATGTAGCCGCtt                                                                                        >  2:2449807/1‑100 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                      <  1:1246756/100‑4 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCt                                                                                         <  2:210233/97‑1 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCtttt                                                                                      <  2:378218/100‑1 (MQ=255)
               aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCtttt                                                                                      >  1:232104/1‑100 (MQ=255)
               aTAGCGCTGCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaan                                                                                      <  2:2538777/100‑4 (MQ=255)
                  gcgcTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCt                                                                                   >  2:3313392/1‑100 (MQ=255)
                  gcgcTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCt                                                                                   >  1:1277766/1‑100 (MQ=255)
                   cgcttcgcttATCAGGCCTACCCGCTCGACATCGTCATTTGTGGTTCGGCAAAAGCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaattgag                                                                                  <  2:1970096/100‑8 (MQ=37)
                   cgcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATAAGCGGCTACATGAGCCGCtaattgag                                                                                  <  2:640360/100‑8 (MQ=37)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaattgagg                                                                                 <  1:1241314/100‑9 (MQ=38)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaattgagg                                                                                 <  1:3269484/100‑9 (MQ=38)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgc                                                                                 <  1:1873422/100‑1 (MQ=255)
                    gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgc                                                                                 <  2:528108/100‑1 (MQ=255)
                     cttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATTCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgcg                                                                                >  1:3991140/1‑100 (MQ=255)
                      ttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgcgc                                                                               <  1:530621/100‑1 (MQ=255)
                       tcgcttATCATGCCTAAGCGCTCCAAACCGTCATTTGTGGTTCGGCAAAATCTTCCATAATGCCTCAATTAGCGGCTACATGAGCCGCtaattgaggc                                                                                <  2:958850/98‑10 (MQ=16)
                       tcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGca                                                                              >  2:3024900/1‑100 (MQ=255)
                       tcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGca                                                                              <  1:645078/100‑1 (MQ=255)
                        cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcac                                                                             >  1:965917/1‑100 (MQ=255)
                        cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcac                                                                             <  2:2333376/100‑1 (MQ=255)
                        cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcac                                                                             >  2:603841/1‑100 (MQ=255)
                         gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcaca                                                                            >  1:3688071/1‑100 (MQ=255)
                         gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcaca                                                                            <  2:1613349/100‑1 (MQ=255)
                         gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcaca                                                                            >  1:2685689/1‑100 (MQ=255)
                           ttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcacaca                                                                          <  1:1588074/100‑1 (MQ=255)
                            tATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCacacac                                                                         >  2:57719/1‑100 (MQ=255)
                            tATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCacacac                                                                         >  2:1085155/1‑100 (MQ=255)
                              tCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGc                                                                       >  1:2473974/1‑100 (MQ=255)
                              tCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCACAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGc                                                                       >  2:3958895/1‑100 (MQ=255)
                               cAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGcc                                                                      >  1:966093/1‑100 (MQ=255)
                                 ggCCTACCCGCGCGACAACGTCATTTGTGCTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGccc                                                                     <  2:3137973/99‑1 (MQ=255)
                                   ccTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGc                                                                  <  2:402213/100‑1 (MQ=255)
                                      actcgcgcgACAACGTCATTTGGGGTTCGGCGAAATATTGCCGAATGACTCAGTTAGCGGCTCATGTAGCCGCTTTGTCTGCGCACACACGCCCAg                                                                   >  2:2563669/4‑96 (MQ=255)
                                      aCCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCa                                                               >  2:2292715/1‑100 (MQ=255)
                                       cccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCa                                                               >  2:3669643/1‑99 (MQ=255)
                                        ccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCaaa                                                             >  1:3031880/1‑100 (MQ=255)
                                        ccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTACAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACAAACGCCCAGCTCaaa                                                             >  2:1040801/1‑100 (MQ=255)
                                            gcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGc                                                         <  1:224063/100‑1 (MQ=255)
                                             cgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCg                                                        <  1:1086338/100‑1 (MQ=255)
                                             cgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCg                                                        <  1:3159100/100‑1 (MQ=255)
                                                cAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGttt                                                     <  1:3775935/100‑1 (MQ=255)
                                                cAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGttt                                                     <  1:797526/100‑1 (MQ=255)
                                                 aaCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGtttt                                                    >  2:421178/1‑100 (MQ=255)
                                                  aCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtc                                                   >  1:1692436/1‑100 (MQ=255)
                                                  aCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtc                                                   <  1:675974/100‑1 (MQ=255)
                                                    gTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtctc                                                 <  1:2431550/100‑1 (MQ=255)
                                                    gTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCCCTTTTTCTGCGCACACACGCCCAGCTCACACGGCGTTTtcac                                                 <  1:2837867/100‑3 (MQ=255)
                                                     tCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                >  2:2257937/1‑100 (MQ=255)
                                                     tCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                >  1:2372477/1‑100 (MQ=255)
                                                       aTTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATc                                              <  1:1318360/100‑1 (MQ=255)
                                                        tttGTGGTTCGGCAAAATCTTCCAGAACGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCg                                             >  2:1118044/1‑100 (MQ=255)
                                                         ttGTGGTTCGGCAAAATCTTCCAGACTGCCTCAATTAGCCGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGc                                            <  2:3929745/100‑1 (MQ=255)
                                                         ttGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGc                                            <  2:2590676/100‑1 (MQ=255)
                                                           gtgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTg                                          <  1:3091720/100‑1 (MQ=255)
                                                            tgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgc                                         >  2:531025/1‑100 (MQ=255)
                                                                  ggCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTa                                   >  1:3350854/1‑100 (MQ=255)
                                                                   gCAACATCTTCCCGCTTGCCTCAGCTAGCGGCTCAGGTAGCCGCTTTTTCTGCGCACGCGCGCCCAGCTTAAACGGCGTTTTCTCATCGCTGCGCTCTa                                   <  2:3810289/99‑1 (MQ=255)
                                                                   gCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTaa                                  >  2:110857/1‑100 (MQ=255)
                                                                     aaaaTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTaaaa                                <  2:2372472/100‑1 (MQ=255)
                                                                     aaaaTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTaaaa                                <  1:3879750/100‑1 (MQ=255)
                                                                       aaTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatc                              <  2:4154778/100‑1 (MQ=255)
                                                                        aTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatca                             <  2:1374187/100‑1 (MQ=255)
                                                                          cTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatcat                            >  1:1503895/1‑99 (MQ=255)
                                                                            tCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAc                         >  1:3565212/1‑100 (MQ=255)
                                                                              ccgAATGTCTCAATTACCGGCTCATGTAGCCGCTTTTTCTCCGCACACACGCCCAGCTCAAACCGCGTTTTCTCATCGCTGCGCTCTAAAATAATCAggc                       >  2:3152178/3‑100 (MQ=255)
                                                                              cAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAcgc                       <  2:3582984/100‑1 (MQ=255)
                                                                               aGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACGCGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAAGCAGCAcgcg                      <  2:820981/100‑1 (MQ=255)
                                                                                  aTGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGtt                   <  1:648763/100‑1 (MQ=255)
                                                                                   gcgcagaaaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTa                  >  2:1736021/11‑100 (MQ=37)
                                                                                     gcagaaaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACg                >  1:2349949/9‑100 (MQ=38)
                                                                                     ccTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACg                <  2:2477914/100‑1 (MQ=255)
                                                                                     ccTCAATTAGCGGCCCATGTAGCCGCTTTTTCTGCGCACACACGCCCCGCTCAACCGGCGCTTTCTCATCCCTGCGCTCTAAAATCATCACGCGGTTCCg                >  1:116529/1‑100 (MQ=255)
                                                                                      cTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACGCGCCCAGCTCAAACGGCGTCTTCTCACCGCTGCGCTCTAAAATCATCACGCGGTTACGc               <  2:281237/100‑1 (MQ=255)
                                                                                      cTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGc               <  1:456509/100‑1 (MQ=255)
                                                                                      cTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGc               >  1:110806/1‑100 (MQ=255)
                                                                                        gaaaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAc             >  1:1976100/6‑100 (MQ=255)
                                                                                        gaaaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAc             >  1:1089156/6‑100 (MQ=255)
                                                                                        cAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCa              <  2:362195/99‑1 (MQ=255)
                                                                                         aaTTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAcc            <  2:2802031/100‑1 (MQ=255)
                                                                                          aaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCg           >  1:4134863/4‑100 (MQ=255)
                                                                                          aaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCg           >  2:779704/4‑100 (MQ=255)
                                                                                             aGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGtt        >  1:841352/1‑100 (MQ=255)
                                                                                               cGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCg      <  2:858419/100‑1 (MQ=255)
                                                                                                ggCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgg     >  1:2536570/1‑100 (MQ=255)
                                                                                                 gCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCggg    >  1:3973943/1‑100 (MQ=255)
                                                                                                  ntCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgggg   >  1:229834/2‑100 (MQ=255)
                                                                                                    |                                                                                                   
AGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGC  >  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: