Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,996,049 T→A L156* (TTG→TAG)  yigA → DUF484 domain‑containing protein YigA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,996,0490TA97.7% 274.3 / ‑2.6 87L156* (TTG→TAG) yigADUF484 domain‑containing protein YigA
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (41/43);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGACCGAT  >  NC_000913/3995805‑3996295
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
ttGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGc                                                                                                                                                                                                                                                    >  1:137392/1‑250 (MQ=255)
   aTGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGGCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAggg                                                                                                                                                                                                                                                >  2:607895/1‑250 (MQ=255)
        aCAGGCTATCGCCAACGTAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGCGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGACGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGGACAGc                                                                                                                                                                                                                                           >  2:99922/1‑251 (MQ=255)
          aGGCTATCGCCAACGAAGGCCTGTTTTAGCGCCTACTCTACCGGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCAc                                                                                                                                                                                                                                         <  2:137392/251‑1 (MQ=255)
            gCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTa                                                                                                                                                                                                                                       <  1:477141/251‑1 (MQ=255)
              tATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCCGCGTTAGGGGCAGGAACCGCAatatt                                                                                                                                                                                                                                     >  2:597513/1‑250 (MQ=255)
                ccgccatagaatgCCTGTTTTAACACTTACTCTACCTGCGGCGCCGGCTCCCCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCtt                                                                                                                                                                                                                                   <  2:631173/239‑1 (MQ=255)
                 cGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCtt                                                                                                                                                                                                                                   <  1:434308/250‑1 (MQ=255)
                   ccAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGAGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGTGGCAGGAACCGCACTATCTTggg                                                                                                                                                                                                                                >  1:606073/1‑251 (MQ=255)
                   ccAACGAAGGCCTGTTTTAGCGCCTCCTCTCCCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTggg                                                                                                                                                                                                                                <  1:527550/251‑1 (MQ=255)
                    cAACGACGGCCTTTTTTTTCTCCTTCTCTCCCTGCAGCTCAGTCTCACCGCCGCCAGCAGTCTCGACGATAGGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGc                                                                                                                                                                                                                               <  2:375976/251‑1 (MQ=255)
                        gAAGGCCTGTTTTATCGCTTACTCTACCGGCAGCGCAGTCTCCCCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTATGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCt                                                                                                                                                                                                                           <  1:563992/251‑1 (MQ=255)
                          aGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCtt                                                                                                                                                                                                                          >  1:355224/1‑250 (MQ=255)
                                     acctcttcctctcccggcaccgcATTCTCACCTCCGCCAGCAGGCTCGACGATATGCTGATGCTCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGTGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCaga                                                                                                                                                                                                               <  2:651660/230‑1 (MQ=255)
                                      tCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgc                                                                                                                                                                                                             >  1:2973/1‑251 (MQ=255)
                                       cGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCGCTATCTTGGGCCGCTTAACGGACCAGAgcc                                                                                                                                                                                                            >  1:393706/1‑250 (MQ=255)
                                        gCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgctg                                                                                                                                                                                                           <  1:478992/251‑1 (MQ=255)
                                        gCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAgctg                                                                                                                                                                                                           >  2:478992/1‑251 (MQ=255)
                                             ctctACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGCGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGGCCCGAGCTGCtggg                                                                                                                                                                                                      >  1:144788/1‑250 (MQ=255)
                                               ctACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGGGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGAGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGGGACGCtggtgg                                                                                                                                                                                                    >  2:97066/1‑251 (MQ=255)
                                                tACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGGGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGGCCAGAGCTGCtggtggg                                                                                                                                                                                                   >  2:41094/1‑250 (MQ=255)
                                                         cgcAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCggg                                                                                                                                                                                          >  1:46458/1‑250 (MQ=255)
                                                          gcAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCCGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGACCTATCTTGGGCCGCTTAAAGGACCAGAGCTGCTGGAGGTGCTACCGGa                                                                                                                                                                                          >  1:356917/1‑250 (MQ=255)
                                                           cAGTCTCACCGCCGCCAGCGGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGa                                                                                                                                                                                          >  1:279789/1‑249 (MQ=255)
                                                            aGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGACGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGGAGc                                                                                                                                                                                       >  1:489192/1‑251 (MQ=255)
                                                             gTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTAACTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTTCCGGAAGCg                                                                                                                                                                                      >  2:120987/1‑251 (MQ=255)
                                                              cctcACCGCCGCCCGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGGGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATGAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGa                                                                                                                                                                                     <  2:157439/250‑1 (MQ=255)
                                                                  aCCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGc                                                                                                                                                                                 >  1:642392/1‑251 (MQ=255)
                                                                   ccgccgCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCTAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGgtggtg                                                                                                                                                                                                  <  1:698277/233‑1 (MQ=255)
                                                                   ccgccgCCAGCAGTCTCGACGATATGCTGATGCGATTTCACCGCTGGGCACGCGATCTCGGCCTGGCAGGGGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGGCGAACCACACTCATCTGGCATTAAGCCGTCCGTCTTTCGGACCGCTGCGTATTCAGCGTTAGGGGCAGGAGCAGCACTATCTTGGGCCGCCTAAAGGACCAGAGCTGAtggtggggaggtcggcagggcgtcgg                                                                                                                                                                                >  2:698277/1‑231 (MQ=255)
                                                                         ccAGCAGGCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGa                                                                                                                                                                          <  1:99922/251‑1 (MQ=255)
                                                                          cagcagTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTAACGGAAGCGAAAGCGGTGGGAt                                                                                                                                                                         >  1:265874/1‑251 (MQ=255)
                                                                           agcagTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATc                                                                                                                                                                        >  1:159807/1‑251 (MQ=255)
                                                                           agcagTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATc                                                                                                                                                                        <  1:636484/251‑1 (MQ=255)
                                                                           agcagTCCCGGCGAGATGCTGATGCGCATTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCAGTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATc                                                                                                                                                                        <  2:472996/251‑1 (MQ=255)
                                                                             cagTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCgg                                                                                                                                                                      <  1:55077/251‑1 (MQ=255)
                                                                             cagTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGGGGGATCgg                                                                                                                                                                      >  1:568930/1‑251 (MQ=255)
                                                                             cagTCTCGACGATATGCTGATGCGCTTTCACCGCTGGCCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCgg                                                                                                                                                                      >  1:418999/1‑251 (MQ=255)
                                                                                    gacgaTATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGGGCGCCGTCGAACCACACTCATCTGTCATTAAGCCGACAGTCTTTCGGACCGGTGCGTATTCAGCGTTAGGTGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCt                                                                                                                                                                                                <  1:335577/218‑1 (MQ=255)
                                                                                      cgaTATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGTGAAAGCGGTGGGATCGGTGGCGATGt                                                                                                                                                             <  1:625217/251‑1 (MQ=255)
                                                                                            gCTGGTGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGTACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGCACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGCAGCGAATGCGGTGGGATCGGAGGCGATGTCGATGc                                                                                                                                                       >  1:545764/1‑251 (MQ=255)
                                                                                             cTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCt                                                                                                                                                      >  1:562593/1‑251 (MQ=255)
                                                                                               gATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGCCCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGGTCGGTGGCGAGGTCGGTGcagg                                                                                                                                                    >  2:511746/1‑248 (MQ=255)
                                                                                               gATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTgg                                                                                                                                                    <  1:446446/251‑1 (MQ=255)
                                                                                                                           ctcGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGAGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGCTGCTGATTTGGGGGGCGTGCt                                                                                                                        >  1:376536/1‑251 (MQ=255)
                                                                                                                           cgcGGCCTGGCAGGTGCGAGTCTGCGCCTGGTGCCGGATCGCGGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGc                                                                                                                         <  2:562593/248‑1 (MQ=255)
                                                                                                                              ggtcttgCAGTTGCGAGGCTGCTCCTGTTTCCGGATCCCGGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCAGTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGATCTGCTGTTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATTTCGATGCTGGGAAGCGATTCTGATTTGGGTGTCGTGCTGtt                                                                                                                     <  2:537200/245‑1 (MQ=255)
                                                                                                                               gCCTGTCAGGCGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTCAACCATCCGTCTTTCGAACCGCTGCGTCTTCAGCGTTAGGGGCAGGAACAGAACTATCTTGGGCCCCCTCACGGCCCAGAGCTGCTGGTGGTGCTACCGGAAACGCAAGCGGTGGGGACGGTTGCCATGTCGCTGCTGAGCAGCGATACTGTTTTGGgtgccgggcagttt                                                                                                                    >  2:51928/1‑240 (MQ=255)
                                                                                                                                ccttgCAGGTGCGAGTCGGCGCCTGGTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTTTTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCGGCTGGTGGTGCTACTGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTa                                                                                                                   <  2:302224/247‑1 (MQ=255)
                                                                                                                                    gCAGGTGAGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGAGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGCACCGCTGCGTATTCAGCGTTCGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGGGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGTTTTGGGTGTCGTGCTGTTTACCa                                                                                                                >  2:567329/1‑250 (MQ=255)
                                                                                                                                          gcgcgTCTGCGCCTGTTTCCGGATCACTGGCGCTTAGGTGCGCCGTCGAAGCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGCACCGCTGCGTATTCAGCGTTAGGGGCAGGAGAAGCAATATCTTGTGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAACGAAGGCAGGGGGATCGGTGGCGATGACGATGCTGGGAAGCGATGCTGATTTGGGTGTCGGGCGGTTTAGCTGTcgcgt                                                                                                         >  2:224533/5‑250 (MQ=255)
                                                                                                                                           cGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGAt                                                                                                        <  1:400148/251‑1 (MQ=255)
                                                                                                                                              gTCTGCGACTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGGGGGGCTACCGGAAGCGAAAGCGGTGGGATCGGGGGGGATGTCGGTGCGTGGAAGCGATGATGATTTGGGGGTCggtatgtttcccagtgtcggtgcc                                                                                                     >  2:257835/1‑228 (MQ=255)
                                                                                                                                                 tGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGt                                                                                                  >  1:510018/1‑251 (MQ=255)
                                                                                                                                                 ggcgcCTGTTTCCGGATCGCTGGCGCTTAGTTGCGCCGTCGAACCACACTCATCTGGCAGTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTTCCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGTTGATTTGGGTGTCGGGCTGTTTACCAGTCGCGATGCCAGt                                                                                                  <  2:393758/250‑1 (MQ=255)
                                                                                                                                                  gcgcttggttCCGGATCGCTGGCGCTTAGGTGCGCCGGCGATCCACACACATCTGGCATTAAGCCGTCAGTCTTTCTAACCGCTGCGTATTCAGCGGTAGGGGCAGGAACAGCACTATCGGGGGACGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGGTCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGGGTCGTGCTGTTTACCAGTCGCGATGCCAg                                                                                                   <  2:418999/241‑1 (MQ=255)
                                                                                                                                                  gcgcCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTGCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAg                                                                                                   <  2:177096/249‑1 (MQ=255)
                                                                                                                                                  gcgcATGTTGCCGGATCGCTGTCGCTTAGGTGCGCCGTCGAACCACGCTCACCTGTCATTAAGCCGTCAGTCGTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGATAGCGGTGGGATCGTTGGCGATGTCGATTCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTTTTTACCAGTCGCGATGCCAGt                                                                                                  <  2:111187/250‑1 (MQ=255)
                                                                                                                                                   cgcCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCa                                                                                                <  1:620743/251‑1 (MQ=255)
                                                                                                                                                   cgcCTGTTTCCGGATCGCTGGCGCGTAGGTGCGCCGTCGAACCACACTCATCTGGCATGAAGCCGTCAGTCGTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCa                                                                                                <  2:38254/251‑1 (MQ=255)
                                                                                                                                                        gTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTAt                                                                                            <  1:595371/250‑1 (MQ=255)
                                                                                                                                                              ggATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAg                                                                                     <  1:451796/251‑1 (MQ=255)
                                                                                                                                                                       gcgcTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGaa                                                                            <  1:212471/251‑1 (MQ=255)
                                                                                                                                                                        cgcTTAGGTGCGCCGGCGAACCGCACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGAGTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAAc                                                                           <  1:354071/251‑1 (MQ=255)
                                                                                                                                                                           ttAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGc                                                                         >  1:597719/1‑250 (MQ=255)
                                                                                                                                                                                tGCGACGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGCAAACGGTGGGATCGGTGGCGATGTCGCTGCTGGGAAGCAATGCTGATTTGGGTGTCATGCTGTTTCCCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TAc                                                                  >  1:445084/1‑251 (MQ=255)
                                                                                                                                                                                    ccGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCa                                                              <  1:252562/251‑1 (MQ=255)
                                                                                                                                                                                      gTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGTGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGG‑TACTTCATg                                                            >  1:200949/1‑251 (MQ=255)
                                                                                                                                                                                        cGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGa                                                           >  1:242580/1‑250 (MQ=255)
                                                                                                                                                                                               cacTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTgcgc                                                   >  1:365200/1‑251 (MQ=255)
                                                                                                                                                                                                    agcTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATGGCGCTGATg                                              <  2:642392/249‑1 (MQ=255)
                                                                                                                                                                                                     tCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCACTGATGt                                             >  1:536213/1‑251 (MQ=255)
                                                                                                                                                                                                             tAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCgg                                       >  1:333941/1‑249 (MQ=255)
                                                                                                                                                                                                                gCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTATGGGCAGGAACAGCCCTACCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTAACGGAAGCGAACGCGGTGGGATCGGGGGCGATGTCGATGCTGGGAAGCGATGCTGATTTTGGGGTCGTGCTGGTTACCCATCGCGGTGCCAGTCACTCTCAACAAGGGCCAGGAAAGCCGTATA‑TTCATGAAATTGCGCAGATGGTGCCGGAggtt                                  >  2:61048/1‑248 (MQ=255)
                                                                                                                                                                                                                gCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCtt                                  <  1:120987/251‑1 (MQ=255)
                                                                                                                                                                                                                  cgtccgtcTTTCGACCTGCTGCGTATTCAGCGTTAGGGGCAGTAACAGCACTTTCGTGGGCCGCTTACCGGGCTCGCGGTGCTGGTGGTTCTACCGGAATCGAAAGCGGTGGGATCGGTGGCGATTTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGGTTAACAGTCGCGATGACAGTCTCTATCAACAAGGGCAAGGAACGCAGT‑TACTACATTAAATTGCGCTGATTTTCCCGGACCtttt                                <  2:151888/251‑3 (MQ=255)
                                                                                                                                                                                                                     cAGTCTTTCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAATCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGa                             <  1:96219/251‑1 (MQ=255)
                                                                                                                                                                                                                            tCGAACCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCGTGGGCCGCTTTACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTg                       <  1:272903/250‑1 (MQ=255)
                                                                                                                                                                                                                               aaCCGCTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGAtt                   <  1:347633/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGCGGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGa                 <  2:365200/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGCGTATTCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGGGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAAcg              >  2:663802/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gcgtgtgcatcGTTAGGGGCAGGATCAGCTCTATCTTGGGCCTCTTAACTGTCCTGAGCTGCTGGTGGTGCTACCGGAAGCGCAAGCGGTGGGATTGGTGGCGATGTCGATGCTGGGAAGCGTTGCTGATTTGGGTGTCGTTCTGTTTACCAGTCGCGATGCCAGTCACTATCAACCAGGGCAAGGAACGCAGT‑GTCTTCATGTAATTGCGCTGATGTTGCCGGAGCTTCTGGAGTGTTGGATTGAAcgc             <  2:545522/240‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATg        >  1:706503/1‑250 (MQ=255)
                                                                                                                                                                                                                                           ttCAGCGATAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGa                             <  1:697964/229‑1 (MQ=255)
                                                                                                                                                                                                                                             cAGCGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCGTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGAcc     <  2:279789/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGACCGAt  >  1:147239/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cGTTAGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGACCGAt  >  1:105383/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
TTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACCGCCGCCAGCAGTCTCGACGATATGCTGATGCGCTTTCACCGCTGGGCGCGCGATCTCGGCCTGGCAGGTGCGAGTCTGCGCCTGTTTCCGGATCGCTGGCGCTTAGGTGCGCCGTCGAACCACACTCATCTGGCATTAAGCCGTCAGTCTTTCGAACCGCTGCGTATTCAGCGTTTGGGGCAGGAACAGCACTATCTTGGGCCGCTTAACGGACCAGAGCTGCTGGTGGTGCTACCGGAAGCGAAAGCGGTGGGATCGGTGGCGATGTCGATGCTGGGAAGCGATGCTGATTTGGGTGTCGTGCTGTTTACCAGTCGCGATGCCAGTCACTATCAACAAGGGCAAGGAACGCAGT‑TACTTCATGAAATTGCGCTGATGTTGCCGGAGCTTCTGGAGCGTTGGATTGAACGCGTATGACCGAT  >  NC_000913/3995805‑3996295

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: ‑