Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,279 C→T 5.8% R46R (CGG→CGA pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,2790CT 5.8% 391.1 / 10.3 121R46R (CGG→CGApdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (50/63);  new base T (5/2);  total (55/66)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 123 total aligned reads displayed.

CCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAG  >  NC_000913/53135‑53427
                                                                                                                                                |                                                                                                                                                    
cccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTCGCGCGACAGGAATTAGCGATAATGTGCCCGTAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGTGTGAGCGGCAAACCGAGCATCGCTGCCCGGt                                                                                                                                                  >  3:578793/1‑149 (MQ=255)
 ccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                                 >  4:167508/1‑149 (MQ=255)
 ccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                                 >  1:151277/1‑149 (MQ=255)
 ccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                                 >  3:102242/1‑149 (MQ=255)
 ccATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                                 >  2:622306/1‑149 (MQ=255)
  cATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGATAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTg                                                                                                                                                >  3:283857/1‑149 (MQ=255)
   aTTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTgg                                                                                                                                               >  2:319925/1‑149 (MQ=255)
    ttttCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGt                                                                                                                                              <  1:142967/149‑1 (MQ=255)
    ttttCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGt                                                                                                                                              >  3:456080/1‑149 (MQ=255)
      ttCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTaa                                                                                                                                            <  3:51264/149‑1 (MQ=255)
      ttCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTaa                                                                                                                                            <  1:322137/149‑1 (MQ=255)
       tCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAg                                                                                                                                           <  2:207761/149‑1 (MQ=255)
       tCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAg                                                                                                                                           <  1:327478/149‑1 (MQ=255)
         aaCCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGga                                                                                                                                         >  2:597819/1‑149 (MQ=255)
         aaCCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGga                                                                                                                                         >  2:329690/1‑149 (MQ=255)
           ccGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGgaga                                                                                                                                       >  2:464653/1‑149 (MQ=255)
            cGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGagag                                                                                                                                      >  2:549141/1‑149 (MQ=255)
             gCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGt                                                                                                                                     <  1:230175/149‑1 (MQ=255)
              cTAGCTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTg                                                                                                                                    <  2:375179/149‑1 (MQ=255)
              cTAACTGCCCCGCAGTGACAGGTGCATGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTACAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCGAACCTAGCATCGCTGCCAGATTGGTAAGGAGAGTg                                                                                                                                    <  4:578793/149‑1 (MQ=255)
              cTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTg                                                                                                                                    <  4:691044/149‑1 (MQ=255)
                   tGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATc                                                                                                                               >  4:129109/1‑149 (MQ=255)
                    gCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCg                                                                                                                              <  3:129109/149‑1 (MQ=255)
                     ccccGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCgg                                                                                                                             >  1:416810/1‑149 (MQ=255)
                            gTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCacaaac                                                                                                                      <  4:310373/149‑1 (MQ=255)
                              gACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCacaaacaa                                                                                                                    <  4:167431/149‑1 (MQ=255)
                               aCAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAaacaac                                                                                                                   <  2:565386/149‑1 (MQ=255)
                               aCAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAaacaac                                                                                                                   <  1:306034/149‑1 (MQ=255)
                                cAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAAcc                                                                                                                  >  1:587138/1‑149 (MQ=255)
                                 aGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCa                                                                                                                 <  1:549141/149‑1 (MQ=255)
                                     gCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTc                                                                                                             <  3:462120/149‑1 (MQ=255)
                                       aCGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGtt                                                                                                              >  3:271272/1‑146 (MQ=255)
                                       aCGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGtt                                                                                                              <  4:271272/146‑1 (MQ=255)
                                       aCGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGa                                                                                                           <  1:130143/149‑1 (MQ=255)
                                        cGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGAc                                                                                                          <  2:180955/149‑1 (MQ=255)
                                            gcgcGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGc                                                                                                      <  3:429707/149‑1 (MQ=255)
                                              gcgACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCa                                                                                                    >  4:433339/1‑149 (MQ=255)
                                                  cAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCTGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGTATCGCTGCCCGGTTGGTACGGAGAGTTGCATCGGCACAAACAACCAGTTCGACCGGCCACTCa                                                                                                <  3:582161/149‑1 (MQ=255)
                                                   aGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCAc                                                                                               >  4:472781/1‑149 (MQ=255)
                                                    ggAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACg                                                                                              >  2:545294/1‑149 (MQ=255)
                                                    ggAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACg                                                                                              >  4:316610/1‑149 (MQ=255)
                                                        gTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgt                                                                                          <  4:456080/149‑1 (MQ=255)
                                                        gTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgt                                                                                          <  3:316610/149‑1 (MQ=255)
                                                         tAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCtgtg                                                                                         <  3:433339/149‑1 (MQ=255)
                                                          aGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGc                                                                                        >  2:548131/1‑149 (MQ=255)
                                                           gCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                                 <  4:1961/91‑1 (MQ=255)
                                                           gCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGtt                                                                                                                                                 >  3:1961/1‑91 (MQ=255)
                                                             gTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAg                                                                                     >  2:658480/1‑149 (MQ=255)
                                                              ttAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCTCTCACGCTGTGCAAGc                                                                                    >  3:601061/1‑149 (MQ=255)
                                                               tAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCt                                                                                   <  2:151277/149‑1 (MQ=255)
                                                                 aTGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTgg                                                                                 >  3:287675/1‑149 (MQ=255)
                                                                  tgtgCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGa                                                                                <  4:257187/149‑1 (MQ=255)
                                                                     gCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACaa                                                                             <  3:167508/149‑1 (MQ=255)
                                                                     gCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACaa                                                                             <  2:178427/149‑1 (MQ=255)
                                                                      cccGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAAc                                                                            >  1:593695/1‑149 (MQ=255)
                                                                       ccGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACt                                                                           >  3:168157/1‑149 (MQ=255)
                                                                       ccGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACt                                                                           >  2:598781/1‑149 (MQ=255)
                                                                         gCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTaa                                                                         <  1:464653/149‑1 (MQ=255)
                                                                               tGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCgg                                                                   <  1:319925/149‑1 (MQ=255)
                                                                                gCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCggg                                                                  <  2:632773/149‑1 (MQ=255)
                                                                                   gTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCa                                                               >  2:218181/1‑149 (MQ=255)
                                                                                    ttGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCaa                                                              <  2:406356/149‑1 (MQ=255)
                                                                                       tgCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATcc                                                           >  3:303305/1‑149 (MQ=255)
                                                                                          aGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCa                                                                                                                 >  2:10587/1‑92 (MQ=255)
                                                                                          aGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCa                                                                                                                 <  1:10587/92‑1 (MQ=255)
                                                                                          aGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcgg                                                        >  3:10179/1‑149 (MQ=255)
                                                                                          aGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCcgg                                                        >  1:536195/1‑149 (MQ=255)
                                                                                                 ttGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTc                                                 >  1:348287/1‑149 (MQ=255)
                                                                                                     gggAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                             <  4:555570/149‑1 (MQ=255)
                                                                                                        aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGa                                          <  3:199731/149‑1 (MQ=255)
                                                                                                        aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGa                                          <  1:658480/149‑1 (MQ=255)
                                                                                                        aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGa                                          <  4:102242/149‑1 (MQ=255)
                                                                                                        aaTAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGa                                          <  4:344958/149‑1 (MQ=255)
                                                                                                            aGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGa                                      <  1:430428/149‑1 (MQ=255)
                                                                                                             gggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGAt                                     <  1:134017/149‑1 (MQ=255)
                                                                                                             gggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGAt                                     >  2:118379/1‑149 (MQ=255)
                                                                                                             gggCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCACGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGAt                                     <  3:472781/149‑1 (MQ=255)
                                                                                                               gCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATca                                   <  2:416810/149‑1 (MQ=255)
                                                                                                                  gAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaa                                <  1:548131/149‑1 (MQ=255)
                                                                                                                   aGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaac                               >  1:400667/1‑149 (MQ=255)
                                                                                                                         aGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTg                         <  3:95344/149‑1 (MQ=255)
                                                                                                                         aGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTg                         <  2:270512/149‑1 (MQ=255)
                                                                                                                           cGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                       <  3:557400/149‑1 (MQ=255)
                                                                                                                             gCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtt                     <  1:218181/149‑1 (MQ=255)
                                                                                                                               caaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTg                         <  3:15692/142‑1 (MQ=255)
                                                                                                                               caaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt                   >  2:195470/2‑149 (MQ=255)
                                                                                                                               caaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt                   >  4:208140/2‑149 (MQ=255)
                                                                                                                               aaaCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGtttt                   >  1:374158/1‑149 (MQ=255)
                                                                                                                                aaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                  >  1:259810/1‑149 (MQ=255)
                                                                                                                                aaCCGGGCAGCGATGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                  >  1:35397/1‑149 (MQ=255)
                                                                                                                                 aCCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                 <  3:469421/149‑1 (MQ=255)
                                                                                                                                  ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc                <  2:168541/149‑1 (MQ=255)
                                                                                                                                  ccGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGCTCGACCGGCCACTCACGCTGTGCAAGCTTGACAACTAAGTCCGTGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTGTTAAc                <  4:283857/149‑1 (MQ=255)
                                                                                                                                   cGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGt                      >  1:275098/1‑142 (MQ=255)
                                                                                                                                   cGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATAACAACACGTTGGGt                      <  2:275098/142‑1 (MQ=255)
                                                                                                                                     aGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCGCGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAt             >  4:399412/1‑149 (MQ=255)
                                                                                                                                       cagcgatGCTCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTa           <  1:195470/143‑1 (MQ=255)
                                                                                                                                            cTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGc      <  4:10179/149‑1 (MQ=255)
                                                                                                                                                ccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAg  >  3:683715/1‑149 (MQ=255)
                                                                                                                                                ccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAg  >  3:384418/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCCATTTTCAACCGCTAACTGCCCCGCAGTGACAGGTGCACGTAGCGCGACAGGAAGTAGCGTTAATGTGCCCGCAGTTTGCGGTTGTGCAGGGGAGTTGGGGGAATAAGGGCGGAGGGTGAGCGGCAAACCGAGCATCGCTGCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAG  >  NC_000913/53135‑53427

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: