Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,151,281 T→A 5.7% V204V (GTT→GTA) ‡ fabG → 3‑oxoacyl‑[acyl‑carrier‑protein] reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,151,2810TA 5.7% 303.3 / 7.7 87V204V (GTT→GTA) ‡fabG3‑oxoacyl‑[acyl‑carrier‑protein] reductase
Reads supporting (aligned to +/- strand):  ref base T (42/40);  new base A (3/2);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

AGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAA  >  NC_000913/1151134‑1151421
                                                                                                                                                   |                                                                                                                                            
aGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTc                                                                                                                                             <  4:629931/149‑1 (MQ=255)
     gCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGAGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCg                                                                                                                                        >  1:775562/1‑149 (MQ=255)
       ttGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCggg                                                                                                                                      <  3:523978/149‑1 (MQ=255)
         gATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTc                                                                                                                                    <  1:735751/149‑1 (MQ=255)
           tCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGc                                                                                                                                  >  4:69696/1‑149 (MQ=255)
           tCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGc                                                                                                                                  >  2:217181/1‑149 (MQ=255)
            cGGCTTCAGTAAATCACTGGCGAGCGAAGTTGCGTCACGCGGTATTACTGTAAACTTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGATGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGcc                                                                                                                                 >  4:21966/1‑149 (MQ=255)
              gCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTc                                                                                                                               <  1:503062/149‑1 (MQ=255)
                       aaTCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGca                                                                                                                      >  4:740626/1‑149 (MQ=255)
                         tCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGcaca                                                                                                                    <  4:247494/149‑1 (MQ=255)
                         tCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGcaca                                                                                                                    >  1:739648/1‑149 (MQ=255)
                           aCTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCTATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAgg                                                                                                                  >  1:767040/1‑149 (MQ=255)
                                cgcgcgAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATc                                                                                                             <  1:577041/149‑1 (MQ=255)
                                  cgcgAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGc                                                                                                           <  4:200023/149‑1 (MQ=255)
                                     gAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTcggc                                                                                                                            <  1:335444/129‑1 (MQ=255)
                                     gAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTcggc                                                                                                                            >  2:335444/1‑129 (MQ=255)
                                     gAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCaa                                                                                                        >  1:235695/1‑149 (MQ=255)
                                            cGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGtt                                                                                                 <  4:538639/149‑1 (MQ=255)
                                                aCGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCAt                                                                                             >  3:114729/1‑149 (MQ=255)
                                                  gcgGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTc                                                                                           <  4:487464/149‑1 (MQ=255)
                                                   cgGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTcc                                                                                          <  1:339446/149‑1 (MQ=255)
                                                   cgGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTcc                                                                                          <  3:224389/149‑1 (MQ=255)
                                                    ggTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCt                                                                                         >  2:395802/1‑149 (MQ=255)
                                                    ggTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCt                                                                                         <  3:288380/149‑1 (MQ=255)
                                                     gTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTg                                                                                        <  3:69696/149‑1 (MQ=255)
                                                     gTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTg                                                                                        <  4:704766/149‑1 (MQ=255)
                                                       aTTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGc                                                                                      >  1:432944/1‑149 (MQ=255)
                                                           cTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATcc                                                                                  >  1:344643/1‑149 (MQ=255)
                                                            tGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGCAATCGCCAACGCGGTTGCATTCCTGGCATCcg                                                                                 >  4:453818/1‑149 (MQ=255)
                                                            tGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCcg                                                                                 >  3:205074/1‑149 (MQ=255)
                                                            tGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCcg                                                                                 >  4:64055/1‑149 (MQ=255)
                                                            tGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTCCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCcg                                                                                 >  4:588686/1‑149 (MQ=255)
                                                             gTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCcga                                                                                >  1:643498/1‑149 (MQ=255)
                                                              tAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCcgac                                                                               >  3:181751/1‑149 (MQ=255)
                                                                aaCGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCg                                                                                                                                        >  2:336971/1‑90 (MQ=255)
                                                                aaCGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCg                                                                                                                                        <  1:336971/90‑1 (MQ=255)
                                                                aaCGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCcgacga                                                                             >  2:633386/1‑149 (MQ=255)
                                                                aaCGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCcgacga                                                                             <  3:100098/149‑1 (MQ=255)
                                                                 aCGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGaa                                                                            <  1:654756/149‑1 (MQ=255)
                                                                  cGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAag                                                                           <  2:635725/149‑1 (MQ=255)
                                                                    ttgttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAagca                                                                         <  3:39459/149‑1 (MQ=255)
                                                                      gttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAagcagc                                                                       <  4:159529/149‑1 (MQ=255)
                                                                      gttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAagcagc                                                                       <  4:198883/149‑1 (MQ=255)
                                                                      gttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAagcagc                                                                       <  3:740626/149‑1 (MQ=255)
                                                                      gttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAagcagc                                                                       >  2:616054/1‑149 (MQ=255)
                                                                      gttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGAACCTGCgcca                                                                                                                                     <  2:338695/87‑4 (MQ=255)
                                                                      gttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGAACCTGCgcca                                                                                                                                     >  1:338695/1‑84 (MQ=255)
                                                                      gttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGCATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAagcagc                                                                       <  3:32586/149‑1 (MQ=255)
                                                                       ttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTc                                                                                                                                             <  4:668292/78‑1 (MQ=255)
                                                                       ttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTc                                                                                                                                             >  3:668292/1‑78 (MQ=255)
                                                                       ttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCt                                                                      >  4:120938/1‑149 (MQ=255)
                                                                         gCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTa                                                                    >  2:354677/1‑149 (MQ=255)
                                                                         gCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTa                                                                    <  2:109770/149‑1 (MQ=255)
                                                                           tCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACa                                                                  >  4:480569/1‑149 (MQ=255)
                                                                                 cTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACgg                                                            >  3:26635/1‑149 (MQ=255)
                                                                                 cTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACgg                                                            <  2:328403/149‑1 (MQ=255)
                                                                                      ttGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGaa                                                       >  1:142619/1‑149 (MQ=255)
                                                                                        gAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCAc                                                              >  1:314276/1‑140 (MQ=255)
                                                                                        gAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCAc                                                              <  2:314276/140‑1 (MQ=255)
                                                                                        gAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAAc                                                     <  2:261519/149‑1 (MQ=255)
                                                                                              gACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCa                                               <  2:592191/149‑1 (MQ=255)
                                                                                               aCATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCAt                                              >  4:57482/1‑149 (MQ=255)
                                                                                                 aTGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCAtgt                                            <  2:236709/149‑1 (MQ=255)
                                                                                                 aTGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCAtgt                                            <  4:181751/149‑1 (MQ=255)
                                                                                                  tGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCAtgtg                                           >  4:30285/1‑149 (MQ=255)
                                                                                                   gACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGa                                          >  2:712545/1‑149 (MQ=255)
                                                                                                      acGTGCGCTGAGCGATGACCAGCGTGAGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAAcg                                       >  4:223457/1‑149 (MQ=255)
                                                                                                           cgcTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCggg                                  <  2:648984/149‑1 (MQ=255)
                                                                                                                aGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTa                             >  3:164003/1‑149 (MQ=255)
                                                                                                                  cGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACa                           >  1:449847/1‑149 (MQ=255)
                                                                                                                   gATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAAcggcgg                                   >  4:321014/1‑140 (MQ=255)
                                                                                                                   gATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAAcggcgg                                   <  3:321014/140‑1 (MQ=255)
                                                                                                                   gATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACAt                          <  3:120938/149‑1 (MQ=255)
                                                                                                                   gATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACAt                          <  3:598205/149‑1 (MQ=255)
                                                                                                                     tGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATgg                        <  2:643498/149‑1 (MQ=255)
                                                                                                                        ccAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCt                     <  2:767040/149‑1 (MQ=255)
                                                                                                                        ccAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCt                     >  4:584541/1‑149 (MQ=255)
                                                                                                                         cAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAAcggcgg                                   <  4:291006/134‑1 (MQ=255)
                                                                                                                         cAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAAcggcgg                                   >  3:291006/1‑134 (MQ=255)
                                                                                                                                   gTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGAttt          >  4:598126/1‑149 (MQ=255)
                                                                                                                                     aTCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGc        <  3:292854/149‑1 (MQ=255)
                                                                                                                                       ccTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGcac      >  1:609152/1‑149 (MQ=255)
                                                                                                                                       ccTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGcac      >  1:635718/1‑149 (MQ=255)
                                                                                                                                         acccgcAGGAACCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCCTTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACaa    >  4:673056/4‑149 (MQ=255)
                                                                                                                                         acccgcAGGAACCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCa                                               >  2:764143/4‑106 (MQ=255)
                                                                                                                                         acccgcAGGAACCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCa                                               <  1:764143/103‑1 (MQ=255)
                                                                                                                                           gcgcAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACaaaa  >  2:40139/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
AGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACATGGTCTGACCGCGATTTGCACAAAA  >  NC_000913/1151134‑1151421

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: