Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,834,958 G→T 43.3% Q48K (CAG→AAG)  nepI ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,834,9580GT43.3% 251.8 / 295.9 269Q48K (CAG→AAG) nepIputative transporter
Reads supporting (aligned to +/- strand):  ref base G (102/48);  new base T (72/45);  total (174/95)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

GAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTAC  >  CP009273/3834809‑3835107
                                                                                                                                                     |                                                                                                                                                     
gagCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTg                                                                                                                                                       <  1:18266/150‑1 (MQ=255)
    gTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGcc                                                                                                                                                   >  1:13585/1‑150 (MQ=255)
    gTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGcc                                                                                                                                                   >  1:2318853/1‑150 (MQ=255)
    gTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGcc                                                                                                                                                   >  1:3999292/1‑150 (MQ=255)
      cAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCAt                                                                                                                                                 <  2:3556713/150‑1 (MQ=255)
       aGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCAtt                                                                                                                                                <  1:60082/150‑1 (MQ=255)
         cAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTgg                                                                                                                                              <  1:100688/150‑1 (MQ=255)
         cAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTgg                                                                                                                                              >  1:976712/1‑150 (MQ=255)
          aaaaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGc                                                                                                                                             >  1:2803691/1‑150 (MQ=255)
           aaaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCg                                                                                                                                            >  2:2151252/1‑150 (MQ=255)
           aaaCGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCg                                                                                                                                            >  2:2151066/1‑150 (MQ=255)
              cGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTca                                                                                                                                         >  1:4247918/1‑150 (MQ=255)
               ggCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTcaa                                                                                                                                        >  1:1149093/1‑150 (MQ=255)
               ggCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTcaa                                                                                                                                        >  1:3513116/1‑150 (MQ=255)
                gCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTcaac                                                                                                                                       <  2:104800/150‑1 (MQ=255)
                gCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTcaac                                                                                                                                       <  2:3714379/150‑1 (MQ=255)
                    acaaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTTCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAAc                                                                                                                                   >  2:4316671/1‑150 (MQ=255)
                     caaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACt                                                                                                                                  >  2:3620056/1‑150 (MQ=255)
                     caaaATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACt                                                                                                                                  <  2:2241715/150‑1 (MQ=255)
                      aaaaTAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTg                                                                                                                                 >  2:818749/1‑150 (MQ=255)
                      aaaaTAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTg                                                                                                                                 >  2:2031381/1‑150 (MQ=255)
                         aTAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACg                                                                                                                              >  2:1891547/1‑150 (MQ=255)
                          tAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACgg                                                                                                                             <  1:3355756/150‑1 (MQ=255)
                             caacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCa                                                                                                                          >  2:3728194/1‑150 (MQ=255)
                              aacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCaa                                                                                                                         >  1:586950/1‑150 (MQ=255)
                              aacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCaa                                                                                                                         >  1:1282841/1‑150 (MQ=255)
                              aacGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCaa                                                                                                                         >  2:2935294/1‑150 (MQ=255)
                                  tAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAAc                                                                                                                     >  2:556592/1‑150 (MQ=255)
                                   aGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACt                                                                                                                    >  1:256392/1‑150 (MQ=255)
                                    gcggcgGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTc                                                                                                                   >  1:579374/1‑150 (MQ=255)
                                     cggcggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCa                                                                                                                  <  2:3296946/150‑1 (MQ=255)
                                      ggcggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCaa                                                                                                                 >  2:370745/1‑150 (MQ=255)
                                       gcggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAAc                                                                                                                <  1:4318210/150‑1 (MQ=255)
                                       gcggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAAc                                                                                                                <  2:2646348/150‑1 (MQ=255)
                                        cggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACg                                                                                                               <  1:598032/150‑1 (MQ=255)
                                         ggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGa                                                                                                              >  1:3411885/1‑150 (MQ=255)
                                         ggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGa                                                                                                              >  1:2254525/1‑150 (MQ=255)
                                         ggTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAACGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGa                                                                                                              >  1:4474083/1‑150 (MQ=255)
                                            cAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATaa                                                                                                           >  1:3323444/1‑150 (MQ=255)
                                               tAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCa                                                                                                        <  2:1549147/150‑1 (MQ=255)
                                               tAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCa                                                                                                        <  2:2492243/150‑1 (MQ=255)
                                                aGCCTGAATTGTCTGGGTAATCAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAg                                                                                                       <  1:3472977/150‑1 (MQ=255)
                                                 gCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGa                                                                                                      >  2:2008043/1‑150 (MQ=255)
                                                  ccTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGAc                                                                                                     <  2:3495316/150‑1 (MQ=255)
                                                    tGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAg                                                                                                   >  1:4127797/1‑150 (MQ=255)
                                                      aaTTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGc                                                                                                 >  1:946307/1‑150 (MQ=255)
                                                       aTTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCg                                                                                                >  2:4103285/1‑150 (MQ=255)
                                                          gTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGaca                                                                                             >  1:1743273/1‑150 (MQ=255)
                                                            cTGGGTAATAAACAAACTTGCAAACATAGCCACAAAGGCGGTAACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGAcaca                                                                                           >  1:242880/1‑150 (MQ=255)
                                                            cTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGAcaca                                                                                           >  1:4353473/1‑150 (MQ=255)
                                                             tGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACaa                                                                                          >  1:3830050/1‑150 (MQ=255)
                                                             tGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACaa                                                                                          >  1:1519212/1‑150 (MQ=255)
                                                                gTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAAc                                                                                       <  2:3317324/150‑1 (MQ=255)
                                                                 taataaACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACg                                                                                      >  2:3332992/1‑150 (MQ=255)
                                                                 taataaACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACg                                                                                      >  1:783313/1‑150 (MQ=255)
                                                                    taaACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCa                                                                                   <  2:2452606/150‑1 (MQ=255)
                                                                      aacaaacTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCAcc                                                                                 >  1:1912947/1‑150 (MQ=255)
                                                                           acTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGaa                                                                            <  2:3999292/150‑1 (MQ=255)
                                                                                 aaaCATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGc                                                                      >  2:415114/1‑150 (MQ=255)
                                                                                 aaaCATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGc                                                                      >  2:3727926/1‑150 (MQ=255)
                                                                                       tGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAATCGCCACCGAGAAAACGGATGACCa                                                                >  2:644242/1‑150 (MQ=255)
                                                                                        gCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACAGAGAAAACGGCTGACCAg                                                               >  2:3275925/1‑150 (MQ=255)
                                                                                             aaaGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcgg                                                          >  2:2374147/1‑150 (MQ=255)
                                                                                             aaaGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcgg                                                          >  1:2503424/1‑150 (MQ=255)
                                                                                              aaGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGAGGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcggt                                                         >  1:3132666/1‑150 (MQ=255)
                                                                                                ggCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTtcggtcg                                                       >  2:2832320/1‑150 (MQ=255)
                                                                                                   ggtcacggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGt                                                    >  1:1033809/1‑150 (MQ=255)
                                                                                                      cacggtcacCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGAt                                                 >  2:784917/1‑150 (MQ=255)
                                                                                                       acggccacCGATTGCCCGGCAACCCCTTCCGCAATCCCCACATCCTCGGCCATTGGCGTCAACAAACCGACGGGCAAAAACTCCACGAGAATCAGCCACGCGACACACAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATg                                                <  2:3104560/150‑1 (MQ=255)
                                                                                                             acCGATTGCCCGGCACCCCCTTCCGAACTGCCCCAATCCTCGGCCATTGGCGTCAACAACCTGACGGGCAAACCCTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATcc                                          <  2:2254525/150‑1 (MQ=255)
                                                                                                             acCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATcc                                          <  2:3020032/150‑1 (MQ=255)
                                                                                                             acCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATcc                                          >  1:717860/1‑150 (MQ=255)
                                                                                                              ccGATTGCCCGGCAACCCCTTCCGAAACGCCTCAATCCTTGGCCATCGGCGACAACCCACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccg                                         <  2:1867632/150‑1 (MQ=255)
                                                                                                               cGATTGCCCGGCAACCCCTTCCGACATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgc                                        <  1:3498031/150‑1 (MQ=255)
                                                                                                                gATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcg                                       >  2:1672510/1‑150 (MQ=255)
                                                                                                                gATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcg                                       >  2:1672971/1‑150 (MQ=255)
                                                                                                                  ttGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcgcc                                     >  2:1695630/1‑150 (MQ=255)
                                                                                                                  ttGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcgcc                                     >  2:2727297/1‑150 (MQ=255)
                                                                                                                  ttGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcgcc                                     >  2:237300/1‑150 (MQ=255)
                                                                                                                     cccGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATccgcgccgcg                                  >  2:3010568/1‑150 (MQ=255)
                                                                                                                       cGGCAACCCCTTCCGAACTGCCTAAATCCTTGGCCATTGGCGTCAACCAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGt                                <  2:3030231/150‑1 (MQ=255)
                                                                                                                        ggCAACCCCTTCCGAAATGCCTCAATCCTGGGCCATTGGCGTCAACACACTGACGGGCAAAAACTCAACGAGAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGtt                               <  2:226994/150‑1 (MQ=255)
                                                                                                                          cAACCCCTTCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAACCTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGtttt                             <  2:2318853/150‑1 (MQ=255)
                                                                                                                          cAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGtttt                             >  1:1402985/1‑150 (MQ=255)
                                                                                                                                cTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCaa                       >  1:2949303/1‑150 (MQ=255)
                                                                                                                                  tCCGAAATGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATa                     <  1:556592/150‑1 (MQ=255)
                                                                                                                                     gAAATGCCGAAATCCTGGGCCAGTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAAt                  <  2:2803691/150‑1 (MQ=255)
                                                                                                                                      aaaTGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAAtt                 <  1:3728194/150‑1 (MQ=255)
                                                                                                                                       aaTGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTc                >  2:3194326/1‑150 (MQ=255)
                                                                                                                                       aaTGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTc                >  2:3195868/1‑150 (MQ=255)
                                                                                                                                        aTGCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCa               <  1:4017077/150‑1 (MQ=255)
                                                                                                                                          gCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGATAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACt             >  1:3690795/1‑150 (MQ=255)
                                                                                                                                          gCCTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACt             <  2:3924186/150‑1 (MQ=255)
                                                                                                                                          gCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGTGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACt             >  1:3669984/1‑150 (MQ=255)
                                                                                                                                           ccTAAATCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTc            >  1:1626950/1‑150 (MQ=255)
                                                                                                                                             tAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCAt          <  1:3176496/150‑1 (MQ=255)
                                                                                                                                             tAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCAt          >  2:1186469/1‑150 (MQ=255)
                                                                                                                                                aTCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGgt       >  2:1107207/1‑150 (MQ=255)
                                                                                                                                                aTCCTTGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGgt       >  2:4439102/1‑150 (MQ=255)
                                                                                                                                                  ccTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGgtgt     >  2:3256017/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAATAAACAAACTGGCAAACATTGCCACAAAGGCGGTCACGGTCACCGATTGCCCGGCAACCCCTTCCGAAATGCCTAAATCCTGGGCCATTGGCGTCAACAAACTGACGGGCAAAAACTCAACGATAATCAGACAGGCGACACAAAACGCCACCGAGAAAACGGCTGACCAGTTCGGTCGGGTGATGGCATCCGCGCCGCGGTTTTCGGCAATAAATTCACTCATGGTGTTAC  >  CP009273/3834809‑3835107

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