Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 74,027 G→A I165I (ATC→ATT)  thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091374,0270GA97.7% 273.2 / NA 86I165I (ATC→ATT) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (29/55);  total (31/55)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCAT‑CCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCC  >  NC_000913/73782‑74261
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
gcCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGtt                                                                                                                                                                                                                                       >  1:691432/1‑251 (MQ=255)
        cgCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGt                                                                                                                                                                                                                                <  1:666347/249‑1 (MQ=255)
                gggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGt                                                                                                                                                                                                                       <  1:96739/251‑1 (MQ=255)
                 ggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGta                                                                                                                                                                                                                      <  1:80427/251‑1 (MQ=255)
                   cAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGGGCAGCAACCGGCGGGAATTGTCCCCGGAACCACGGCCCTTCGACGAAGCGGAAAAAAAGCCAACCGCGCAA‑CCCAAGCTGGGCGGAAAGCTGGCGATTTTCGCCGGGGATGTTTTCCAATGCCTGGGGtaat                                                                                                                                                                                                                    >  2:528968/1‑251 (MQ=255)
                    gggATCGTTGTCGTAACCCAGCCCCTTATTGATTGCCAGCACGATAGTGGTTGCCGGCGGACCGCCGCCCAGGGAGAGCACGGGGGCGAAGCTGGCGAAACAGAGCAGAAAGATAAGTGCAGCAACCGGCGGGATTTGCCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGTGAATGTTTTCCAGTGCCTGAAGtaata                                                                                                                                                                                                                   <  2:90608/250‑1 (MQ=255)
                    aGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCATGAGGACCGCCGCCCAACGATAGCACGGTGGCGAAGCTGGCGAAACAGAGAATAAAGATAAGTGCAGCAACCGGCGCGAAATGACGCCGTAACCACGGACATTCGACGAAGCGGAAAAAATGCCAACCGCGCCT‑CCCAAGATCGGCAGCAAGCTGACGTTGTACGCCGGGAATGTTTTCCCGTGCCTGCAGtcaa                                                                                                                                                                                                                    >  1:103455/1‑247 (MQ=255)
                               cgtcaCTCAGCGCCTGTTAGATTGCCAGACCGGTAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGGGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGCGTAATAAGCGGCTCGc                                                                                                                                                                                                         <  2:602324/246‑1 (MQ=255)
                               cGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGcc                                                                                                                                                                                                        <  1:56112/251‑1 (MQ=255)
                                gTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGGTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGcc                                                                                                                                                                                                        >  1:290057/1‑250 (MQ=255)
                                  aaCTCAGCGCCTGATAGATTGCCAGCTCGATAGGGGTTGCCTGCGGACCGCCGCACAGCGATATCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGccgtc                                                                                                                                                                                                     <  2:522136/251‑4 (MQ=255)
                                      cAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCggta                                                                                                                                                                                                 <  1:237307/251‑3 (MQ=255)
                                       aGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCggtag                                                                                                                                                                                                <  1:319798/251‑4 (MQ=255)
                                         cgcCTGATAGACTGCCAGCTCTATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGAt                                                                                                                                                                                              <  1:230826/251‑1 (MQ=255)
                                           ccTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTa                                                                                                                                                                                            <  1:403253/251‑1 (MQ=255)
                                           ccTGATAGATGGCCAGCGCGAGAGTGGTTGCCGGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGATACAGAGCAGAAAGATAAGTGCAGCAACCGGCGGGATGTGGCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTa                                                                                                                                                                                            <  2:91573/251‑1 (MQ=255)
                                              gatttaTTGCCAGCCCGATGGTGGGTGCCTGCGGACCGTCGTCCAGCGAAAGCACGGTGGCGAAGCTGGCGAAACAAAGCATAAAGAGAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCGGCTCAT‑CCCAAGCGGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCAACGGTAGATTaa                                                                                                                                                                                           <  2:210020/244‑1 (MQ=255)
                                                 aGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaa                                                                                                                                                                                       <  1:276400/250‑1 (MQ=255)
                                                             cGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatacgtgcgcca                                                                                                                                                                            <  1:529171/249‑12 (MQ=255)
                                                                       tgCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCt                                                                                                                                                                 <  1:618329/249‑1 (MQ=255)
                                                                        gCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGGCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCtt                                                                                                                                                                <  2:691432/250‑1 (MQ=255)
                                                                             cGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGc                                                                                                                                                          <  1:591415/251‑1 (MQ=255)
                                                                              ggACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAgggc                                                                                                                                                         >  1:453057/1‑249 (MQ=255)
                                                                                   gccgcccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAgg                                                                                                                                                    <  1:552105/247‑1 (MQ=255)
                                                                                     cgcccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGc                                                                                                                                                  >  1:137681/3‑251 (MQ=255)
                                                                                        ccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGaa                                                                                                                                               <  1:245428/251‑1 (MQ=255)
                                                                                                aGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCa                                                                                                                                       <  1:230898/251‑1 (MQ=255)
                                                                                                 gCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCAc                                                                                                                                      <  1:436139/251‑1 (MQ=255)
                                                                                                     ggTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGGGCGCCAGCAAAAAACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTcc                                                                                                                                   >  1:389005/1‑250 (MQ=255)
                                                                                                         gCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGAcc                                                                                                                              <  1:433258/251‑1 (MQ=255)
                                                                                                         gCGAAGCTGGAGAAACAGAGCAAAAAGCTAAGTGCAGCACCCGCTGGGACTTGTATCCGTAACTACGGCCATTCGACGAAGCGGAAAAAATGCCAGACGCTCAT‑CCCAAGCTGGGCGGAAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCCGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGGGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGAcc                                                                                                                              <  2:547762/251‑1 (MQ=255)
                                                                                                          cGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCg                                                                                                                             >  1:629466/1‑251 (MQ=255)
                                                                                                           gAAGCGGGCGAAACAGAGTATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGa                                                                                                                            <  2:1662/251‑1 (MQ=255)
                                                                                                            aaGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAACGGGAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGAGGCAAGCTGGCGTTGGTTGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGGAGATTAAAAAATACGGGGGCCAGAAAACTACCTTGCAGGCCGTAGGGAGAAAAGGACCACTCCAGACCgc                                                                                                                            >  2:530677/1‑249 (MQ=255)
                                                                                                               cTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCg                                                                                                                         >  1:438177/1‑250 (MQ=255)
                                                                                                               cTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAACCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCg                                                                                                                         >  1:677832/1‑250 (MQ=255)
                                                                                                                     aaaCAGAGCATAAAGATAAGTGCAGCAACCGACGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTACAAAATACGTGCGACAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGc                                                                                                                   >  1:253347/1‑250 (MQ=255)
                                                                                                                      aaCATAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCg                                                                                                                                                 <  2:565007/219‑1 (MQ=255)
                                                                                                                      aaCAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCg                                                                                                                                                 >  1:565007/1‑219 (MQ=255)
                                                                                                                      aaCAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCCCACCAGACCGAGCGATTGGc                                                                                                                   >  1:419645/1‑249 (MQ=255)
                                                                                                                         agagCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAgt                                                                                                              <  1:574338/251‑1 (MQ=255)
                                                                                                                           agCATAAAGATAAGGGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAgtgt                                                                                                            <  2:417088/251‑1 (MQ=255)
                                                                                                                             cATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAACGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCCCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGttt                                                                                                          >  1:324989/1‑250 (MQ=255)
                                                                                                                                    aTAAGTGCAGCAACCGGCGGGATTGTTCGCCGTACCCACGGCCCGTCGAAGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCCGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcca                                                                                                   <  2:281537/251‑1 (MQ=255)
                                                                                                                                      aagtgcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGccagccag                                                                                                  <  1:75404/246‑1 (MQ=255)
                                                                                                                                       agtgcagcaACCGGCGGGATTTGTCGCCGTAACCACGGACATTCGAAGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAAGGCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCagccagcc                                                                                                >  1:568662/4‑251 (MQ=255)
                                                                                                                                        gtgcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CACAAGCCGGGCGGCAAGCTGAAGGTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCCCCGGTAGATTAAAAAATACGTGCGCCAACAAAATAACTTGCAGGGCGGAGGGCCGAAAGGGCCACCCCAGACCGAGAGATTGGCAGAGGGTTGCCAGCCAAccc                                                                                               >  2:471904/3‑251 (MQ=255)
                                                                                                                                         tgcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCt                                                                                              >  1:337047/2‑251 (MQ=255)
                                                                                                                                          gcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCt                                                                                              <  1:687187/250‑1 (MQ=255)
                                                                                                                                          gcagcaACCGGCGGGATTGGGCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACGCCAGACCGAGAGATTGGCAGAGTGTTGCCAGCCAGCCCTg                                                                                             <  2:337047/251‑1 (MQ=255)
                                                                                                                                                   ggcggGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCaa                                                                                     >  1:418326/1‑249 (MQ=255)
                                                                                                                                                      gggATTTGTCGCCGTAACCACGGCCATTCGAAGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGCCCATAgcc                                                                                 >  1:177262/1‑249 (MQ=255)
                                                                                                                                                       ggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACg                                                                                <  1:309359/251‑1 (MQ=255)
                                                                                                                                                           ttGTCGCCTTAACCACGGTCAGGCGAGGAAGCGGAAAAAATGCCAGCCGCGGTA‑CCCAAGGTGGGCGGCAAGCTGTCGCTGTGAGCCGGGAATGTTTTCCAGTTCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATCCCTGGCAGGCCGGAGGGCGAAAAGGTCCACTCCAGACCGAGCGATGGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTa                                                                             <  2:309732/250‑1 (MQ=255)
                                                                                                                                                             gTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAAAGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGGGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGa                                                                          <  2:320175/251‑1 (MQ=255)
                                                                                                                                                                gCCGTAACCACGGACAGGCGACGAAGCGGAAAAAATGCCAGCCGC‑CGTGCCCAAGCAGGGCGGGAAGATGACGTTGGTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAAGAAGCGGCTCGCCATCGGTAGATTAAAAAAGGCGTGCGCAAGCAAAATACCTTGCAGGCCGGAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATg                                                                       <  2:174620/251‑1 (MQ=255)
                                                                                                                                                                gCCGAAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGGGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGGTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCGTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATg                                                                       <  2:103960/251‑1 (MQ=255)
                                                                                                                                                                cccGTACCCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATg                                                                       <  2:569945/250‑1 (MQ=255)
                                                                                                                                                                  cGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGcc                                                                     >  1:21613/1‑251 (MQ=255)
                                                                                                                                                                         ccGGCCATTCTGCGAAGCGTAAAAAATGCCAGCCGCGCCT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACGCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACa                                                              <  2:583716/250‑1 (MQ=255)
                                                                                                                                                                         aCGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAAc                                                               >  1:256970/1‑250 (MQ=255)
                                                                                                                                                                                 tCGACGAAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCa                                                       <  1:708495/250‑1 (MQ=255)
                                                                                                                                                                                 tCGACGAAGCGGAAAAAAGGCCAGCCGCGCAT‑CCCAAGCGGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCa                                                       <  2:438177/250‑1 (MQ=255)
                                                                                                                                                                                   gacgaAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGACTGTTTTCCAGTGCCTGGCGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGTCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCGGCCCCGGCGACCAACGCCGCTACGAATGCCGCCAACAGCAACCagcc                                                    >  2:243545/1‑248 (MQ=255)
                                                                                                                                                                                     cgaAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAACAATACGTGCGCCAGCAAAATACCTTGCAGGCCGAAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCCGGCGACCATAGACGCGAAGAATGCCGAAAACAGccaccagcccc                                                  >  2:597309/1‑246 (MQ=255)
                                                                                                                                                                                      gaAGCGGAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCATCACCg                                                 <  1:201602/251‑1 (MQ=255)
                                                                                                                                                                                        aGCGGAAAAAATGCCAGCCGCGCAT‑ACAAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGAAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCggg                                               >  1:513726/1‑251 (MQ=255)
                                                                                                                                                                                        aGCGGAAAAAAAGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTAGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTGAGAATGCCGAAAACAGCAACCAGCACCggg                                               <  2:316927/251‑1 (MQ=255)
                                                                                                                                                                                           ggAAAAAATGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGggagg                                            <  1:258924/251‑1 (MQ=255)
                                                                                                                                                                                             aaaaaaTGCCAGCCGCCCAC‑CCCAAGCTGGGGGGAAAGCCGACTGTGTTCGCCGGGAAGGTCTTCCTGGTCCTGGAGTAATAAGCTGCTCGCCATCGGTAGATTAAAAAATACGTGTGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACGCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGTAGGAt                                          <  2:442278/251‑1 (MQ=255)
                                                                                                                                                                                              aaaaaTGCCAGCCGCGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTGCGCCGGGAATGTTTTCCAGTGCCTGGAGTAAGAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGGAGGGCGAAAAGGACCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGAtt                                         <  2:256970/251‑2 (MQ=255)
                                                                                                                                                                                                  aTGCCAGCCGCGCAT‑CCCATGCTGGGCGGCACGCGGACGGTGTTCGCCGGGAATGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCATAGGTAGATTAAAAAATACGGGCGCCAGCAAAATACCGTGCAGGCCGTAGGGCGAAAAGGTCCACCCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAgg                                     <  2:5534/251‑1 (MQ=255)
                                                                                                                                                                                                     ccagccgcGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTc                                   >  1:238672/8‑250 (MQ=255)
                                                                                                                                                                                                       agccgcGCAT‑CCCAAGCTGGGCGGCATGCTGACGTTGTTTGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCTTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCAtt                                <  2:568662/246‑1 (MQ=255)
                                                                                                                                                                                                         ccgcGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGCTTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGCTTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGCCGCTCAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGc                              >  2:315983/4‑251 (MQ=255)
                                                                                                                                                                                                          cgcGCAT‑CCCAAGCTGGGCGGCAAGCTGCCGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAg                                                                                                                           >  1:108543/3‑157 (MQ=255)
                                                                                                                                                                                                          cgcGCAT‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGcg                             <  1:322462/249‑2 (MQ=255)
                                                                                                                                                                                                          cgcGCAG‑CCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAAGACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAg                                                                                                                           <  2:108543/155‑1 (MQ=255)
                                                                                                                                                                                                              cAT‑CCCAAGCTGGGCGGCAAGCGGACGTTGTTCGCTGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATGGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGcgcaca                         <  1:415966/251‑1 (MQ=255)
                                                                                                                                                                                                                    cAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGc                    <  1:16420/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGTGTAATAAGCGGCTCGCCCTCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTATTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGCGGATTAAGGACATTGCGCACAGACGcc                   >  1:326587/1‑250 (MQ=255)
                                                                                                                                                                                                                     aaGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCa                   <  1:385625/251‑1 (MQ=255)
                                                                                                                                                                                                                             cggcAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGcc           >  2:220754/1‑251 (MQ=255)
                                                                                                                                                                                                                               gcAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGGGCCTGGAGTAAGAAGCGGCTCGCCATCGGTAGATTAAAAAATACGGGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAg         <  2:374385/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aatctgaggtggtccGCCGGGAATGTTTTCCATGGCCCGGCGTAACAAGCGGCTCGCCATCGGTAGATTAAAAAATCCGCGCGCCAGCAAAATCCCTTGCAGGCCGTAGGGCGAAAAGTTCCCCCCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGATTAAGGTCATTGCGCACAGACGCAACAGCGCCAGcc       <  2:509117/237‑1 (MQ=255)
                                                                                                                                                                                                                                      gACGTTGTTCGCCGGGAATGGTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGCATGGCAGAGTGATGCCAGCACGCCCCGGCGACCATAGACGCTAGGAATGCCGAACACAGCAACCAGCACCCGGGGGAGTTAGGGCATTGCGCACAGACGCAACCGCGCCAGCCGAccc  >  2:221399/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
GCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCAT‑CCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCC  >  NC_000913/73782‑74261

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