Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP034943 607,838 A→G 5.3% T227T (ACA→ACG EO946_RS03185 → collagen‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP034943607,8380AG 5.3% 1778.7 / 49.3 544T227T (ACA→ACGEO946_RS03185collagen‑like protein
Reads supporting (aligned to +/- strand):  ref base A (270/236);  new base G (13/16);  total (289/255)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 562 total aligned reads displayed.

AAGCCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCACCCCAGGAATTACAC  >  NZ_CP034943/607699‑607977
                                                                                                                                           |                                                                                                                                           
acgtcGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACa                                                                                                                                             >  1:1043913/5‑140 (MQ=255)
 cgtcGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCCGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGCAGCGACAg                                                                                                                                            >  2:670010/4‑140 (MQ=255)
  gtcGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAgg                                                                                                                                           >  1:5647814/3‑140 (MQ=255)
  gtcGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGCGAAGCCGGACCACCCGGAGCGCCAgg                                                                                                                                           >  2:2854500/3‑140 (MQ=255)
  gCCGGACCAGCCGGAGCCACAGGAGCTACCGGGGAAGCCGGACCACCCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACggg                                                                                                                                           <  1:797095/140‑1 (MQ=255)
    cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGACGCCGGACCAACCGGAGCGACAGGAg                                                                                                                                         >  2:1396717/1‑140 (MQ=255)
    cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCACCCGGAGCGCCGGGAg                                                                                                                                         >  2:3546296/1‑140 (MQ=255)
    cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAg                                                                                                                                         >  1:3703525/1‑140 (MQ=255)
    cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAg                                                                                                                                         >  1:1792523/1‑140 (MQ=255)
    cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAg                                                                                                                                         >  2:4446556/1‑140 (MQ=255)
    cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAg                                                                                                                                         >  2:4636646/1‑140 (MQ=255)
    cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCACCCGGAGCGACAGGAGCTACCGAGGAAGCCGGACCAACCGCAGCGAcagcag                                                                                                                                         >  2:5272621/1‑140 (MQ=255)
            cggtaGCGAGTGGAGCTACGGTGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGa                                                                                                                                 <  1:3352476/136‑1 (MQ=17)
            ccGGAGCGACAGGAGCTACCGGCGAAGCCGGACCAACCGGCGCGACGGGAGCCACAGGGGCTACTGGGGCCGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGa                                                                                                                                 <  1:3914623/140‑1 (MQ=255)
             cGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACGGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAg                                                                                                                                <  1:4450349/140‑1 (MQ=255)
                 gCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcac                                                                                                                            >  1:2406232/1‑140 (MQ=255)
                     cAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGCCGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCgactggg                                                                                                                        >  2:3821726/1‑139 (MQ=14)
                      aGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcacaggca                                                                                                                       <  2:5454463/140‑1 (MQ=255)
                      aGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcacaggca                                                                                                                       >  1:1305543/1‑140 (MQ=255)
                          tctcccGGCGAAGCCTGACCAACCGCTGCGACGTGAGCGCCAGGGTCTCCTGGGGACCTGGGACCCGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcacaggcacagg                                                                                                                   <  2:2570455/136‑1 (MQ=18)
                               cGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcaca                                                                                                                           <  1:5672548/127‑1 (MQ=255)
                               cGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcaca                                                                                                                           >  2:5672548/1‑127 (MQ=255)
                                  ggAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACCGGCACAGGAGCAGCTg                                                                                                           >  1:5382316/1‑140 (MQ=255)
                                   gAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGCGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTgg                                                                                                          >  2:882492/1‑140 (MQ=255)
                                   gAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCCCGGGAGCCACAGGCAGAGGAGCAGCTgg                                                                                                          >  2:3425014/1‑140 (MQ=255)
                                   gAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTgg                                                                                                          >  1:4631968/1‑140 (MQ=255)
                                        cgccccACCCGGAGCGCCGGGAGCCACCGGTGCTACCGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAt                                                                                                     <  1:1417613/136‑1 (MQ=21)
                                        cGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAt                                                                                                     <  1:439094/140‑1 (MQ=255)
                                        cGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAt                                                                                                     >  2:293247/1‑140 (MQ=255)
                                        cGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAt                                                                                                     >  2:3918109/1‑140 (MQ=255)
                                           aCCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCCCGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCg                                                                                                  >  2:1664365/1‑140 (MQ=255)
                                           aCCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCg                                                                                                  >  2:591900/1‑140 (MQ=255)
                                               aCCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATa                                                                                              <  1:2668669/140‑1 (MQ=255)
                                                    cgccaCGGGAGCCACAGGGGCTACTGGGGAAGTCGTACCAGCCGGAGCGACAGGAGCTACCGGCGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         <  2:2189976/136‑1 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGCGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         >  1:2928653/1‑140 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         >  2:322207/1‑140 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         >  2:4363670/1‑140 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         <  1:2973624/140‑1 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         >  2:493589/1‑140 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         <  1:1715581/140‑1 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         <  1:591900/140‑1 (MQ=255)
                                                    aGCGACGGGAGCCACAGGGGCTACTGGGGACGACGGACCAGCCGGAGCGACAGGCGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGCAGCCACGGGAGCCACAGCCACAGCAGCAGCGGGGTTATCCGAATACgggg                                                                                         >  2:4768111/1‑139 (MQ=21)
                                                        aCGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcac                                                                                                                            <  1:4120701/101‑1 (MQ=255)
                                                                   aGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGCCAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTAGATATCTAACTTaa                                                                          >  1:5080151/1‑140 (MQ=255)
                                                                   aGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTaa                                                                          >  1:1146844/1‑140 (MQ=255)
                                                                   aGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTaa                                                                          >  1:2276123/1‑140 (MQ=255)
                                                                         tACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         >  1:4912257/1‑119 (MQ=255)
                                                                          accggggACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCa                                                                   <  1:3543485/137‑1 (MQ=255)
                                                                          aCTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATCTAACTTAAGCGCTCa                                                                   >  2:1807703/1‑140 (MQ=255)
                                                                          aCTGGGGACGGCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCa                                                                   <  1:735785/140‑1 (MQ=255)
                                                                           ccggggAAGCCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGCAGCAGCTGGGGTATCCGAAGACGGGGAGCGATCTAACTTATGCGCTcga                                                                  >  1:2823422/3‑138 (MQ=18)
                                                                             ggggACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCACCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAAc                                                                <  2:4124153/140‑1 (MQ=255)
                                                                               ggACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACgg                                                              >  1:4685333/1‑140 (MQ=255)
                                                                                gACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGt                                                             <  1:2008334/140‑1 (MQ=255)
                                                                                gACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGt                                                             <  1:3119517/140‑1 (MQ=255)
                                                                                gACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGt                                                             <  1:4976359/140‑1 (MQ=255)
                                                                                gACGTCGGACCACCTGGAGCGACATGAGCTACCGGGGAAGCCGCACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGt                                                             <  1:3414776/140‑1 (MQ=255)
                                                                                 aaggcGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCACCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTa                                                            <  1:898983/136‑1 (MQ=255)
                                                                                 aCGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAt                                                                                                     <  2:3677592/99‑1 (MQ=255)
                                                                                   gccGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGcc                                                                                                                              <  2:2021929/70‑1 (MQ=255)
                                                                                   gTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACCGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGc                                                          >  2:4690836/1‑140 (MQ=255)
                                                                                   gTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGc                                                          >  2:2251628/1‑140 (MQ=255)
                                                                                     cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCg                                                                                                  <  2:92588/98‑1 (MQ=255)
                                                                                     cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATa                                                                                              <  1:169411/102‑1 (MQ=255)
                                                                                     cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         >  2:5590294/1‑107 (MQ=255)
                                                                                     cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCta                                                        >  1:5291871/1‑140 (MQ=255)
                                                                                     cGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCta                                                        >  1:3280848/1‑140 (MQ=255)
                                                                                           aGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAgg                                                  >  2:4078684/1‑140 (MQ=255)
                                                                                           aGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAgg                                                  >  2:3993855/1‑140 (MQ=255)
                                                                                             tgGGAGCCACAGGAGCTACCGGGGAAGCCGGACCAGCCGGGGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCt                                                <  2:113054‑M2/138‑1 (MQ=255)
                                                                                             cgGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAGCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAt                                                                                                     >  2:1229440‑M2/3‑87 (MQ=255)
                                                                                             cgGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAGCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAt                                                                                                     <  1:1229440‑M2/85‑1 (MQ=255)
                                                                                             ccGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGt                                                                                         <  2:1505019/99‑1 (MQ=255)
                                                                                                 aGCGACAGGAGCTACCGGGGAAGCCGGACCAGCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATa                                                                                              >  2:1924111/1‑90 (MQ=255)
                                                                                                 aGCGACAGGAGCTACCGGGGAAGCCGGACCACCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATa                                                                                              <  1:1924111/90‑1 (MQ=255)
                                                                                                  gCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGt                                           >  2:5360417/1‑140 (MQ=255)
                                                                                                       aGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCAttt                                      >  1:87641/1‑140 (MQ=255)
                                                                                                       aGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCAttt                                      >  1:363722/1‑140 (MQ=255)
                                                                                                           gCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGa                                  <  1:993764/140‑1 (MQ=255)
                                                                                                              aCCGGGGAAGCCGGACCAGCCGGGGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTggg                                                                                                         <  1:2541707/66‑1 (MQ=255)
                                                                                                              aCCGGGGAAGCCGGACCACCCGGAGCGACAGGAGCCACGGGAGCcacaggcacagg                                                                                                                   <  1:2141162/56‑1 (MQ=255)
                                                                                                                cGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGtatata                                                                                    <  1:4026217/85‑1 (MQ=255)
                                                                                                                cGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCaa                             >  1:5234657/1‑140 (MQ=255)
                                                                                                                cGGGGAAGCCGGAAGACCCGGAGCGACTGTAGTCACGGGAGCCACAGGCACAGGAGCAGATGGGTTATCCGAATACGTGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCaa                             <  2:2531304/140‑1 (MQ=255)
                                                                                                                    gAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCcaca                                                                                                                           <  2:2351909/42‑1 (MQ=255)
                                                                                                                    gAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATa                                                                                              >  2:4203286/1‑71 (MQ=255)
                                                                                                                    gAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATa                                                                                              <  1:4203286/71‑1 (MQ=255)
                                                                                                                       gCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGAt                      >  2:30912/1‑140 (MQ=255)
                                                                                                                        ccGGACCAACCGGAGCGACAGGAGCCACGGGAGCcac                                                                                                                            <  1:5491698/37‑1 (MQ=255)
                                                                                                                            aCCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAacca                 >  1:4771336/1‑140 (MQ=255)
                                                                                                                                aCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCAcccc             <  2:5061659/140‑1 (MQ=255)
                                                                                                                                aCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCAcccc             <  1:4251833/140‑1 (MQ=255)
                                                                                                                                aCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCAcccc             <  2:764321/140‑1 (MQ=255)
                                                                                                                                aCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCAcccc             <  1:5181513/140‑1 (MQ=255)
                                                                                                                                aCCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTAACCGAATACGGGGATATATATAACTTAAGCCATCAAACGGTAGATATAGAGGCTGATGTAATTTTTGATTCAACAGGGATAACCAccca             >  2:4831348/1‑139 (MQ=255)
                                                                                                                                 ccGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCAcccc             >  1:3841826/1‑139 (MQ=255)
                                                                                                                                 ccGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCACCCCa            >  1:3067841/1‑140 (MQ=255)
                                                                                                                                     acccccAGGCGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCACCCCAGGaa        <  2:2712369/135‑1 (MQ=255)
                                                                                                                                      gCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCACCCCAGGAAt       >  1:1631692/1‑140 (MQ=255)
                                                                                                                                      gCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCACCCCAGGAAt       <  1:1587787/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AAGCCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACGGGAGCCACAGGGGCTACTGGGGACGTCGGACCAGCCGGAGCGACAGGAGCTACCGGGGAAGCCGGACCAACCGGAGCGACAGGAGCCACGGGAGCCACAGGCACAGGAGCAGCTGGGTTATCCGAATACGGGTATATATATAACTTAAGCGCTCAAACGGTAGCTATAGAGGCTGATGTCATTTTTGATTCAACAGGGATAACCACCCCAGGAATTACAC  >  NZ_CP034943/607699‑607977

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