Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,430 T→C A120A (GCA→GCG)  ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,4300TC98.7% 227.7 / ‑6.4 74A120A (GCA→GCG) ytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (38/35);  total (39/35)
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.91e-01

GCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGAT‑CGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGC  >  NC_000913/4624187‑4624676
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
gCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCttttt                                                                                                                                                                                                                                                  <  1:578115/251‑1 (MQ=255)
     cGGTTAACAGGAAGGGCGATTGCTGCCAGCTGGTACGTTTTCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCTCCTGTTGGGCTTCGGTGTCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCg                                                                                                                                                                                                                                             <  2:551545/251‑1 (MQ=255)
      ggTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATTTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGt                                                                                                                                                                                                                                            <  1:409990/251‑1 (MQ=255)
       gTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGGGTTATCCGCCTGTGGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTGTTCCCGG‑CAATTGGTTCGACAATCGGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTc                                                                                                                                                                                                                                           <  2:260627/251‑1 (MQ=255)
           aCAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTTCCCTGTAGCAGCGTGCGGGTCAGCAATACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAgc                                                                                                                                                                                                                                       <  1:363806/251‑1 (MQ=255)
           aCAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAgc                                                                                                                                                                                                                                       >  1:339228/1‑251 (MQ=255)
            cAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAgcc                                                                                                                                                                                                                                      >  1:245846/1‑250 (MQ=255)
               gAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCATCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCa                                                                                                                                                                                                                                   <  2:341059/251‑1 (MQ=255)
                aatggCGATTGCTGCCAGCGGGTACGTTTACCCAGTAGCAGCGTGCGGTGCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCTCCGGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCGGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAg                                                                                                                                                                                                                                  <  2:44229/248‑1 (MQ=255)
                   ggCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCg                                                                                                                                                                                                                               <  1:94577/251‑1 (MQ=255)
                       aTTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGTTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTcg                                                                                                                                                                                                                           <  1:70545/251‑1 (MQ=255)
                             gcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCAc                                                                                                                                                                                                                      <  1:538844/250‑1 (MQ=255)
                                  cGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATGGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGAc                                                                                                                                                                                                                <  1:420143/251‑1 (MQ=255)
                                   gggTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCAAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGCGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCgccg                                                                                                                                                                                                               >  2:593832/1‑248 (MQ=255)
                                      tacGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACTCtt                                                                                                                                                                                                            >  2:31038/3‑251 (MQ=255)
                                        cGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTc                                                                                                                                                                                                          >  2:4114/1‑251 (MQ=255)
                                           ttACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGcc                                                                                                                                                                                                       >  2:417815/1‑251 (MQ=255)
                                                   tAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACggg                                                                                                                                                                                               >  2:315439/1‑249 (MQ=255)
                                                              gagTCAGCAATACGCGGATGGCCAGGGAAAGCAGCAAAAACAGGCGTAAAATGTTAGTGGTGTTAGCCGCCTGTTGGGCTAGGGTGGCGAGCGTATGGGTGGCGAGCGTCAGGCGCAGAAAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTGGACAATCTGCGGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGGCGCTTTTGCCAGAGCGGGCGATAAGAT‑Cg                                                                                                                                                                                   <  2:15285/249‑1 (MQ=255)
                                                                gTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGcc                                                                                                                                                                                 >  1:564608/1‑251 (MQ=255)
                                                                  cAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCTCCTGTTGGGCTTCGGTGGCGTGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCt                                                                                                                                                                                <  1:291877/250‑1 (MQ=255)
                                                                    gCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGt                                                                                                                                                                             <  1:450294/251‑1 (MQ=255)
                                                                         acgccGATTGCCAGTGAACGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTc                                                                                                                                                                        <  2:148211/248‑1 (MQ=255)
                                                                            ccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCAt                                                                                                                                                                      >  2:556579/1‑250 (MQ=255)
                                                                            ccGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATa                                                                                                                                                                     <  1:538445/251‑1 (MQ=255)
                                                                               aTTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGGGGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCTCTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATacac                                                                                                                                                                  <  2:299650/251‑1 (MQ=255)
                                                                                              agcaACATCACGCGTAAAATGTCAGTGGTGTTATCCACCTGTTGGCCTTCCGTGGCGAGCGTTTGGCTGTCGAGCGACAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTCAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGCACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGCGGac                                                                                                                                                   >  2:284902/1‑250 (MQ=255)
                                                                                                 aaCATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACg                                                                                                                                                >  1:340058/1‑251 (MQ=255)
                                                                                                   cATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGc                                                                                                                                              <  2:197630/251‑1 (MQ=255)
                                                                                                    atcaggcgtatgctattAGTGGTGTTAACCACATGTTGGTCTTCGGTGGCGAGCAAATAGGTGCCGCGTGCCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTACAATAGCCGCCTGCTTTTTTACAGTCGAGCGCCAGCCGGTCGCGCACTTCCACGCTTTAGCCAGAACGGGCGGTGAGAT‑CGCCTTGTTCGTCATACTCGCCCGCGTCGAgaaaaaggct                                                                                                                                             >  2:233558/16‑242 (MQ=255)
                                                                                                    aTCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCt                                                                                                                                             >  1:369458/1‑251 (MQ=255)
                                                                                                             aaaaTGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGCTAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGt                                                                                                                                    >  1:318046/1‑251 (MQ=255)
                                                                                                              aaaTGTTAGAGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CACTTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGCTCCGGCTTTCATCCGtt                                                                                                                                   >  2:303640/1‑251 (MQ=255)
                                                                                                                 tGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAAc                                                                                                                                >  1:289731/1‑251 (MQ=255)
                                                                                                                    tAGTGGTGTTATCCTCCTGTGGGGCGTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGTTCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCGTGCGCCAGATCGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTg                                                                                                                              <  1:592049/250‑1 (MQ=255)
                                                                                                                     aGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTGCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTGCGCCAGAACGGTCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGGCGAGGATACGGCTTTCATCCGTTAACTGa                                                                                                                             <  2:592770/250‑1 (MQ=255)
                                                                                                                       tggtgTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACCCGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCg                                                                                                                          >  1:81033/1‑251 (MQ=255)
                                                                                                                         gtgtTATCCACCTGTTGGGCTTCGGTGGCGAGCGTGTGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGCCCGTTTTTCCCGG‑CAATTGGTTCGCCCATCTGCTGGTTAAAACAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATCCGGCTTTCATCCGTTAGCTGATCg                                                                                                                          >  2:443354/1‑249 (MQ=255)
                                                                                                                           gtTATCCACCTGTTGGGCGTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCGCTTCGACGCTTTCGCCAGAACGGGCGAGAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGa                                                                                                                      <  1:47182/251‑1 (MQ=255)
                                                                                                                                 cACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCCGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGTCAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGTGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAAAT‑CGCCTTGTTCGTCATACACGCCAGCGGCGAGGATACGGCTTTCATCAGTTAACTGATCGAGGATCGCCt                                                                                                                >  1:347891/1‑251 (MQ=255)
                                                                                                                                  aCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTg                                                                                                               >  1:422245/1‑251 (MQ=255)
                                                                                                                                      gttgagctgcggggacGAGCGTATGGGTTTCGAGCGTCAGGCTCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGACCGGGCGATAAGAT‑CACCTTTTTCGTCATACACGACCGAGTCGAGGATATGGCTTTCATCCGTAAATTGATTGAGGGTAGCCTGaagc                                                                                                           >  2:246983/16‑249 (MQ=255)
                                                                                                                                      gTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGAGCGAGGATCGCCTGAAt                                                                                                            <  1:130765/250‑1 (MQ=255)
                                                                                                                                      gTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGGCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATg                                                                                                           <  2:536734/251‑1 (MQ=255)
                                                                                                                                          ggCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCCGAACGGGCGCTAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATg                                                                                                           >  2:300392/1‑247 (MQ=255)
                                                                                                                                          ggCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGtt                                                                                                       <  1:344735/251‑1 (MQ=255)
                                                                                                                                           gCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGttt                                                                                                      >  2:582957/1‑251 (MQ=255)
                                                                                                                                             ttCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTACCTGATCGAGGATCGCCTGAATGCGttttt                                                                                                    >  1:213907/1‑251 (MQ=255)
                                                                                                                                                      ggtcgtctGGGGGTCGAGCGGCAGGCGCGGATAGCCGAGCGGTCCGTTTT‑GCCCGTCATTTGGGTCGACAATCGGCTGGTTAAAAGAGCCTCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGGCGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCGTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTATCTGATCGAGGATCGCCTGTATGCGTTTTTCAGCCggcg                                                                                           <  2:437300/244‑3 (MQ=255)
                                                                                                                                                               ggTGTCGAGAGTCTGGCGCAGACAACCGCACGGTCCGTTTTTTCCGG‑CAATTGGTTCGACAATCTGCTGGTCAAAATAGCCGCCCGATTTTTTACAGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCGTTTGCCCGTACGGGCGGTA‑GATTCGACTTGTTGGTCATACACGCCCGCGTCACGGATACGGCATTCAACCGTTAACTGCTCGCGGGTCGCCGGACAGCGGTTTTCTCCCGCTGAGTCgggac                                                                                  >  2:405410/1‑248 (MQ=255)
                                                                                                                                                                 tgtCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGc                                                                                >  1:183419/1‑251 (MQ=255)
                                                                                                                                                                     gagcgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATc                                                                             >  1:108660/5‑250 (MQ=255)
                                                                                                                                                                        cgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCg                                                                         >  1:493896/2‑251 (MQ=255)
                                                                                                                                                                        cgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCg                                                                         <  1:36977/250‑1 (MQ=255)
                                                                                                                                                                                     tAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTc                                                              >  1:241993/1‑249 (MQ=255)
                                                                                                                                                                                        ccGAGCCGTCCGTTTTTCCCGG‑CAATTGGTTCGGCAATCTCCTTGTTAAAATAGCCGCCCTCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTGCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTGCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAGCGTTCAGCGt                                                         <  2:81033/251‑1 (MQ=255)
                                                                                                                                                                                          gAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATACGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGtt                                                        >  1:589507/1‑249 (MQ=255)
                                                                                                                                                                                                   gTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCCCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGAcgcg                                              <  2:371126/251‑1 (MQ=255)
                                                                                                                                                                                                     tttttcccgacAATTGGTTCGACAATCTGATGGCTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACgcgc                                             >  1:492503/11‑251 (MQ=255)
                                                                                                                                                                                                       ttcccgg‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCCCCAGCCGGTCGCGCACTTCGACTCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGTTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTTCCTGACGCGCCa                                           <  2:245846/244‑1 (MQ=255)
                                                                                                                                                                                                            gg‑CAATTGGTTCGACCATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGGACGCATCAGCGGTGCAACGTTAAGCGTTACCTGACGCGCCAGGGTGc                                     >  1:532940/2‑251 (MQ=255)
                                                                                                                                                                                                                aaTTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCggg                                  <  1:180999/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggttCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCCTCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGcc                                >  2:217185/4‑249 (MQ=255)
                                                                                                                                                                                                                       tCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGttc                           >  1:494026/1‑251 (MQ=255)
                                                                                                                                                                                                                        cGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGttct                          >  1:416659/1‑251 (MQ=255)
                                                                                                                                                                                                                          aCAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGttcttc                        >  2:512120/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gctgGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCCGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCtg                  <  2:13780/250‑1 (MQ=255)
                                                                                                                                                                                                                                    ggTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCtgtg                <  2:289731/249‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaaTAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTAc          >  2:237603/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaaaGAGCCTCCCGCTTTTTTCCCGTCGAGCGCCAGCTGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTGAACTGATCGAGGATCGCCTGAATGCGTTTGTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGGTACCTGACGCGCCATGGTGCGGGCCAGTTCTTCCAGCTGTGGATTAc          <  2:527475/251‑1 (MQ=255)
                                                                                                                                                                                                                                             atccgcccgcTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCAGCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCTGGCCAGTTCTTCCAGCTGTGGATTACgttg      <  1:433992/248‑4 (MQ=255)
                                                                                                                                                                                                                                             aGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACTCTCTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACgttg      <  2:339228/250‑4 (MQ=255)
                                                                                                                                                                                                                                             aGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGt     >  1:245968/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ccgcccgcTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCg   <  1:303640/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cgccCGCTTTTTTACCGTCGAGCTCCAGCCGGGCGCGCACTTCGACGCTTGCGCCAGAACGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGc  <  2:518113/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cgccCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGACCGGGCGATAAGAT‑CGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGCCGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTAAGTTGTCGc  >  2:397869/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
GCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGAT‑CGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGC  >  NC_000913/4624187‑4624676

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