Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,344,033 A→G 15.0% L105P (CTG→CCG)  yciT ← global regulator of transcription; DeoR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,344,0330AG15.0% 298.5 / 43.2 132L105P (CTG→CCG) yciTglobal regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand):  ref base A (47/65);  new base G (11/9);  total (58/74)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

Only 100 of 139 total aligned reads displayed.

ACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCA  >  NC_000913/1343934‑1344132
                                                                                                   |                                                                                                   
 cGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCag                                                                                                    <  1:3230900/100‑1 (MQ=255)
 cGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCag                                                                                                    <  1:367401/100‑1 (MQ=255)
   ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCagag                                                                                                  <  1:2489156/100‑1 (MQ=255)
   ccACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGATGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCagag                                                                                                  >  2:1165084/1‑100 (MQ=255)
    cACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGt                                                                                                 <  1:803780/100‑1 (MQ=255)
           cAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                          <  2:885282/100‑1 (MQ=255)
           cAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                          <  2:3435123/100‑1 (MQ=255)
           cAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                          <  1:1489290/100‑1 (MQ=255)
            aaaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                          >  2:874209/1‑99 (MQ=255)
            aaaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGa                                                                                         >  2:3403515/1‑100 (MQ=255)
            aaaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGa                                                                                         >  1:1914045/1‑100 (MQ=255)
            aaaaTAACTTCACAAGGCGCGTCTTTCAGCAAATAAGCGATGTAGCTACTGACCGTGCTGATAGTGACATTATTCTTCTGCTGGCAAAGAGTACggg                                                                                            >  2:932269/1‑97 (MQ=255)
              aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAgc                                                                                       <  1:2291133/100‑2 (MQ=255)
              aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGt                                                                                       <  2:1486678/100‑1 (MQ=255)
              aaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGt                                                                                       <  2:777635/100‑1 (MQ=255)
               aTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAgca                                                                                      <  1:2663727/100‑3 (MQ=255)
               aTAACTGCACGTGTCGCGGGTTTCAGCGAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                          <  2:3530666/95‑1 (MQ=255)
                tAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAGCAg                                                                                     <  2:1604343/100‑1 (MQ=255)
                tAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAGCAg                                                                                     <  2:1377756/100‑1 (MQ=255)
                 aaCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCGGACTCTGGGCGAGCAga                                                                                    <  1:1767861/100‑2 (MQ=255)
                 aaCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                    >  2:2271269/1‑100 (MQ=255)
                  aCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGaggagg                                                                                    >  2:3110242/1‑99 (MQ=255)
                  aCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAggg                                                                                   <  1:3344838/100‑1 (MQ=255)
                  aCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAggg                                                                                   <  1:1423487/100‑1 (MQ=255)
                  aCTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAGCAgaa                                                                                   <  1:879639/100‑3 (MQ=37)
                   cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTTTTTTTTTTTCCCAGAGTCGGGGCGAGTAGGGc                                                                                  >  2:875714/1‑100 (MQ=255)
                   cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                  >  1:292886/1‑100 (MQ=255)
                   cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAgcaagat                                                                                  >  1:641683/1‑94 (MQ=38)
                    ttCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCa                                                                                 >  2:3617034/1‑100 (MQ=255)
                    tacacgAGGCTCGTCTTTAAGCAATTGAGCGATGTAGCTGCTGCCCGTGGTGAGAGTGACATTTTTCTTCTGCTCGCCCGGACGCTGGGCGAgcagaag                                                                                  <  2:1669810/97‑7 (MQ=255)
                     tCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCAt                                                                                <  2:3298249/100‑1 (MQ=255)
                      cacaAGGCGCGTCTTTCAGCAAATGCGCGACGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCAtt                                                                               <  1:3410736/100‑1 (MQ=255)
                        caAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgc                                                                             <  1:1345622/100‑1 (MQ=255)
                        caAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgc                                                                             <  1:1596075/100‑1 (MQ=255)
                           ggCGAGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTGCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgc                                                                          <  2:254522/100‑1 (MQ=255)
                             cgcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                        <  1:3806229/100‑1 (MQ=255)
                             cgcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                        >  1:449949/1‑100 (MQ=255)
                              gcgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctgc                                                                       <  1:3597724/100‑1 (MQ=255)
                               cgTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGcc                                                                      <  2:3438280/100‑1 (MQ=255)
                                  cTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtt                                                                   >  1:1121054/1‑100 (MQ=255)
                                    ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtttt                                                                 <  1:3534105/100‑1 (MQ=255)
                                    ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtttt                                                                 >  1:1599469/1‑100 (MQ=255)
                                     tCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTc                                                                <  2:1262324/100‑1 (MQ=255)
                                     tCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCGGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTc                                                                <  2:371202/100‑1 (MQ=255)
                                     tCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGGGACATTTTTCTTCTGCTCGCCCAGAGTCCAGGCGAGTAGGGCATTGCTGCTGCCATTTTc                                                                <  2:976600/100‑1 (MQ=255)
                                      cAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                               <  2:2767341/100‑1 (MQ=255)
                                       aGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGTTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGa                                                              <  1:236130/100‑1 (MQ=255)
                                        gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                             >  1:2902017/1‑100 (MQ=255)
                                        gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                             <  2:3204245/100‑1 (MQ=255)
                                        gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                             >  2:1277035/1‑100 (MQ=255)
                                         ccaattgccgATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                            <  2:75411/92‑1 (MQ=255)
                                         cAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATa                                                            >  1:2244608/1‑100 (MQ=255)
                                          aaaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATaa                                                           <  2:2918202/100‑1 (MQ=255)
                                              gcgcGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGAt                                                       <  1:2474081/100‑1 (MQ=255)
                                               cgcgATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATg                                                      >  1:126003/1‑100 (MQ=255)
                                                gcgATGTAGCTGCTGGCGGTGATGAGAGTGACATTTTTCTTCTGCTCGCCTAGTGTCCGGGCTAGTGGGGCAGGGCTGCTGCCATTTTCTATATAGAttg                                                     <  2:1114560/100‑3 (MQ=255)
                                                    tGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                 <  2:475837/100‑1 (MQ=255)
                                                    tGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTc                                                 >  1:2667524/1‑100 (MQ=255)
                                                     gTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCg                                                <  2:3677428/100‑1 (MQ=255)
                                                     gTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCg                                                <  1:3257598/100‑1 (MQ=255)
                                                       aGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGcc                                              >  1:2614579/1‑100 (MQ=255)
                                                       aGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGcc                                              <  1:970894/100‑1 (MQ=255)
                                                          tgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATCCAGCTGGTTTAGCCCgg                                           >  1:1652179/1‑100 (MQ=255)
                                                           gctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGt                                          >  1:1568058/1‑100 (MQ=255)
                                                           gctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGt                                          <  2:1933341/100‑1 (MQ=255)
                                                            ctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGtt                                         <  1:2206188/100‑1 (MQ=255)
                                                              gACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGa                                       >  1:2920999/1‑100 (MQ=255)
                                                                ccGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAAc                                     >  2:1988709/1‑100 (MQ=255)
                                                                ccGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAAc                                     >  2:815252/1‑100 (MQ=255)
                                                                 cGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAAcc                                    <  2:2105459/100‑1 (MQ=255)
                                                                  gTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTGGCCCGGTTg                                        >  2:2245365/1‑95 (MQ=255)
                                                                   tgatgaTAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAg                                  >  1:3367902/1‑100 (MQ=255)
                                                                    gatgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGc                                 <  1:605730/100‑1 (MQ=255)
                                                                      tgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGa                               >  2:2772961/1‑100 (MQ=255)
                                                                         tAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAcgc                            <  2:1895850/100‑1 (MQ=255)
                                                                           gTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCgg                          <  2:2616723/100‑1 (MQ=255)
                                                                           gTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCgg                          <  2:3005804/100‑1 (MQ=255)
                                                                           gTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCgg                          <  2:3070224/100‑1 (MQ=255)
                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGc                         >  2:1543058/1‑99 (MQ=255)
                                                                             gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCa                        <  2:15449/100‑1 (MQ=255)
                                                                               cATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCaaa                      >  2:483242/1‑100 (MQ=255)
                                                                               cATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCaaa                      >  1:2788210/1‑100 (MQ=255)
                                                                                aTTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAAc                     <  1:2220068/100‑1 (MQ=255)
                                                                                      cttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCg               >  1:1644535/1‑100 (MQ=255)
                                                                                      cttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCg               <  2:3787241/100‑1 (MQ=255)
                                                                                       ttctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGa              <  2:81944/100‑1 (MQ=255)
                                                                                       ccctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGa              >  1:2331212/6‑100 (MQ=37)
                                                                                        tctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCGGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAg             >  2:1512129/1‑100 (MQ=255)
                                                                                        cctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAg             >  1:2898121/5‑100 (MQ=38)
                                                                                         ctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt            >  1:1132670/4‑100 (MQ=255)
                                                                                         ctacTCGCCCGGACTCTGGGCGAGCAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt            >  1:1959111/4‑100 (MQ=38)
                                                                                           gCTCGCCCACAGTCCGTGCTAGTACGGCATTGCTGCTGCCATTTTAGATAAAGATGGTTTCGCCCGGTTGAGCCAGCGACGCGGCAAACTCGGCGAGTTc          <  2:23841/100‑1 (MQ=255)
                                                                                            cTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCa         >  2:8051/1‑100 (MQ=255)
                                                                                             tCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCAc        <  1:1186044/100‑1 (MQ=255)
                                                                                             tCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCAc        <  1:2388769/100‑1 (MQ=255)
                                                                                              cGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACg       >  1:249592/1‑100 (MQ=255)
                                                                                               gCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGc      <  1:1294887/100‑1 (MQ=255)
                                                                                                cccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCt     >  2:1630505/1‑100 (MQ=255)
                                                                                                 ccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTGTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGtt           <  2:1900318/93‑1 (MQ=255)
                                                                                                   agagTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa  <  1:387005/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCA  >  NC_000913/1343934‑1344132

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