Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,934,224 T→A 5.3% L127L (CTA→CTT edd ← 6‑phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,934,2240TA 5.3% 390.6 / 6.6 109L127L (CTA→CTTedd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (55/48);  new base A (3/3);  total (58/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 112 total aligned reads displayed.

GCACAAACACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTGGTGCGCGGAGA  >  NC_000913/1934076‑1934372
                                                                                                                                                    |                                                                                                                                                    
gCACAAACACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCtt                                                                                                                                                      <  4:419427/149‑1 (MQ=255)
  acaaacaaCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACATAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTag                                                                                                                                                    >  4:434251/1‑149 (MQ=255)
  acaAACACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTag                                                                                                                                                    >  2:313828/1‑149 (MQ=255)
  acaAACACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTag                                                                                                                                                    >  1:69144/1‑149 (MQ=255)
  acaAACACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTag                                                                                                                                                    >  3:119710/1‑149 (MQ=255)
  acaAACACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTag                                                                                                                                                    <  4:230585/149‑1 (MQ=255)
    aaaCACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTagca                                                                                                                                                  >  1:273476/1‑149 (MQ=255)
    aaaCACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTagca                                                                                                                                                  >  2:561331/1‑149 (MQ=255)
     aaCACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGACCAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTGCGCGGCTTagcag                                                                                                                                                 >  4:127058/1‑149 (MQ=255)
      acacCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTagcagc                                                                                                                                                <  3:299138/149‑1 (MQ=255)
         ccGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGAc                                                                                                                                             <  3:422363/149‑1 (MQ=255)
          cGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCTGCTAAGCCGCGaag                                                                                                                                            >  1:644162/1‑147 (MQ=255)
          cGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACa                                                                                                                                            <  1:489947/149‑1 (MQ=255)
            cAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCTGCTAAGCCGCGaagtg                                                                                                                                          <  2:644162/149‑5 (MQ=255)
             aGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGATAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAatt                                                                                                                                         <  3:487833/149‑1 (MQ=255)
             aGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAatt                                                                                                                                         <  1:616824/149‑1 (MQ=255)
                 aaaTGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCTGCTAAGCCGCGaagt                                                                                                                                           >  2:380729/1‑140 (MQ=255)
                 aaaTGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCTGCTAAGCCGCGaagt                                                                                                                                           <  1:380729/143‑4 (MQ=255)
                 aaaTGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAattccat                                                                                                                                     >  2:514749/1‑149 (MQ=255)
                  aaTGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAattccatt                                                                                                                                    <  2:483946/149‑1 (MQ=255)
                  aaTGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAattccatt                                                                                                                                    <  4:54964/149‑1 (MQ=255)
                   aTGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAattccattc                                                                                                                                   <  3:292176/149‑1 (MQ=255)
                   aTGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAattccattc                                                                                                                                   >  3:581183/1‑149 (MQ=255)
                         aaaCGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCt                                                                                                                             >  3:234793/1‑149 (MQ=255)
                         aaaCGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCt                                                                                                                             >  2:49787/1‑149 (MQ=255)
                             gACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGccc                                                                                                                         >  4:631639/1‑149 (MQ=255)
                                 gggCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTgg                                                                                                                     <  3:403431/149‑1 (MQ=255)
                                 gggCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTgg                                                                                                                     <  2:606988/149‑1 (MQ=255)
                                 gggCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTgg                                                                                                                     <  3:670476/149‑1 (MQ=255)
                                   gCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGt                                                                                                                   <  1:549086/149‑1 (MQ=255)
                                    cTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTg                                                                                                                  >  3:572513/1‑149 (MQ=255)
                                     tGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGa                                                                                                                 <  3:195145/149‑1 (MQ=255)
                                      gCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGGCCCTGGGTGac                                                                                                                <  1:514749/149‑1 (MQ=255)
                                           cGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCAt                                                                                                           <  1:633261/149‑1 (MQ=255)
                                            gTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATc                                                                                                          <  1:561331/149‑1 (MQ=255)
                                             tCAGACCCGGGTTAGTCTTGTGGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCTGCTTAGCAGCGACAATTCCATTCCATCCTGCGCCTGGGTGACACCATca                                                                                                         >  4:123217/1‑149 (MQ=255)
                                             tCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATca                                                                                                         <  3:119322/149‑1 (MQ=255)
                                               aGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATcaca                                                                                                       >  4:541849/1‑149 (MQ=255)
                                                     gggACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCTGCTAAGCcgcgtagt                                                                                                                                           >  3:450361/1‑103 (MQ=255)
                                                     gggACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCTGCTAAGCcgcgtagt                                                                                                                                           <  4:450361/107‑5 (MQ=255)
                                                       gACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCg                                                                                               >  2:510340/1‑149 (MQ=255)
                                                       gACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCg                                                                                               <  4:383815/149‑1 (MQ=255)
                                                        aCAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCgg                                                                                              >  2:672245/1‑149 (MQ=255)
                                                            tCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAaca                                                                                          >  4:577078/1‑149 (MQ=255)
                                                                         acCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCt                                                                             <  1:49787/149‑1 (MQ=255)
                                                                           cGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGa                                                                           <  2:610670/149‑1 (MQ=255)
                                                                           cGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGa                                                                           >  3:674505/1‑149 (MQ=255)
                                                                             aGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCg                                                                                                  <  1:223682/124‑1 (MQ=255)
                                                                             aGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCg                                                                                                  >  2:223682/1‑124 (MQ=255)
                                                                                aaCAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACAccgccc                                                                                    <  2:30060/135‑1 (MQ=255)
                                                                                aaCAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACAccgccc                                                                                    >  1:30060/1‑135 (MQ=255)
                                                                                 aCAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAAcc                                                                     <  1:170741/149‑1 (MQ=255)
                                                                                   agagCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATca                                                                                                         <  1:306370/111‑1 (MQ=255)
                                                                                   agagCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATca                                                                                                         >  2:306370/1‑111 (MQ=255)
                                                                                    gagCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCa                                                                  >  3:131676/1‑149 (MQ=255)
                                                                                     agCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCAt                                                                 >  1:14618/1‑149 (MQ=255)
                                                                                         ccATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGc                                                             >  4:602039/1‑149 (MQ=255)
                                                                                         ccATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGc                                                             <  2:438597/149‑1 (MQ=255)
                                                                                           aTCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCAtt                                                                >  4:243236/1‑144 (MQ=255)
                                                                                           aTCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCAtt                                                                <  3:243236/144‑1 (MQ=255)
                                                                                            tCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTc                                                          <  4:572513/149‑1 (MQ=255)
                                                                                                  aTGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAg                                                    <  3:631639/149‑1 (MQ=255)
                                                                                                   tGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAgg                                                   >  4:473297/1‑149 (MQ=255)
                                                                                                    gTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAggg                                                  >  1:559993/1‑149 (MQ=255)
                                                                                                     ttATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGc                                                 >  4:423168/1‑149 (MQ=255)
                                                                                                        tGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCttt                                              <  3:602039/149‑1 (MQ=255)
                                                                                                         gggACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTa                                             <  1:368397/149‑1 (MQ=255)
                                                                                                         gggACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTa                                             >  2:210395/1‑149 (MQ=255)
                                                                                                           gACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACg                                           >  4:114678/1‑149 (MQ=255)
                                                                                                           gACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACg                                           >  1:371502/1‑149 (MQ=255)
                                                                                                            aCAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCAt                                                                                                           >  2:107558/1‑84 (MQ=255)
                                                                                                            aCAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCAt                                                                                                           <  1:107558/84‑1 (MQ=255)
                                                                                                             cAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCAt                                                        >  4:444368/1‑134 (MQ=255)
                                                                                                             cAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCAt                                                        <  3:444368/134‑1 (MQ=255)
                                                                                                              aGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAAt                                        <  1:672245/149‑1 (MQ=255)
                                                                                                                ccccACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGa                                      >  1:617759/1‑149 (MQ=255)
                                                                                                                ccccACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGa                                      >  1:662481/1‑149 (MQ=255)
                                                                                                                   cACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGGttt                                   >  1:289956/1‑149 (MQ=255)
                                                                                                                        ccgcAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTTCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGa                              <  3:434251/149‑1 (MQ=255)
                                                                                                                        ccgcAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGa                              <  2:674387/149‑1 (MQ=255)
                                                                                                                         cgcAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGAt                             >  2:60705/1‑149 (MQ=255)
                                                                                                                          gcAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATa                            >  3:584379/1‑149 (MQ=255)
                                                                                                                               cATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGtt                       <  2:617759/149‑1 (MQ=255)
                                                                                                                                aTCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTc                      >  2:415809/1‑149 (MQ=255)
                                                                                                                                aTCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTc                      >  1:324876/1‑149 (MQ=255)
                                                                                                                                aTCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTc                      >  3:590961/1‑149 (MQ=255)
                                                                                                                                aTCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTc                      <  1:514633/149‑1 (MQ=255)
                                                                                                                                   gCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATa                   >  3:159749/1‑149 (MQ=255)
                                                                                                                                   gCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATa                   >  2:295743/1‑149 (MQ=255)
                                                                                                                                    cTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGGCACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTAAGAATGATTTCTGGATAGTGTTCATaa                  <  2:273476/149‑1 (MQ=255)
                                                                                                                                    cTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGACTTCTGGATAGTGTTCAtga                  >  4:567588/1‑147 (MQ=255)
                                                                                                                                     tATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAg                 <  3:423168/149‑1 (MQ=255)
                                                                                                                                     tATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAg                 >  2:538133/1‑149 (MQ=255)
                                                                                                                                         aCTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTg             <  4:301537/149‑1 (MQ=255)
                                                                                                                                         aCTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTg             >  3:56285/1‑149 (MQ=255)
                                                                                                                                         aCTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTg             >  3:344447/1‑149 (MQ=255)
                                                                                                                                               cgGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTGGTgcgc       <  4:581183/149‑1 (MQ=255)
                                                                                                                                                  cTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTGGTGCGCGga    <  3:650879/149‑1 (MQ=255)
                                                                                                                                                  cTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTGGTGCGCGga    <  3:577078/149‑1 (MQ=255)
                                                                                                                                                  cTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTGGTGCGCGga    <  3:567588/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACAAACACCGCAGGCAAATGACCAAACGACAGGGCTGCCATCGTCAGACCCGGGACAATCTTGTCGCACACACCGAGGAACAGAGCACCATCAAACATGTTATGGGACAGCCCCACCGCCGCAGACATCGCTATCACTTCGCGGCTTAGCAGCGACAATTCCATTCCATCCTGCCCCTGGGTGACACCATCACACATCGCCGGAACACCGCCCGCAACCTGACCAACCGCATTCGCTTCATGCAGGGCTTTACGAATGATTTCTGGATAGTGTTCATAAGGCTGGTGCGCGGAGA  >  NC_000913/1934076‑1934372

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