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

ATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCC  >  NC_000913/3822922‑3823407
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aTCCACCCCATCAAAACCGAAATCGAAGAGCTGGGTTGTGAGGCGCTGTAACCCAACCGGTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTTCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGAGCAt                                                                                                                                                                                                                                             <  2:142599/251‑1 (MQ=255)
       acaTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATc                                                                                                                                                                                                                                      >  1:529451/1‑251 (MQ=255)
                ccGAACTCGAAGAGCTGGGATTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACAGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGGGCTCAAAGAGCAGCGTTTTTACTCGATCAGGGACATCTACGCTTTc                                                                                                                                                                                                                             >  2:675961/1‑251 (MQ=255)
                  gAACTCGAAGAGCTGGGTTTTGAGGCGATGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCcg                                                                                                                                                                                                                           >  1:593425/1‑251 (MQ=255)
                      tCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGATTTCCGCGTg                                                                                                                                                                                                                       >  1:122992/1‑251 (MQ=255)
                      ccGAAGAGCGGGGTGTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTg                                                                                                                                                                                                                       <  2:593425/250‑1 (MQ=255)
                       cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGATTGCAGGAAGCGGGCATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTg                                                                                                                                                                                                                       >  2:207856/1‑250 (MQ=255)
                          agagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGTAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGGTCg                                                                                                                                                                                                                   >  2:483017/1‑251 (MQ=255)
                           gagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGAAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGt                                                                                                                                                                                                                  <  1:335514/251‑1 (MQ=255)
                             gCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCa                                                                                                                                                                                                                <  1:112298/251‑1 (MQ=255)
                                    tttGAGGCGCGGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAtt                                                                                                                                                                                                          <  2:477368/250‑1 (MQ=255)
                                      tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCt                                                                                                                                                                                                        <  1:284741/250‑1 (MQ=255)
                                          gcgcTGTATCCCAACCGTTATCGCGTCATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGacac                                                                                                                                                                                                   >  1:194086/1‑251 (MQ=255)
                                                    ccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATGCAGAAGATTCTTTCGGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACTTGTTAgcgc                                                                                                                                                                                         <  2:231268/251‑1 (MQ=255)
                                                     cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCATGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCACAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGGTCGGCAATGATTCTGACCCCTGTTCTcgcg                                                                                                                                                                                        >  2:268040/1‑251 (MQ=255)
                                                     cAACCGTTATCGCGTAATCAAAGAAGGGGTGAAAGACGCGCGCGGCAACCGTAAAGAGATGAGCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcgcg                                                                                                                                                                                        <  2:5208/251‑1 (MQ=255)
                                                         cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAACTGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATgg                                                                                                                                                                                                                                           >  2:282115/1‑196 (MQ=255)
                                                         cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATgg                                                                                                                                                                                                                                           <  1:282115/196‑1 (MQ=255)
                                                          gttatcttgTAATCAAAGAAGGGGTGAAAGCCGCTGGCGGATACCGTAAAGAAATAATCCAGAAGTTTCTTTCAGAAATCGAAGGGCGTTTGCAGGAAGCTGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACGCCTGTTATCGCGTgttg                                                                                                                                                                                   <  2:464929/243‑4 (MQ=255)
                                                            tATCTCGTAATCAAAGTATTGGGTACAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTGCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTGTCCTCGTGATCGTCAATGATTCTGACCCCTTTTAGCGCGTGCTgg                                                                                                                                                                                  <  2:122992/250‑1 (MQ=255)
                                                               cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCa                                                                                                                                                                              <  1:683935/251‑1 (MQ=255)
                                                               cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCa                                                                                                                                                                              <  1:390586/251‑1 (MQ=255)
                                                                    aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGc                                                                                                                                                                         >  1:519460/1‑251 (MQ=255)
                                                                     aTCAAAGATGTGGGGAAAGCCGCGCGGGGCAACAGTATAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATTCCGGGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGGGCTCAAAGAGCAGCGTTTTTAATCGATCATGGACATCTATGCTTTCCGCGGGATCGTCAAGGATTCTGACACCTGTTATCGCGTGCTGGGCCTGATGca                                                                                                                                                                        <  2:15499/251‑1 (MQ=255)
                                                                        aaaGAAATGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGATTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGCGACGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAACGAGCAGCGTTTTTACTCGATCATGGCCATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTCTCGCGTGCTTGGCCAGATGcac                                                                                                                                                                       >  2:224640/1‑249 (MQ=255)
                                                                           gaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCt                                                                                                                                                                  <  2:184608/251‑1 (MQ=255)
                                                                                  tgAAAGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACaa                                                                                                                                                            >  1:194599/1‑250 (MQ=255)
                                                                                  tgAAAGCCGCGCGCGGAAGCCGTAAAGAGATAATCCAGCAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTAGTCGCGCGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGCTCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCCCAGCCTGTACaa                                                                                                                                                            >  2:210407/1‑250 (MQ=255)
                                                                                      aGCCGCGCGCGGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCcgc                                                                                                                                                       >  2:691684/1‑251 (MQ=255)
                                                                                            gcgcgGCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCgg                                                                                                                                                 <  1:144638/251‑1 (MQ=255)
                                                                                               cggaaaCCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAGGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                                                                              <  2:547834/247‑1 (MQ=255)
                                                                                                 gCAACCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgc                                                                                                                                            >  1:552393/1‑251 (MQ=255)
                                                                                                    aCCGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTg                                                                                                                                         >  1:774/1‑251 (MQ=255)
                                                                                                     ccGTAAAGAGATAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTg                                                                                                                                         >  1:116827/1‑250 (MQ=255)
                                                                                                             agaTAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGACAGACtat                                                                                                                                >  1:201966/1‑251 (MQ=255)
                                                                                                             agaTAATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGAGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCCCAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAACGACtat                                                                                                                                >  2:521170/1‑251 (MQ=255)
                                                                                                               ataatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTCCTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACCCCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                                               <  2:452210/231‑1 (MQ=255)
                                                                                                               ataatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGc                                                                                                                                               >  1:452210/4‑234 (MQ=255)
                                                                                                                   tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGt                                                                                                                                                                                                                        >  2:542114/1‑157 (MQ=255)
                                                                                                                   tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGt                                                                                                                                                                                                                        <  1:542114/157‑1 (MQ=255)
                                                                                                                           aTTCTTTATGAAATCGAAGGGCGTTGGCAGGCAGCGGGAATAACTTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGGGCTCAATGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCTGCCATGATTCTGACAACTGTTATCACGTGCTGGGCCAGATGAACGGCCTGTACAAGCCGCGTCCGCGCCGCGTCAAACACGATCTCGCCATTcacc                                                                                                                   >  2:522884/1‑247 (MQ=255)
                                                                                                                             tctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaaaa                                                                                                                 <  1:591829/250‑1 (MQ=255)
                                                                                                                             tctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGCTTTACTGCAAAATGGTGCTCAAAGCGCAGCGTTTTTACTCGATAATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATAGCGTGCTGGGCCAGATGCACAGCCTGTACACGCCGCGTCCGGGCCGCATGAAAGACTATATCGCCATTCCaaa                                                                                                                  >  2:480394/1‑249 (MQ=255)
                                                                                                                                 tatGAAATCGAAGGGCGTTCGCAGATGGCGCGAATCACGTGCCGCGTCAGTGGGCGCGAGAAGCATCTTTATTCGCTTTACTCCAAAATGGGGTTCAAAGAGCAGCGTTCTTCCTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCTATGATTCTGACACCTGTTATCGCGTGCTGGGCAAGATGACCAGCCTGTACGAGCTGCGCCCGGGACGCGTGAACGACTATATCGCCATTCCAAACGCGa                                                                                                             >  2:237753/3‑250 (MQ=255)
                                                                                                                                      aaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCt                                                                                                       >  1:315960/1‑251 (MQ=255)
                                                                                                                                          gAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATc                                                                                                    >  2:333111/1‑250 (MQ=255)
                                                                                                                                                    tGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAGGACTCTATCGCCATTCCAAAAGAGAACGGCTATAAGTCTTTGcac                                                                                         >  2:183754/1‑251 (MQ=255)
                                                                                                                                                      cAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCacac                                                                                       >  1:591669/1‑251 (MQ=255)
                                                                                                                                                              gggAATTCCGTGCCGCGGCAGTGGTCGCGAGAAGCCTCTTTATGCGATTTACTTCAAAATGGGGCTCAAAGAGCAGCGTTTGTACTCGATCATGGACATCTTCGCTTCCCGCGTGAGCGTCAAGGATTCTGACCCCTGTTAGCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCCTTGCACACCTCgatga                                                                               <  2:411869/251‑1 (MQ=255)
                                                                                                                                                               ggAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatgat                                                                              >  1:675634/1‑251 (MQ=255)
                                                                                                                                                                 aaTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                                                                            <  1:196418/251‑1 (MQ=255)
                                                                                                                                                                  aTACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCgg                                                                           <  2:102141/251‑1 (MQ=255)
                                                                                                                                                                            gcgTCAGTGGTCGCGAGAAGCATTTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAGTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGgt                                                                 <  2:675634/251‑1 (MQ=255)
                                                                                                                                                                                cAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                                              <  1:691684/250‑1 (MQ=255)
                                                                                                                                                                                 agtgggctcgAGAAGCATCTTTTTTCGATTTACGGCAAAATGGTTCACAAAGAGCAGCGGTTTTACGTGCACATGGGCATCTACGGGTTACGGGTGATCGGTAATGTCTCTATCACCTGTTATCGCGTGCTGTGCCAGAGGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGACTTTGCACCCCTCGATGATCGGCCCGCACGGTGTGca                                                            <  2:434146/243‑2 (MQ=255)
                                                                                                                                                                                 agtgggatcgAGAAGCATCTTTATTCGATTTATTGCAAAAGGGTGCTCAAAGAGCAGTGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTAGCGCGTGCTGGGCCAGATGCACAGCTGGTACAAGCCGCTTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                                            <  2:434172/243‑1 (MQ=255)
                                                                                                                                                                                 aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                                            >  1:275344/1‑251 (MQ=255)
                                                                                                                                                                                    ggTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGGGCCGGt                                                         >  1:127029/1‑251 (MQ=255)
                                                                                                                                                                                        gcgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGCAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGGGTGCCGGTTg                                                       >  2:477136/1‑249 (MQ=255)
                                                                                                                                                                                         cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                                                                                                                                                                                >  2:668891/1‑111 (MQ=255)
                                                                                                                                                                                         cgAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                                                                                                                                                                                <  1:668891/111‑1 (MQ=255)
                                                                                                                                                                                          gagaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAggc                                                   >  1:134617/1‑250 (MQ=255)
                                                                                                                                                                                               gcaattttATTCGAATTACTGCAAAGTGGCGCTCAAAGAGCAGCGTTTTTCCTCGATCCTGGACATATACGCTTTCCGCGTGATCGTCAATGATTCCGCCACCTGTTATCGCGTGCTGGGCCAGAAGCACAGCCTGTACAAGCAGCGTCCGTGCTGCGTGAAAGACTATATCACCATTCACCAAGCGAGCGGCTATCAGTCTTTGCACACCTCGATGATCGACCAGAACGGTGTGCCGGTTGAggaacagg                                              >  2:208867/6‑245 (MQ=255)
                                                                                                                                                                                                cacttttatcagATTTACTGAAAAATGGTGCTCAAATAGCAGCGTTTTTACTCGATCATGGACAGCTACGGTTGCCGCGTGATCGGAAAAGATTCTGCCACTTGTTATCGCGGGCGGGGCCATATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTGTGCACACCTCGATGAGCGGCCCGCACGGTGTGACGGTTGAGGTCCAGAt                                             <  2:529451/240‑1 (MQ=255)
                                                                                                                                                                                                   cTTTATTTGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATc                                            <  1:247671/249‑1 (MQ=255)
                                                                                                                                                                                                   cTTTATTCGATTTACTGCAAAATGGGGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTACGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCCAAAGCGAACGGCTATCCGTCTTTGCACCCCTCGATGATCGGCCCGCCAGGTGAGCCGGTTGAGGTCCCGATCCg                                          >  2:570942/1‑251 (MQ=255)
                                                                                                                                                                                                    tttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                                         >  1:157575/1‑251 (MQ=255)
                                                                                                                                                                                                            aTTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa                                    >  1:637828/1‑248 (MQ=255)
                                                                                                                                                                                                             tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAgct                                >  2:65043/1‑249 (MQ=255)
                                                                                                                                                                                                                aCTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGCTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGCCACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGCAGATATg                             >  2:220833/1‑251 (MQ=255)
                                                                                                                                                                                                                 cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGACAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGCTCCGTACCGACGATATgg                            >  1:147362/1‑251 (MQ=255)
                                                                                                                                                                                                                  tGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCATTGCTTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATAACGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGAATGCCGGTTGAGGTCCAGATCCGACCCGAAGATATgg                            >  2:406975/1‑250 (MQ=255)
                                                                                                                                                                                                                  tGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAACCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGGTGAGGTCCAGCTCCGTACCGAAGATATggg                           >  2:250676/1‑250 (MQ=255)
                                                                                                                                                                                                                       aaTGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTCCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGAGGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTCGAGGTCCAGATCCGTACCGAAGATATGGACCAGa                      <  1:107916/251‑1 (MQ=255)
                                                                                                                                                                                                                         tGGTGCTCAAAGCGCAGCGTTTTTACTCGTTCATGGACATCTACGCTTTCCGCGAGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTAAATCGCCATTCCAAAAGCGAACGGCTCTCAGTCTTTGCACACCCCGATGATCGGCCCCCACGGTGTGCCGGTTGATGTCCAAATCCGTACCGAAGATAGGGACCAggtg                    >  2:493510/1‑248 (MQ=255)
                                                                                                                                                                                                                          ggtgtaaaaaGAGCAGCGTTTTTACTCGATCATGGACATCGACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTACCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATg                    <  2:172475/243‑1 (MQ=255)
                                                                                                                                                                                                                           gagCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTATGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGGTGAGGTCCAGATCCGTACCGAGGATATGGACCAGATGGc                  >  2:525815/3‑251 (MQ=255)
                                                                                                                                                                                                                                cAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCgg                >  1:221137/1‑248 (MQ=255)
                                                                                                                                                                                                                                cAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGgaga             <  1:268040/251‑1 (MQ=255)
                                                                                                                                                                                                                                   agagCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATgg          <  1:615760/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gTTTTTACTCGATCATGGACATCTGCGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGCTGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCAAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGCGCTCCGGTTGCGGGCCAGATCCGAACCGAAGATATGACCCAGATGGCGGAGATGGGTGTTGcc  >  2:127567/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGAACCGAAGATATGGACCAGATGGCGGAGATGGGGGTTGcc  >  1:196357/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
ATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCC  >  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: