Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,365,253 G→T 34.8% I227I (ATC→ATA)  nanT ← sialic acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,365,2530GT34.8% 541.8 / 302.2 316I227I (ATC→ATA) nanTsialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (75/131);  new base T (51/59);  total (126/190)
Fisher's exact test for biased strand distribution p-value = 9.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 321 total aligned reads displayed.

CCCGTTGGCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAA  >  CP009273/3365104‑3365402
                                                                                                                                                     |                                                                                                                                                     
  cGTTGGCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGat                                                                                                                                                     <  1:3207426/150‑1 (MQ=255)
   gTTGGCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGata                                                                                                                                                    <  1:1335087/150‑1 (MQ=255)
   gTTGGCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGata                                                                                                                                                    <  1:4106091/150‑1 (MQ=255)
   gTTGGCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGata                                                                                                                                                    <  1:844156/150‑1 (MQ=255)
      ggCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTg                                                                                                                                                 >  2:3642950/1‑150 (MQ=255)
      ggCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTg                                                                                                                                                 <  1:2477643/150‑1 (MQ=255)
        ccAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGc                                                                                                                                               <  1:3389289/150‑1 (MQ=255)
            cGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATg                                                                                                                                           <  1:2396263/150‑1 (MQ=255)
             gTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGc                                                                                                                                          <  2:2849621/150‑1 (MQ=255)
              tttGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCg                                                                                                                                         <  2:2761784/150‑1 (MQ=255)
                  ccTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGt                                                                                                                                     <  2:854362/150‑1 (MQ=255)
                   cTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGtt                                                                                                                                    >  1:2148456/1‑150 (MQ=255)
                    tGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTc                                                                                                                                   <  2:3467333/150‑1 (MQ=255)
                        ctctGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTcacca                                                                                                                               <  1:3629211/150‑1 (MQ=255)
                         tctGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTcaccac                                                                                                                              <  1:151015/150‑1 (MQ=255)
                            gCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGt                                                                                                                           <  2:3455306/150‑1 (MQ=255)
                            gCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGt                                                                                                                           <  2:3369676/150‑1 (MQ=255)
                             cACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTa                                                                                                                          <  2:3584678/150‑1 (MQ=255)
                              aCCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTag                                                                                                                         >  1:4502429/1‑150 (MQ=255)
                               ccATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTaga                                                                                                                        <  2:2640470/150‑1 (MQ=255)
                                  tttcgCTGATCACGATTTCGGAGAATCACAGGCCAAGAACCGCGACGATCGCGGCATTTTGCTGGTTTCCGGCGAAGCTGAACATCAGCGCAGTTGCGTGCGCCTTTGTCATTACTTTATTGGCTATGCGATGTTCACCTCTGTagagtt                                                                                                                     <  2:1304776/146‑3 (MQ=255)
                                   aaaGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAat                                                                                                                    >  2:489716/1‑150 (MQ=255)
                                    aaGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAata                                                                                                                   <  2:2864290/150‑1 (MQ=255)
                                    aaGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAata                                                                                                                   <  1:1317748/150‑1 (MQ=255)
                                     aGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAatat                                                                                                                  >  2:2617426/1‑150 (MQ=255)
                                       cTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATcc                                                                                                                <  2:1602181/150‑1 (MQ=255)
                                       cTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATcc                                                                                                                >  2:2187258/1‑150 (MQ=255)
                                             aaGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCAtt                                                                                                          >  1:364424/1‑150 (MQ=255)
                                                aTTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATtgtg                                                                                                       <  2:1295159/150‑1 (MQ=255)
                                                  tGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCg                                                                                                     <  1:4091817/150‑1 (MQ=255)
                                                  tGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCg                                                                                                     <  1:78428/150‑1 (MQ=255)
                                                    cggcgCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTa                                                                                                   <  2:2940152/150‑1 (MQ=255)
                                                     ggcgCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTAc                                                                                                  >  1:3134053/1‑150 (MQ=255)
                                                         cATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGt                                                                                              >  1:2970504/1‑150 (MQ=255)
                                                          aTAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTg                                                                                             >  2:3572394/1‑150 (MQ=255)
                                                          aTAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTg                                                                                             >  1:1956915/1‑150 (MQ=255)
                                                           tAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGc                                                                                            <  2:2732810/150‑1 (MQ=255)
                                                             aCAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCtt                                                                                          >  1:2078919/1‑150 (MQ=255)
                                                                 cccAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTAcc                                                                                      <  1:3841957/150‑1 (MQ=255)
                                                                  ccAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCt                                                                                     <  1:2882850/150‑1 (MQ=255)
                                                                  ccAAGAACAGCGACGATCGCGGCATTTGGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCt                                                                                     <  2:2816420/150‑1 (MQ=255)
                                                                    aaGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGc                                                                                   >  1:2286464/1‑150 (MQ=255)
                                                                      gAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCgt                                                                                 >  1:3239924/1‑150 (MQ=255)
                                                                          aGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGttt                                                                             <  2:1911983/150‑1 (MQ=255)
                                                                                   cgcgGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCa                                                                    <  1:676551/150‑1 (MQ=255)
                                                                                   cgcgGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCa                                                                    >  2:1614943/1‑150 (MQ=255)
                                                                                   cgcgGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCa                                                                    <  2:4572698/150‑1 (MQ=255)
                                                                                    gcgGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAg                                                                   <  2:4260192/150‑1 (MQ=255)
                                                                                    gcgGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAg                                                                   <  1:4242968/150‑1 (MQ=255)
                                                                                    gcgGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAg                                                                   >  2:3901730/1‑150 (MQ=255)
                                                                                     cgGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGt                                                                  <  1:100113/150‑1 (MQ=255)
                                                                                      ggCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTc                                                                 <  2:3214873/150‑1 (MQ=255)
                                                                                       gCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTct                                                                >  2:1245852/1‑150 (MQ=255)
                                                                                       gCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTct                                                                <  1:1578523/150‑1 (MQ=255)
                                                                                        cATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTctt                                                               <  2:3385103/150‑1 (MQ=255)
                                                                                        cATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTctt                                                               <  1:3729780/150‑1 (MQ=255)
                                                                                         aTTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTcttc                                                              <  2:4320119/150‑1 (MQ=255)
                                                                                             tGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTcttccgct                                                          <  1:2789615/150‑1 (MQ=255)
                                                                                                 ggTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTcttccgcttccg                                                      >  2:3105868/1‑150 (MQ=255)
                                                                                                 ggTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTcttccgcttccg                                                      <  2:2505159/150‑1 (MQ=255)
                                                                                                   ttACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCggg                                                    >  1:3997176/1‑150 (MQ=255)
                                                                                                    tACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGa                                                   >  1:894809/1‑150 (MQ=255)
                                                                                                    tACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGa                                                   <  2:4042709/150‑1 (MQ=255)
                                                                                                      ccGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATg                                                 <  1:3933014/150‑1 (MQ=255)
                                                                                                           gAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGttttt                                            >  2:3716386/1‑150 (MQ=255)
                                                                                                               cAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGc                                        >  2:3913734/1‑150 (MQ=255)
                                                                                                                aGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCa                                       >  1:23554/1‑150 (MQ=255)
                                                                                                                aGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCa                                       <  1:149940/150‑1 (MQ=255)
                                                                                                                aGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCa                                       <  2:238472/150‑1 (MQ=255)
                                                                                                                aGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCa                                       <  2:3167109/150‑1 (MQ=255)
                                                                                                                 gAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAg                                      <  1:4054152/150‑1 (MQ=255)
                                                                                                                 gAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAg                                      >  1:4536811/1‑150 (MQ=255)
                                                                                                                  aaCCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGc                                     >  1:1069912/1‑150 (MQ=255)
                                                                                                                    ccACAGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCa                                   <  2:864954/150‑1 (MQ=255)
                                                                                                                    ccACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCa                                   >  1:1855879/1‑150 (MQ=255)
                                                                                                                    ccACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCa                                   <  1:2003899/150‑1 (MQ=255)
                                                                                                                      acaGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCaga                                 >  2:2407792/1‑150 (MQ=255)
                                                                                                                      acaGCGCAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCaga                                 >  2:1806534/1‑150 (MQ=255)
                                                                                                                      acaGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCaga                                 >  2:331854/1‑150 (MQ=255)
                                                                                                                       caGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCagag                                <  2:4504799/150‑1 (MQ=255)
                                                                                                                         gcgcAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCagagag                              <  2:817642/150‑1 (MQ=255)
                                                                                                                          cgcAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGc                             <  1:124328/150‑1 (MQ=255)
                                                                                                                          cgcAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGc                             <  1:3600/150‑1 (MQ=255)
                                                                                                                          cgcAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGc                             >  2:1094411/1‑150 (MQ=255)
                                                                                                                            cAGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCaa                           <  2:353316/150‑1 (MQ=255)
                                                                                                                             aGTAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCaaa                          >  1:43925/1‑150 (MQ=255)
                                                                                                                             aGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCaaa                          <  2:4083537/150‑1 (MQ=255)
                                                                                                                             aGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCaaa                          <  2:275482/150‑1 (MQ=255)
                                                                                                                              gTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAg                         <  1:489716/150‑1 (MQ=255)
                                                                                                                               tAGCCGCCGCCAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAga                        >  2:709346/1‑150 (MQ=255)
                                                                                                                                 gccgccgccAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAgatg                      <  1:1456800/150‑1 (MQ=255)
                                                                                                                                   cgccgccAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAgatgat                    >  2:1541267/1‑150 (MQ=255)
                                                                                                                                     ccgccAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTg                  <  2:4373011/150‑1 (MQ=255)
                                                                                                                                      cgccAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTgg                 <  2:4357317/150‑1 (MQ=255)
                                                                                                                                       gccAGTGTCATTACTATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGc                <  1:1742827/150‑1 (MQ=255)
                                                                                                                                         cAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCaa              >  1:2322831/1‑150 (MQ=255)
                                                                                                                                          aGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCaaa             >  1:2615134/1‑150 (MQ=255)
                                                                                                                                             gtCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATg          <  2:1379054/150‑1 (MQ=255)
                                                                                                                                             gtCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATg          <  2:3533725/150‑1 (MQ=255)
                                                                                                                                                  tACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGAt     <  1:951804/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCGTTGGCCAGCGTTTGCCTGCACTCTGCACCATAAAGCTGATAAAGATTGCGGCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTCATTACGATATTGGCAATGCGATGTTCACCACGGTAGAGAATATCCACCATTGTGCGTACTGGTGCTTTACCTGCGTGTTTCTCTTTCCAGTCTTCCGCTTCCGGGATGTTTTTACGCAGCCAGAGAGCAAAGATGATTGGCAAAATGCCGATAAA  >  CP009273/3365104‑3365402

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