Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,162 C→T H255Y (CAC→TAC)  spoT → bifunctional (p)ppGpp synthase/hydrolase SpoT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,1620CT97.6% 267.2 / ‑4.3 82H255Y (CAC→TAC) spoTbifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (47/33);  total (48/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

ATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCC  >  NC_000913/3822922‑3823407
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aTCCACCCCATCAAAACCGAAATCGAAGAGCTGGGTTGTGAGGCGCTGTAACCCAACCGGTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTTCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGAGCAt                                                                                                                                                                                                                                             <  2:142599/251‑1 (MQ=255)
       acaTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATc                                                                                                                                                                                                                                      >  1:529451/1‑251 (MQ=255)
                ccGAACTCGAAGAGCTGGGATTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACAGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGGGCTCAAAGAGCAGCGTT‑TTTACTCGATCAGGGACATCTACGCTTTc                                                                                                                                                                                                                             >  2:675961/1‑251 (MQ=255)
                  gAACTCGAAGAGCTGGGTTTTGAGGCGATGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCcg                                                                                                                                                                                                                           >  1:593425/1‑251 (MQ=255)
                      tCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGATTTCCGCGTg                                                                                                                                                                                                                       >  1:122992/1‑251 (MQ=255)
                      ccGAAGAGCGGGGTGTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTg                                                                                                                                                                                                                       <  2:593425/250‑1 (MQ=255)
                       cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGATTGCAGGAAGCGGGCATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTg                                                                                                                                                                                                                       >  2:207856/1‑250 (MQ=255)
                          agagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGTAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGGTCg                                                                                                                                                                                                                   >  2:483017/1‑251 (MQ=255)
                           gagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGAAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGt                                                                                                                                                                                                                  <  1:335514/251‑1 (MQ=255)
                             gCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCa                                                                                                                                                                                                                <  1:112298/251‑1 (MQ=255)
                                    tttGAGGCGCGGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAtt                                                                                                                                                                                                          <  2:477368/250‑1 (MQ=255)
                                      tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCt                                                                                                                                                                                                        <  1:284741/250‑1 (MQ=255)
                                          gcgcTGTATCCCAACCGTTATCGCGTCATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGacac                                                                                                                                                                                                   >  1:194086/1‑251 (MQ=255)
                                                    ccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATGCAGAAGATTCTTTCGGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACTTGTTAgcgc                                                                                                                                                                                         <  2:231268/251‑1 (MQ=255)
                                                     cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCATGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCACAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGGTCGGCAATGATTCTGACCCCTGTTCTcgcg                                                                                                                                                                                        >  2:268040/1‑251 (MQ=255)
                                                     cAACCGTTATCGCGTAATCAAAGAAGGGGTGAAAGACGCGCGCGGCAACCGTAAAGAGATGAGCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcgcg                                                                                                                                                                                        <  2:5208/251‑1 (MQ=255)
                                                         cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAACTGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATgg                                                                                                                                                                                                                                           >  2:282115/1‑196 (MQ=255)
                                                         cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATgg                                                                                                                                                                                                                                           <  1:282115/196‑1 (MQ=255)
                                                          gttatcttgTAATCAAAGAAGGGGTGAAAGCCGCTGGCGGATACCGTAAAGAAATAATCCAGAAGTTTCTTTCAGAAATCGAAGGGCGTTTGCAGGAAGCTGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACGCCTGTTATCGCGTgttg                                                                                                                                                                                   <  2:464929/243‑4 (MQ=255)
                                                            tATCTCGTAATCAAAGTATTGGGTACAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTGCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTGTCCTCGTGATCGTCAATGATTCTGACCCCTTTTAGCGCGTGCTgg                                                                                                                                                                                  <  2:122992/250‑1 (MQ=255)
                                                               cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCa                                                                                                                                                                              <  1:683935/251‑1 (MQ=255)
                                                               cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCa                                                                                                                                                                              <  1:390586/251‑1 (MQ=255)
                                                                    aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGc                                                                                                                                                                         >  1:519460/1‑251 (MQ=255)
                                                                     aTCAAAGATGTGGGGAAAGCCGCGCGGGGCAACAGTATAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATTCCGGGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGGGCTCAAAGAGCAGCGTT‑TTTAATCGATCATGGACATCTATGCTTTCCGCGGGATCGTCAAGGATTCTGACACCTGTTATCGCGTGCTGGGCCTGATGca                                                                                                                                                                        <  2:15499/251‑1 (MQ=255)
                                                                        aaaGAAATGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGATTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGCGACGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAACGAGCAGCGTT‑TTTACTCGATCATGGCCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTCTCGCGTGCTTGGCCAGATGcac                                                                                                                                                                       >  2:224640/1‑249 (MQ=255)
                                                                           gaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCt                                                                                                                                                                  <  2:184608/251‑1 (MQ=255)
                                                                                  tgAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACaa                                                                                                                                                            >  1:194599/1‑250 (MQ=255)
                                                                                  tgAAAGCCGCGCGCGGAAGCCGTAAAGAGATAATCCAGCAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTAGTCGCGCGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGCTCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCCCAGCCTGTACaa                                                                                                                                                            >  2:210407/1‑250 (MQ=255)
                                                                                      aGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCcgc                                                                                                                                                       >  2:691684/1‑251 (MQ=255)
                                                                                            gcgcgGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCgg                                                                                                                                                 <  1:144638/251‑1 (MQ=255)
                                                                                               cggaaaCCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAGGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                                                                              <  2:547834/247‑1 (MQ=255)
                                                                                                 gCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgc                                                                                                                                            >  1:552393/1‑251 (MQ=255)
                                                                                                    aCCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTg                                                                                                                                         >  1:774/1‑251 (MQ=255)
                                                                                                     ccGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTg                                                                                                                                         >  1:116827/1‑250 (MQ=255)
                                                                                                             agaTAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGACAGACtat                                                                                                                                >  1:201966/1‑251 (MQ=255)
                                                                                                             agaTAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGAGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCCCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAACGACtat                                                                                                                                >  2:521170/1‑251 (MQ=255)
                                                                                                               ataatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTCCTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACCCCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                                               <  2:452210/231‑1 (MQ=255)
                                                                                                               ataatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                                               >  1:452210/4‑234 (MQ=255)
                                                                                                                   tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGt                                                                                                                                                                                                                        >  2:542114/1‑157 (MQ=255)
                                                                                                                   tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGt                                                                                                                                                                                                                        <  1:542114/157‑1 (MQ=255)
                                                                                                                           aTTCTTTATGAAATCGAAGGGCGTTGGCAGGCAGCGGGAATAACTTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGGGCTCAATGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCTGCCATGATTCTGACAACTGTTATCACGTGCTGGGCCAGATGAACGGCCTGTACAAGCCGCGTCCGCGCCGCGTCAAACACGATCTCGCCATTcacc                                                                                                                   >  2:522884/1‑247 (MQ=255)
                                                                                                                             tctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                                                 <  1:591829/250‑1 (MQ=255)
                                                                                                                             tctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGCTTTACTGCAAAATGGTGCTCAAAGCGCAGCGTT‑TTTACTCGATAATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATAGCGTGCTGGGCCAGATGCACAGCCTGTACACGCCGCGTCCGGGCCGCATGAAAGACTATATCGCCATTCCaaa                                                                                                                  >  2:480394/1‑249 (MQ=255)
                                                                                                                                 tatGAAATCGAAGGGCGTTCGCAGATGGCGCGAATCACGTGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGCTTTACTCCAAAATGGGGTTCAAAGAGCAGCGTTCTTC‑CTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCTATGATTCTGACACCTGTTATCGCGTGCTGGGCAAGATGACCAGCCTGTACGAGCTGCGCCCGGGACGCGTGAACGACTATATCGCCATTCCAAACGCGa                                                                                                             >  2:237753/3‑250 (MQ=255)
                                                                                                                                      aaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCt                                                                                                       >  1:315960/1‑251 (MQ=255)
                                                                                                                                          gAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATc                                                                                                    >  2:333111/1‑250 (MQ=255)
                                                                                                                                                    tGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAGGACTCTATCGCCATTCCAAAAGAGAACGGCTATAAGTCTTTGcac                                                                                         >  2:183754/1‑251 (MQ=255)
                                                                                                                                                      cAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCacac                                                                                       >  1:591669/1‑251 (MQ=255)
                                                                                                                                                              gggAATTCCGTGCCGCGGCAGTGGTCGCGAGAAGCCTCTTTATGCGATTTACTTCAAAATGGGGCTCAAAGAGCAGCGTT‑TGTACTCGATCATGGACATCTTCGCTTCCCGCGTGAGCGTCAAGGATTCTGACCCCTGTTAGCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCCTTGCACACCTCgatga                                                                               <  2:411869/251‑1 (MQ=255)
                                                                                                                                                               ggAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatgat                                                                              >  1:675634/1‑251 (MQ=255)
                                                                                                                                                                 aaTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                                            <  1:196418/251‑1 (MQ=255)
                                                                                                                                                                  aTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCgg                                                                           <  2:102141/251‑1 (MQ=255)
                                                                                                                                                                            gcgTCAGTGGTCGCGAGAAGCATTTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAGTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGgt                                                                 <  2:675634/251‑1 (MQ=255)
                                                                                                                                                                                cAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                                              <  1:691684/250‑1 (MQ=255)
                                                                                                                                                                                 agtgggctcgAGAAGCATCTTTTTTCGATTTACGGCAAAATGGTTCACAAAGAGCAGCGGT‑TTTACGTGCACATGGGCATCTACGGGTTACGGGTGATCGGTAATGTCTCTATCACCTGTTATCGCGTGCTGTGCCAGAGGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGACTTTGCACCCCTCGATGATCGGCCCGCACGGTGTGca                                                            <  2:434146/243‑2 (MQ=255)
                                                                                                                                                                                 agtgggatcgAGAAGCATCTTTATTCGATTTATTGCAAAAGGGTGCTCAAAGAGCAGTGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTAGCGCGTGCTGGGCCAGATGCACAGCTGGTACAAGCCGCTTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                                            <  2:434172/243‑1 (MQ=255)
                                                                                                                                                                                 aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                                            >  1:275344/1‑251 (MQ=255)
                                                                                                                                                                                    ggTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGGGCCGGt                                                         >  1:127029/1‑251 (MQ=255)
                                                                                                                                                                                        gcgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGCAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGGGTGCCGGTTg                                                       >  2:477136/1‑249 (MQ=255)
                                                                                                                                                                                         cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                                                                                                                                                                                >  2:668891/1‑111 (MQ=255)
                                                                                                                                                                                         cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                                                                                                                                                                                <  1:668891/111‑1 (MQ=255)
                                                                                                                                                                                          gagaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAggc                                                   >  1:134617/1‑250 (MQ=255)
                                                                                                                                                                                               gcaattttATTCGAATTACTGCAAAGTGGCGCTCAAAGAGCAGCGTT‑TTTCCTCGATCCTGGACATATACGCTTTCCGCGTGATCGTCAATGATTCCGCCACCTGTTATCGCGTGCTGGGCCAGAAGCACAGCCTGTACAAGCAGCGTCCGTGCTGCGTGAAAGACTATATCACCATTCACCAAGCGAGCGGCTATCAGTCTTTGCACACCTCGATGATCGACCAGAACGGTGTGCCGGTTGAggaacagg                                              >  2:208867/6‑245 (MQ=255)
                                                                                                                                                                                                cacttttatcagATTTACTGAAAAATGGTGCTCAAATAGCAGCGTT‑TTTACTCGATCATGGACAGCTACGGTTGCCGCGTGATCGGAAAAGATTCTGCCACTTGTTATCGCGGGCGGGGCCATATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTGTGCACACCTCGATGAGCGGCCCGCACGGTGTGACGGTTGAGGTCCAGAt                                             <  2:529451/240‑1 (MQ=255)
                                                                                                                                                                                                   cTTTATTTGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATc                                            <  1:247671/249‑1 (MQ=255)
                                                                                                                                                                                                   cTTTATTCGATTTACTGCAAAATGGGGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTACGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCCAAAGCGAACGGCTATCCGTCTTTGCACCCCTCGATGATCGGCCCGCCAGGTGAGCCGGTTGAGGTCCCGATCCg                                          >  2:570942/1‑251 (MQ=255)
                                                                                                                                                                                                    tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                                         >  1:157575/1‑251 (MQ=255)
                                                                                                                                                                                                            aTTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa                                    >  1:637828/1‑248 (MQ=255)
                                                                                                                                                                                                             tttACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAgct                                >  2:65043/1‑249 (MQ=255)
                                                                                                                                                                                                                aCTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGCTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGCCACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGCAGATATg                             >  2:220833/1‑251 (MQ=255)
                                                                                                                                                                                                                 cTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGACAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGCTCCGTACCGACGATATgg                            >  1:147362/1‑251 (MQ=255)
                                                                                                                                                                                                                  tGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCATTGCTTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATAACGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGAATGCCGGTTGAGGTCCAGATCCGACCCGAAGATATgg                            >  2:406975/1‑250 (MQ=255)
                                                                                                                                                                                                                  tGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAACCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGGTGAGGTCCAGCTCCGTACCGAAGATATggg                           >  2:250676/1‑250 (MQ=255)
                                                                                                                                                                                                                       aaTGGTGCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTCCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGAGGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTCGAGGTCCAGATCCGTACCGAAGATATGGACCAGa                      <  1:107916/251‑1 (MQ=255)
                                                                                                                                                                                                                         tGGTGCTCAAAGCGCAGCGTT‑TTTACTCGTTCATGGACATCTACGCTTTCCGCGAGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTAAATCGCCATTCCAAAAGCGAACGGCTCTCAGTCTTTGCACACCCCGATGATCGGCCCCCACGGTGTGCCGGTTGATGTCCAAATCCGTACCGAAGATAGGGACCAggtg                    >  2:493510/1‑248 (MQ=255)
                                                                                                                                                                                                                          ggtgtaaaaaGAGCAGCGTT‑TTTACTCGATCATGGACATCGACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTACCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATg                    <  2:172475/243‑1 (MQ=255)
                                                                                                                                                                                                                           gagCTCAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTATGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGGTGAGGTCCAGATCCGTACCGAGGATATGGACCAGATGGc                  >  2:525815/3‑251 (MQ=255)
                                                                                                                                                                                                                                cAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCgg                >  1:221137/1‑248 (MQ=255)
                                                                                                                                                                                                                                cAAAGAGCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGgaga             <  1:268040/251‑1 (MQ=255)
                                                                                                                                                                                                                                   agagCAGCGTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATgg          <  1:615760/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gTT‑TTTACTCGATCATGGACATCTGCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGCTGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCAAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGCGCTCCGGTTGCGGGCCAGATCCGAACCGAAGATATGACCCAGATGGCGGAGATGGGTGTTGcc  >  2:127567/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gTT‑TTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGAACCGAAGATATGGACCAGATGGCGGAGATGGGGGTTGcc  >  1:196357/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
ATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTT‑TTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCC  >  NC_000913/3822922‑3823407

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