Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,576,325 C→A 94.3% E171D (GAG→GAT)  nudK ← GDP‑mannose pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,576,3250CA94.3% 645.5 / 25.8 193E171D (GAG→GAT) nudKGDP‑mannose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (8/3);  new base A (102/80);  total (110/83)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 202 total aligned reads displayed.

GTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCC  >  CP009273/2576178‑2576471
                                                                                                                                                   |                                                                                                                                                  
gTATCTTCGTTGCGCAGGTAGCGCGGCTCACTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtatc                                                                                                                                                  >  2:1529557/1‑150 (MQ=255)
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtctc                                                                                                                                                  >  2:839581/1‑150 (MQ=255)
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtctc                                                                                                                                                  >  2:4570316/1‑150 (MQ=255)
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtatc                                                                                                                                                  >  2:1046875/1‑150 (MQ=255)
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtatc                                                                                                                                                  >  2:1258838/1‑150 (MQ=255)
 tATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtatca                                                                                                                                                 <  2:2933704/150‑2 (MQ=255)
     ttCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCgg                                                                                                                                             >  1:877094/1‑150 (MQ=255)
        gTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACCACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGttt                                                                                                                                          >  2:868561/1‑150 (MQ=255)
         ttGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGtttt                                                                                                                                         <  2:3870253/150‑1 (MQ=255)
          tGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTg                                                                                                                                        <  1:1602308/150‑1 (MQ=255)
            cgcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGat                                                                                                                                      <  2:1787289/150‑1 (MQ=255)
            cgcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGat                                                                                                                                      >  2:4470459/1‑150 (MQ=255)
             gcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGatc                                                                                                                                     <  1:3253827/150‑1 (MQ=255)
               aGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatcat                                                                                                                                   >  1:2956096/1‑150 (MQ=255)
                      gcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAgc                                                                                                                            >  2:593520/1‑150 (MQ=255)
                          cTCCATCGGATTTGTCGGATTGAGCCAATAATACACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCt                                                                                                                        <  2:2610044/150‑1 (MQ=255)
                          cTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCt                                                                                                                        <  2:34740/150‑1 (MQ=255)
                          cTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCt                                                                                                                        <  2:3320663/150‑1 (MQ=255)
                              aTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCt                                                                                                                    >  2:18277/1‑150 (MQ=255)
                               tCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTg                                                                                                                   <  2:716506/150‑1 (MQ=255)
                                 ggATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGaa                                                                                                                 >  2:2292611/1‑150 (MQ=255)
                                     ttGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATAATCTCCAGCGCCTGGCTGAACGGc                                                                                                             >  2:4164642/1‑150 (MQ=255)
                                          ggATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTc                                                                                                        <  1:2700536/150‑1 (MQ=255)
                                          ggATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTc                                                                                                        >  1:3065767/1‑150 (MQ=255)
                                             ttGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAg                                                                                                     <  2:1312126/150‑1 (MQ=255)
                                              tGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTTATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGc                                                                                                    >  1:3318140/1‑150 (MQ=255)
                                              tGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGc                                                                                                    >  1:2681548/1‑150 (MQ=255)
                                              tGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCGCGAGc                                                                                                    >  1:2680551/1‑150 (MQ=255)
                                                aGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCAc                                                                                                  >  1:1411393/1‑150 (MQ=255)
                                                aGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCAc                                                                                                  >  1:3804127/1‑150 (MQ=255)
                                                         taaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatc                                                                                         >  2:2156359/1‑150 (MQ=255)
                                                          aaaCCACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatct                                                                                        <  1:1906905/150‑1 (MQ=255)
                                                          aaaCAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatct                                                                                        >  2:1340578/1‑150 (MQ=255)
                                                           aacaacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatctt                                                                                       >  1:1147377/1‑150 (MQ=255)
                                                            acaacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttc                                                                                      >  1:4414000/1‑150 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCCGCTCGAGCACTTCAATatcttcat                                                                                    >  1:784690/1‑150 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                    <  2:1398553/150‑1 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                    <  2:1398603/150‑1 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                    >  1:1035304/1‑150 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                    >  1:783338/1‑150 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                    >  1:1034636/1‑150 (MQ=255)
                                                                cGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatct                                                                                  >  2:4409329/1‑150 (MQ=255)
                                                                 gAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatctt                                                                                 >  1:1716940/1‑150 (MQ=255)
                                                                   aaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCg                                                                               >  1:637113/1‑150 (MQ=255)
                                                                   aaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCg                                                                               <  1:2280309/150‑1 (MQ=255)
                                                                    atatTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTTCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATATCAAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTAGa                                                                              >  2:1003468/1‑150 (MQ=255)
                                                                    atatTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTAATCTTCGa                                                                              >  1:2190660/1‑150 (MQ=255)
                                                                     tattatTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGac                                                                             >  2:3639977/1‑150 (MQ=255)
                                                                      attattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGaca                                                                            <  1:4692614/150‑1 (MQ=255)
                                                                        tattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACAcc                                                                          >  1:400325/1‑150 (MQ=255)
                                                                        tattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACAcc                                                                          >  2:1279852/1‑150 (MQ=255)
                                                                           tGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGcc                                                                       <  1:2876911/150‑1 (MQ=255)
                                                                            gTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCa                                                                      >  2:3883324/1‑150 (MQ=255)
                                                                             ttATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAc                                                                     >  1:857046/1‑150 (MQ=255)
                                                                               aTTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCg                                                                   <  2:1822423/150‑1 (MQ=255)
                                                                                     cAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTa                                                             >  1:1513223/1‑150 (MQ=255)
                                                                                        tCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcg                                                          >  1:3700669/1‑150 (MQ=255)
                                                                                         ccATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgc                                                         >  2:4601160/1‑150 (MQ=255)
                                                                                             tAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTg                                                     >  2:2154775/1‑150 (MQ=255)
                                                                                              aaaTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGa                                                    >  2:2467663/1‑150 (MQ=255)
                                                                                                 tgtgACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTg                                                 >  1:3384441/1‑150 (MQ=255)
                                                                                                  gtgACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGt                                                >  2:819876/1‑150 (MQ=255)
                                                                                                    gACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCa                                              >  2:3045179/1‑150 (MQ=255)
                                                                                                    gACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCa                                              <  1:2161394/150‑1 (MQ=255)
                                                                                                       gTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTg                                           <  2:255866/150‑1 (MQ=255)
                                                                                                           gCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTAtt                                       <  2:427647/150‑1 (MQ=255)
                                                                                                            cAAATAGTTAAGCAATAACACCGTCTGACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTc                                      <  1:1860206/150‑1 (MQ=255)
                                                                                                             aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCCGCTCGCGCCCGTCAATATCTTCAACTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCCCTGTATTCg                                     >  1:2385098/1‑150 (MQ=255)
                                                                                                             aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                     <  1:4523831/150‑1 (MQ=255)
                                                                                                              aaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCgg                                    >  1:734642/1‑150 (MQ=255)
                                                                                                                  gTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGAt                                >  1:3116383/1‑150 (MQ=255)
                                                                                                                   ttAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATa                               <  1:4127317/150‑1 (MQ=255)
                                                                                                                   ttAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATa                               <  1:1741983/150‑1 (MQ=255)
                                                                                                                      aGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaa                            >  2:687803/1‑150 (MQ=255)
                                                                                                                        cAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaaaa                          >  2:3346168/1‑150 (MQ=255)
                                                                                                                          aTAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGt                        <  2:400325/150‑1 (MQ=255)
                                                                                                                           tAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTg                       >  1:2225799/1‑150 (MQ=255)
                                                                                                                           tAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTg                       <  1:4316673/150‑1 (MQ=255)
                                                                                                                           tAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTAATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTAGGCGCTaaaaaattt                       >  1:3585960/1‑147 (MQ=255)
                                                                                                                            aaCACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCCGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCCCCGGCGTTAGCGCGCTGATTGTCACGGTATTCGGCGATAAAAAAGTgg                      >  2:2034162/1‑150 (MQ=255)
                                                                                                                              cacCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAt                    <  2:264814/150‑1 (MQ=255)
                                                                                                                                 cGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAg                 >  1:2282737/1‑150 (MQ=255)
                                                                                                                                 cGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGGTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAg                 >  2:619607/1‑150 (MQ=255)
                                                                                                                                    cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc              <  1:1927936/150‑1 (MQ=255)
                                                                                                                                    cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc              <  1:165949/150‑1 (MQ=255)
                                                                                                                                    cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCGGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc              <  1:2939197/150‑1 (MQ=255)
                                                                                                                                      tACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAg            >  2:2404767/1‑150 (MQ=255)
                                                                                                                                        ccGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTc          <  1:1046875/150‑1 (MQ=255)
                                                                                                                                        ccGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTc          <  2:333772/150‑1 (MQ=255)
                                                                                                                                        ccGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTc          >  2:3329609/1‑150 (MQ=255)
                                                                                                                                        ccGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTc          <  1:1634415/150‑1 (MQ=255)
                                                                                                                                         cGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTca         >  1:2596442/1‑150 (MQ=255)
                                                                                                                                          gTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcac        <  1:1472952/150‑1 (MQ=255)
                                                                                                                                           tCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcaca       >  2:3832986/1‑150 (MQ=255)
                                                                                                                                           tCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcaca       >  2:3508498/1‑150 (MQ=255)
                                                                                                                                            cACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCacac      >  1:4271350/1‑150 (MQ=255)
                                                                                                                                              cGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACAGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCg    >  1:2130092/1‑150 (MQ=255)
                                                                                                                                               gTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGc   >  2:2745720/1‑150 (MQ=255)
                                                                                                                                               gTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACAGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGc   >  2:3610564/1‑150 (MQ=255)
                                                                                                                                                tATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGAGAAAAAAGTGGATTAGCTCAGTCACACCGcc  >  2:3558167/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCC  >  CP009273/2576178‑2576471

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