Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,020,520 G→T 18.9% D212Y (GAT→TAT)  fre → NAD(P)H‑flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,020,5200GT18.9% 834.0 / 130.6 311D212Y (GAT→TAT) freNAD(P)H‑flavin reductase
Reads supporting (aligned to +/- strand):  ref base G (131/121);  new base T (32/27);  total (163/148)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 318 total aligned reads displayed.

TGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTAT  >  CP009273/4020371‑4020669
                                                                                                                                                     |                                                                                                                                                     
tGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCcgcg                                                                                                                                                       >  1:4616229/1‑150 (MQ=255)
   aaGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATc                                                                                                                                                    >  1:555250/1‑150 (MQ=255)
   aaGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATc                                                                                                                                                    >  1:852360/1‑150 (MQ=255)
         tgCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGtttt                                                                                                                                              <  2:4088837/150‑1 (MQ=255)
         tgCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGtttt                                                                                                                                              <  1:4940223/150‑1 (MQ=255)
           ccGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGc                                                                                                                                            <  2:2873068/150‑1 (MQ=255)
            tggtggtTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCa                                                                                                                                           >  1:4026321/2‑150 (MQ=255)
              gtggtTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATCTCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGt                                                                                                                                         >  2:4816429/1‑150 (MQ=255)
              gtggtTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTAATGCAGGATCACGGTACGCTGGCAGAGCATGATCTCTATATTGCCGGACGTTTTGAGATTGCGAAAATTGCCCTCGATCAGATTTGCAGt                                                                                                                                         >  1:1749455/1‑150 (MQ=255)
               tggtTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTg                                                                                                                                        <  1:3633413/150‑1 (MQ=255)
               tggtTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTg                                                                                                                                        <  2:3353195/150‑1 (MQ=255)
                  ttGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGc                                                                                                                                     <  2:2341054/150‑1 (MQ=255)
                   tGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCg                                                                                                                                    <  2:1069244/150‑1 (MQ=255)
                   tGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCg                                                                                                                                    >  2:3568700/1‑150 (MQ=255)
                   tGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCg                                                                                                                                    >  1:3801680/1‑150 (MQ=255)
                      acaacCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTaa                                                                                                                                 >  1:4174494/1‑150 (MQ=255)
                        aacCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATg                                                                                                                               <  2:3755670/150‑1 (MQ=255)
                          ccGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATgcg                                                                                                                             <  1:77892/150‑1 (MQ=255)
                           cGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATgcgc                                                                                                                            <  2:163048/150‑1 (MQ=255)
                            ggAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTGTTGCAGTGAGCGTAATGcgcg                                                                                                                           >  2:4326159/1‑150 (MQ=255)
                             gAAGCGGGCTGGCGTGGGCGTACTGGCACCTTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATAATGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTTCAGTGAGCGTAATGCGCgg                                                                                                                          >  1:4290253/1‑150 (MQ=255)
                             gAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCgg                                                                                                                          <  2:3193860/150‑1 (MQ=255)
                                  gggCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGa                                                                                                                     <  1:165535/150‑1 (MQ=255)
                                   ggCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCATAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCTCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGAt                                                                                                                    <  2:348696/150‑1 (MQ=255)
                                    gCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATc                                                                                                                   <  2:873083/150‑1 (MQ=255)
                                      tGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGGAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGc                                                                                                                 >  1:2932749/1‑150 (MQ=255)
                                      tGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGc                                                                                                                 >  2:2733451/1‑150 (MQ=255)
                                         cGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGACGATCGCCTg                                                                                                              >  2:3162240/1‑150 (MQ=255)
                                            gggCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGttt                                                                                                           >  2:2070993/1‑150 (MQ=255)
                                            gggCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGttt                                                                                                           >  2:2433804/1‑150 (MQ=255)
                                            gggCGGACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGttt                                                                                                           <  2:3376367/150‑1 (MQ=255)
                                                gTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCg                                                                                                       >  2:4656513/1‑150 (MQ=255)
                                                   cTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATg                                                                                                    <  1:4548903/150‑1 (MQ=255)
                                                    tGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGc                                                                                                   <  1:2238665/150‑1 (MQ=255)
                                                    tGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGACGAGCGCCTGTTGGGCGATGc                                                                                                   >  2:5101155/1‑150 (MQ=255)
                                                     ggCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCg                                                                                                  >  1:2683393/1‑150 (MQ=255)
                                                      gCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGt                                                                                                 <  2:2549505/150‑1 (MQ=255)
                                                        aCCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGttt                                                                                               <  2:549660/150‑1 (MQ=255)
                                                        aCCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGttt                                                                                               <  1:158687/150‑1 (MQ=255)
                                                           gtgtTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCa                                                                                            <  1:3240200/150‑1 (MQ=255)
                                                              ttAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCAttt                                                                                         >  1:2789751/1‑150 (MQ=255)
                                                              ttAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGACCGCGATCTTTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGTCGATGAGTTTGCAttt                                                                                         >  1:3921609/1‑150 (MQ=255)
                                                                aaCGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTAt                                                                                       >  2:1926707/1‑150 (MQ=255)
                                                                aaCGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTAt                                                                                       >  1:3511236/1‑150 (MQ=255)
                                                                 aCGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATc                                                                                      >  1:1022626/1‑150 (MQ=255)
                                                                 aCGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATc                                                                                      <  1:1977149/150‑1 (MQ=255)
                                                                      ggTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTgaga                                                                                 <  2:2999488/150‑1 (MQ=255)
                                                                      ggTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTgaga                                                                                 <  2:1289898/150‑1 (MQ=255)
                                                                      ggTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTgaga                                                                                 >  2:1571888/1‑150 (MQ=255)
                                                                      ggTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTgaga                                                                                 <  2:3841786/150‑1 (MQ=255)
                                                                       gTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGat                                                                                >  1:2075608/1‑150 (MQ=255)
                                                                         aTTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGatat                                                                              >  1:479886/1‑150 (MQ=255)
                                                                          ttGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGAtata                                                                             >  1:942692/1‑150 (MQ=255)
                                                                           tGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATaa                                                                            >  1:2122882/1‑150 (MQ=255)
                                                                           tGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATaa                                                                            >  2:4250516/1‑150 (MQ=255)
                                                                            gCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATaaa                                                                           >  2:449134/1‑150 (MQ=255)
                                                                              aGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATaaaaa                                                                         >  2:1752161/1‑150 (MQ=255)
                                                                              aGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATaaaaa                                                                         >  1:115929/1‑150 (MQ=255)
                                                                                   cACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAccc                                                                    <  1:1927276/150‑1 (MQ=255)
                                                                                    aCGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCg                                                                   >  2:2136509/1‑150 (MQ=255)
                                                                                     cGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGAGATAAAAAAACCCGc                                                                  >  2:4395657/1‑150 (MQ=255)
                                                                                        tACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGcccc                                                               >  2:2796367/1‑150 (MQ=255)
                                                                                          cGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTg                                                             <  1:2293005/150‑1 (MQ=255)
                                                                                                    gCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGga                                                   <  1:3564132/150‑1 (MQ=255)
                                                                                                     cATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGgaa                                                  >  2:1973252/1‑150 (MQ=255)
                                                                                                      aTGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGgaag                                                 <  2:4374437/150‑1 (MQ=255)
                                                                                                      aTGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGgaag                                                 <  2:2642587/150‑1 (MQ=255)
                                                                                                      aTGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGgaag                                                 <  1:2755712/150‑1 (MQ=255)
                                                                                                         atatCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAAc                                              <  1:787664/150‑1 (MQ=255)
                                                                                                           atCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACgg                                            <  2:3129368/150‑1 (MQ=255)
                                                                                                            tCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGc                                           <  1:3442413/150‑1 (MQ=255)
                                                                                                               atatTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAAc                                        >  1:1803016/1‑150 (MQ=255)
                                                                                                                  ttGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTaa                                     >  2:222848/1‑150 (MQ=255)
                                                                                                                   tGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTaaa                                    >  1:4241682/1‑150 (MQ=255)
                                                                                                                   tGCCGGACGATTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTaaa                                    <  1:1461143/150‑1 (MQ=255)
                                                                                                                    gCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAAc                                   >  2:4504599/1‑150 (MQ=255)
                                                                                                                     ccGGACGTTTTGAGATGGCGAAAATAGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACt                                  >  1:4859194/1‑150 (MQ=255)
                                                                                                                      cGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTg                                 >  2:4747453/1‑150 (MQ=255)
                                                                                                                      cGGAAGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACtt                                 <  2:619405/150‑2 (MQ=255)
                                                                                                                           gTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCTTAAAGCGAGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCACTGACAGGCGTGAATAACGGCAACTAAACTGTTAtt                            >  2:2485219/1‑150 (MQ=255)
                                                                                                                           gTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTAtt                            <  1:2601677/150‑1 (MQ=255)
                                                                                                                            ttttGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTc                           <  1:3619097/150‑1 (MQ=255)
                                                                                                                             tttGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCa                          >  2:345646/1‑150 (MQ=255)
                                                                                                                              ttGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAg                         >  1:710934/1‑150 (MQ=255)
                                                                                                                              ttGAGATGGCGAAAATTGCCCGCTATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAg                         >  2:2960254/1‑150 (MQ=255)
                                                                                                                                gagaTGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTg                       <  1:3368249/150‑1 (MQ=255)
                                                                                                                                 agaTGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTgg                      >  1:3065035/1‑150 (MQ=255)
                                                                                                                                    tGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCAt                   <  2:4323318/150‑1 (MQ=255)
                                                                                                                                    tGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCAt                   >  2:397116/1‑150 (MQ=255)
                                                                                                                                    tGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCAt                   >  2:4875772/1‑150 (MQ=255)
                                                                                                                                      gCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCAttt                 <  1:1729711/150‑1 (MQ=255)
                                                                                                                                             ttGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTa          >  1:2504286/1‑150 (MQ=255)
                                                                                                                                             ttGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTa          <  1:137417/150‑1 (MQ=255)
                                                                                                                                               gCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATg        <  1:1312150/150‑1 (MQ=255)
                                                                                                                                                cccGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGa       >  2:2111296/1‑150 (MQ=255)
                                                                                                                                                cccGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGa       >  2:1524393/1‑150 (MQ=255)
                                                                                                                                                 ccGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGAc      <  1:372137/150‑1 (MQ=255)
                                                                                                                                                  cgcgATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACg     >  1:1072922/1‑150 (MQ=255)
                                                                                                                                                  cgcgATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACg     >  2:4993307/1‑150 (MQ=255)
                                                                                                                                                     gATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTAt  >  1:3784616/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCGAAAATTGCCCGCGATCTGTTTTGCAGTGAGCGTAATGCGCGGGAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTAT  >  CP009273/4020371‑4020669

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: