Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,912,074 G→A S532F (TCC→TTC)  glmS ← L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,912,0740GA93.3% 223.8 / 4.7 76S532F (TCC→TTC) glmSL‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (42/28);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

ATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTGAGAC  >  NC_000913/3911834‑3912322
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
aTTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGttt                                                                                                                                                                                                                                                <  1:160604/251‑1 (MQ=255)
      gacccaaccttcaCCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACAGGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGGTTTTCCa                                                                                                                                                                                                                                           <  2:771892/238‑1 (MQ=255)
         tCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCa                                                                                                                                                                                                                                        >  1:77271/1‑250 (MQ=255)
             ccGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCg                                                                                                                                                                                                                                   >  1:673888/1‑251 (MQ=255)
                     gATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTAGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCgg                                                                                                                                                                                                                           >  2:183059/1‑251 (MQ=255)
                       ttttGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTg                                                                                                                                                                                                                         >  1:403268/1‑251 (MQ=255)
                        tttGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCCTCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAAATGACCGCCACGCGCGCGCACTTCTTCAATGTGGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGc                                                                                                                                                                                                                        >  2:707266/1‑251 (MQ=255)
                           gCCAGGTTAAGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGCAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaac                                                                                                                                                                                                                     >  2:60213/1‑251 (MQ=255)
                             cAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaaca                                                                                                                                                                                                                    >  2:173883/1‑250 (MQ=255)
                                   aCGCGTCTGGTAAAAGTCGGTGCCTTTGATCAGCGAGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACATCTTCCACATACGGCATCTCGATGATGTGCATGTTATCGCTACTTACACAACCACCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGGGCGACCTTCTTCACTGTTGAATTTCAGTTTATCCAGCAATTCGTTGTACGGGGCAACAACAATa                                                                                                                                                                                                              >  2:792902/1‑250 (MQ=255)
                                   aCGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATa                                                                                                                                                                                                              <  1:173883/250‑1 (MQ=255)
                                           ggTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCGGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCata                                                                                                                                                                                                     <  1:280847/251‑1 (MQ=255)
                                                   tCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACATCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATc                                                                                                                                                                                             >  2:517559/1‑251 (MQ=255)
                                                        gCCTTTGATAAGCGCGAAATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCTCCTCTTCCACATGCGGCATCTCGATGAGGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTa                                                                                                                                                                                        <  2:370500/251‑1 (MQ=255)
                                                                       gACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCgg                                                                                                                                                                         >  2:275909/1‑250 (MQ=255)
                                                                                aGCCAGCAGCTGCAGCGGCAAGCTGTAGCAGATCGGTGCAATCACCTCTTCCACATGCGGCATATCGATGATGTGCATGATATAGCTATTTACAAAAACCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTTTTCGGTGCAACAAAAATAACCGGCATATCTGCATCACTTAGCGCCAGCGGACCGTGTTTAAGtt                                                                                                                                                                >  2:597181/1‑251 (MQ=255)
                                                                                     gcagcTGCAGCGGAACGGTGTAGAAGACCGGTGCAATCACCTCTTCCACATGGGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCca                                                                                                                                                           <  2:463216/251‑1 (MQ=255)
                                                                                         cTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTAGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCcagca                                                                                                                                                        >  2:627910/1‑250 (MQ=255)
                                                                                             aGCGGAACGGTGTAGAAGATCGGTGAAATTACACCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAGCCCGCATCCTGATCGTCGAAGACATACAAATGACCGCCACGCGCGAGCACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTAGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATACATTAGAGCCAGCGGACCGTGTTTCCGTTCGCCTGCAGCGtc                                                                                                                                                   >  2:508400/1‑250 (MQ=255)
                                                                                                 gAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCt                                                                                                                                               >  2:755389/1‑251 (MQ=255)
                                                                                                 gAACGGTGTAGAAGAACGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGGGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGATCTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATGTCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCt                                                                                                                                               <  1:822309/251‑1 (MQ=255)
                                                                                                    gggTGTAGAAGATCGGTTCAATCACCCTGTCCCCATGTGGCATCCCTATGATGTGCATGTTATCGCTTCTTACAAAACCCGAATCCTGTTCGGCGAAGACATACAACGGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCa                                                                                                                                            <  2:610796/250‑1 (MQ=255)
                                                                                                       tgtAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCg                                                                                                                                         >  1:844281/1‑251 (MQ=255)
                                                                                                        gtAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCg                                                                                                                                         <  1:436558/250‑1 (MQ=255)
                                                                                                        gtAGAAGATCGGGGCAATCACCTCTTCCACATGCGGCATCTCGATGAGGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCg                                                                                                                                         <  2:404260/250‑1 (MQ=255)
                                                                                                        gtAGAAGAGCGGTGCAATCACCTCTTCCACAGTCGGCATCTCGATGATGTGCATGTAATCGCTCCTTCCAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCg                                                                                                                                         <  2:630004/250‑1 (MQ=255)
                                                                                                          agaagaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTg                                                                                                                                       <  1:626596/250‑1 (MQ=255)
                                                                                                             agaTCGGTGCAATCACCTCTCCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATg                                                                                                                                   <  2:620695/251‑1 (MQ=255)
                                                                                                               aTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGt                                                                                                                                  >  1:527075/1‑250 (MQ=255)
                                                                                                                       aaTCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGAtct                                                                                                                         >  1:692886/1‑251 (MQ=255)
                                                                                                                          cACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCtt                                                                                                                      <  1:802173/251‑1 (MQ=255)
                                                                                                                                   ccACATGCGGCCTCCCGAGGATGTGCATGTTATCGCTACTTACAAAACCCGCAGCCTGATCGGCGAAGATATTCAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATACGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGGTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTc                                                                                                              <  1:294341/251‑1 (MQ=255)
                                                                                                                                      cctGCGGCATCTCGATGATGGTCATGTTATCTCTACTTACAAAAACCGCATCATGATCGGCGAATACATACATCTGACCGCCACGCGCTCTAACTTCGTCAATGTTGAATTTCAGTTTTGCCAGCAATACTTTGGTCGGTGCAACAACAATAACCGGCATATCGGCAGCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTGCAGCGTGAATGTAAGAGATCTCTTTCAACTTCaaa                                                                                                           <  2:40131/249‑2 (MQ=255)
                                                                                                                                       aTGCGGCATCTCGCTGATGTGCATGTTATCGATACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAACGCCAGCGGACCGAGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAATAGATCTATTTCAACTTAAATg                                                                                                          >  2:631416/1‑251 (MQ=255)
                                                                                                                                        ggcgcactttcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACACCTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATg                                                                                                          <  2:485518/241‑1 (MQ=255)
                                                                                                                                            gCATCTAGATGATGAGCATGTTATCGCCACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGCGATCTATTTCAACTTCAATgcgctt                                                                                                     >  2:801177/1‑249 (MQ=255)
                                                                                                                                                 ccGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAg                                                                                                <  1:275909/250‑1 (MQ=255)
                                                                                                                                                   gatgatGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAgcg                                                                                              >  1:553236/1‑251 (MQ=255)
                                                                                                                                                         gtgCATGTTATCGCTACTTACCACACCCGCATCCTGATCGGCGCAGACATACACCTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCATGAATGTAAGAGATCTATTTCAACTTCAATGCGCCTTCCAGCGCGATTg                                                                                        >  2:488475/1‑251 (MQ=255)
                                                                                                                                                           gCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTggg                                                                                      >  1:673367/1‑251 (MQ=255)
                                                                                                                                                                    cGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATc                                                                              <  1:517559/250‑1 (MQ=255)
                                                                                                                                                                     gCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCCCCTTCCAGCGCGATTGGGTACTGATCGc                                                                            >  2:378628/1‑251 (MQ=255)
                                                                                                                                                                      cTACTTACAAAACCCGCATCCTGATCGGCGACGACATACCACTGACCGCCACGCGCGCGAACTTCTTCACTGCTGCATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACACCAATAACCGGCATATCGGCATCAATTAGCGTCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCCACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGcc                                                                           >  1:831761/1‑251 (MQ=255)
                                                                                                                                                                           tACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCCGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATACCCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCCATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGc                                                                      >  1:65343/1‑251 (MQ=255)
                                                                                                                                                                              aaaaCCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGcccc                                                                   >  2:767999/1‑250 (MQ=255)
                                                                                                                                                                               aaaCCCGAATCCTGATCGGAGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGACTGTAAGCGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGccc                                                                    >  2:506421/1‑249 (MQ=255)
                                                                                                                                                                                  cccGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGa                                                                >  1:873891/1‑250 (MQ=255)
                                                                                                                                                                                    cGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAAc                                                              >  1:3098/1‑250 (MQ=255)
                                                                                                                                                                                        tCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAgcgc                                                         >  1:41653/1‑251 (MQ=255)
                                                                                                                                                                                         ccTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGcgcg                                                        >  1:786663/1‑251 (MQ=255)
                                                                                                                                                                                           tGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCAGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGACCAGCGCGTg                                                      >  1:226110/1‑251 (MQ=255)
                                                                                                                                                                                            gATCGGCGAAGAAATACAACTTACCGCCACGCGCGCGAATTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGAGTGGGGACTGATCGCCACGGCCCAGGAACAGCGCGTg                                                      <  2:795777/250‑1 (MQ=255)
                                                                                                                                                                                                ggCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGTCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCACTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGtt                                                 >  1:98285/1‑251 (MQ=255)
                                                                                                                                                                                                    aaGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATATCCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTGCCAGCGCGATTGGGGACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTc                                             <  1:506421/251‑1 (MQ=255)
                                                                                                                                                                                                         aTACACCTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAgaga                                        <  2:77271/251‑1 (MQ=255)
                                                                                                                                                                                                              aCTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACAGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGCTCGCCACGGCCCAGGAACAGCGAGTGATGTTTGTCAGGGAAAtct                                   >  1:550629/1‑251 (MQ=255)
                                                                                                                                                                                                               cTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAAtctt                                  >  1:530619/1‑251 (MQ=255)
                                                                                                                                                                                                                tGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAAtcttc                                 >  1:89059/1‑251 (MQ=255)
                                                                                                                                                                                                                cgtacggcactggcgcgAACTTATTCAAAGTTAAATTGCATTTTTGCCAGCAATTCGTTGTCCGGTGCAACAACAATAAGCGGCATAACGGCATCAATTAGCGCCAGCGGTCCGTGTTTCAGTGCGCCAGCAGCGTAGGCTTCAGCGGGAATGGAAGAGATCTCTTTCAACGGCAATGCTCCTTCCAGCGTGATTGGGTACTGAACGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAAtctta                                 <  2:241473/239‑2 (MQ=255)
                                                                                                                                                                                                                   ccGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGc                              <  1:105178/251‑1 (MQ=255)
                                                                                                                                                                                                                       cACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAgc                          >  1:115874/1‑251 (MQ=255)
                                                                                                                                                                                                                       cACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCAGCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAgc                          <  1:253927/251‑1 (MQ=255)
                                                                                                                                                                                                                          gcgcgcgAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCt                       <  1:183059/251‑1 (MQ=255)
                                                                                                                                                                                                                            gcgcgAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTc                     <  1:304918/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgcgAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTc                     >  1:724428/1‑250 (MQ=255)
                                                                                                                                                                                                                               cgAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGAAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCCGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGACTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCCGCGCTTCact                  >  2:365883/1‑249 (MQ=255)
                                                                                                                                                                                                                                gAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGAGCTTCAATg                 >  2:731489/1‑251 (MQ=255)
                                                                                                                                                                                                                                gAACTTCTTAAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGGGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACAGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATg                 >  2:144460/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tcttcAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGtt             <  1:91437/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cttcAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTg           >  1:443642/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ttcAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGt          >  1:526445/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ttcAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGAGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGt          >  2:151508/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tcAAGGTTGGATTTCAGGTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGTCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCACTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGttt            <  2:98285/248‑1 (MQ=255)
                                                                                                                                                                                                                                           aTGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTg      <  1:707266/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGATTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAATTCTTCTGCCAGCGCTTCAATGCGTTTGTCCTgg     >  2:442826/1‑250 (MQ=255)
                                                                                                                                                                                                                                             gTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGCTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGGGTTTCAGTTCGTCAGCAGCGTAGGCTTCATCGTGAATGTAAGAGATCTATTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTCCTGATCGCCACGGACCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTgag    <  1:485055/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ggATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTGAGAc  <  2:443642/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
ATTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCCTGAGAC  >  NC_000913/3911834‑3912322

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