Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 591,864 C→A 90.1% A270E (GCG→GAG)  cusC → copper/silver efflux system, outer membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273591,8640CA90.1% 426.4 / 43.9 172A270E (GCG→GAG) cusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base C (11/6);  new base A (75/80);  total (86/86)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 174 total aligned reads displayed.

GAAAGTACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGAATATCGACCGC  >  CP009273/591719‑592005
                                                                                                                                                 |                                                                                                                                             
gAAAGTACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACcggagg                                                                                                                                             <  1:3240987/148‑1 (MQ=255)
  aaGTACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGa                                                                                                                                           <  1:1032085/148‑1 (MQ=255)
  aaGTACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGa                                                                                                                                           <  1:2315588/148‑1 (MQ=255)
    gTACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGAtt                                                                                                                                         >  1:166310/1‑148 (MQ=255)
     tACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTg                                                                                                                                        >  2:2937607/1‑148 (MQ=255)
     tACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTg                                                                                                                                        >  2:669096/1‑148 (MQ=255)
     tACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTg                                                                                                                                        >  1:3303306/1‑148 (MQ=255)
      aCCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTgt                                                                                                                                       >  2:190907/1‑148 (MQ=255)
       cccGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTgtc                                                                                                                                      >  2:1788333/1‑148 (MQ=255)
        ccGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTgtcg                                                                                                                                     <  2:214362/148‑1 (MQ=255)
        ccGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTgtcg                                                                                                                                     <  2:214392/148‑1 (MQ=255)
           cAGCGACATCGCTACAAGTCAGGGGGAACTGGCGCAGGAGAATTATGCTTTGCAACTGTTGCTGTGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAATCAGCGACAGCCTGCAAAGCGTTTAATTACCGGAGGGATTGtcgtcg                                                                                                                                  <  2:1572500/148‑1 (MQ=255)
           cAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGtcgtcg                                                                                                                                  <  2:1572125/148‑1 (MQ=255)
              cGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCaa                                                                                                                               >  1:2832602/1‑148 (MQ=255)
               gACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCaaa                                                                                                                              <  1:2471186/148‑1 (MQ=255)
               gACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCaaa                                                                                                                              <  1:2471151/148‑1 (MQ=255)
                aCATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAAt                                                                                                                             <  2:1393819/148‑1 (MQ=255)
                 cATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATc                                                                                                                            >  2:518438/1‑148 (MQ=255)
                   tCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCtt                                                                                                                          >  1:2818669/1‑148 (MQ=255)
                   tCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCtt                                                                                                                          >  1:2818624/1‑148 (MQ=255)
                    cGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTa                                                                                                                         >  2:1000589/1‑148 (MQ=255)
                     gCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTAt                                                                                                                        >  1:265149/1‑148 (MQ=255)
                      cTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTAtt                                                                                                                       >  2:1527072/1‑148 (MQ=255)
                       tAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTg                                                                                                                      <  1:435597/148‑1 (MQ=255)
                         aaCGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCa                                                                                                                    >  1:3149003/1‑148 (MQ=255)
                             tCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAgcgc                                                                                                                >  1:646843/1‑148 (MQ=255)
                                    gAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGata                                                                                                         >  1:1837677/1‑148 (MQ=255)
                                      aCTGGCGCAGGCTAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACATTAAACAGCGACAGCCTGCAAAGCGTTAAATTCCCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGAtatt                                                                                                       <  1:3072785/148‑1 (MQ=255)
                                      aCTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGAtatt                                                                                                       <  2:1199476/148‑1 (MQ=255)
                                       cTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGAtatta                                                                                                      <  2:2387169/148‑1 (MQ=255)
                                       cTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGAtatta                                                                                                      <  2:1303366/148‑1 (MQ=255)
                                          gcgcCGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGAAAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATgg                                                                                                   <  1:1172101/148‑1 (MQ=255)
                                          gcgcAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATgg                                                                                                   >  2:600/1‑148 (MQ=255)
                                           cgcAGGCTAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGa                                                                                                  <  1:997484/148‑1 (MQ=255)
                                           cgcAGGCTAATAATGCATTGCAACTGTTGCTGGGAAGCTACGCCAAGCTGCCGCAAGCGCATACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGa                                                                                                  <  1:2518374/148‑1 (MQ=255)
                                             cAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCACAGCGTTAAATTACCGCCGGGATTGTCTTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAg                                                                                                >  1:578786/1‑148 (MQ=255)
                                             cAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAg                                                                                                <  1:190907/148‑1 (MQ=255)
                                               ggCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCt                                                                                              >  1:1282662/1‑148 (MQ=255)
                                               ggCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCt                                                                                              >  1:1282648/1‑148 (MQ=255)
                                                    ataatGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAaca                                                                                         >  2:3141443/1‑148 (MQ=255)
                                                    ataatGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAaca                                                                                         <  1:1317504/148‑1 (MQ=255)
                                                    ataatGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAaca                                                                                         >  2:436760/1‑148 (MQ=255)
                                                    ataatGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAaca                                                                                         >  2:1007616/1‑148 (MQ=255)
                                                      aatGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACAcg                                                                                       >  1:2725499/1‑148 (MQ=255)
                                                       atGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACAcgc                                                                                      >  2:2206311/1‑148 (MQ=255)
                                                       atGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACAcgc                                                                                      >  1:535251/1‑148 (MQ=255)
                                                       atGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACAcgc                                                                                      >  2:2943525/1‑148 (MQ=255)
                                                           aTTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTa                                                                                  <  2:900740/148‑1 (MQ=255)
                                                           aTTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTa                                                                                  <  2:1838370/148‑1 (MQ=255)
                                                           aTTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTa                                                                                  <  2:1408014/148‑1 (MQ=255)
                                                            ttGCAACTGTTGCTGGGAAGCTACGGCATGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTaa                                                                                 <  1:838648/148‑1 (MQ=255)
                                                            ttGCAACTGGTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTaa                                                                                 <  1:2885857/148‑1 (MQ=255)
                                                             tGCAAGTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAATTTACCGGAGGGTTTGTCGTCGAAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAAt                                                                                <  1:2799393/148‑1 (MQ=255)
                                                             tGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAAt                                                                                >  1:3156351/1‑148 (MQ=255)
                                                             tGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCAGGATATTATGGAAGCTGAACACGCGTTAAt                                                                                <  1:2026405/148‑1 (MQ=255)
                                                                   tGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATggcggc                                                                          <  2:109972/148‑1 (MQ=255)
                                                                    gTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCt                                                                         >  1:2231816/1‑148 (MQ=255)
                                                                     ttGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTa                                                                        >  1:833876/1‑148 (MQ=255)
                                                                       gCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAAt                                                                      <  2:824609/148‑1 (MQ=255)
                                                                          gggAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGcc                                                                   >  1:2976958/1‑148 (MQ=255)
                                                                            gAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCaa                                                                 <  2:2425841/148‑1 (MQ=255)
                                                                             aaGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAat                                                                <  1:600/148‑1 (MQ=255)
                                                                             aaGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAat                                                                >  1:3036492/1‑148 (MQ=255)
                                                                               gCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAatat                                                              <  1:3130000/148‑1 (MQ=255)
                                                                                cTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATAtt                                                             >  1:337993/1‑148 (MQ=255)
                                                                                         gCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCa                                                    >  1:514595/1‑148 (MQ=255)
                                                                                           tGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCAcg                                                  <  2:146577/148‑1 (MQ=255)
                                                                                                   gcgcAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCAt                                          <  1:1002034/148‑1 (MQ=255)
                                                                                                   gcgcAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCAt                                          <  1:1000589/148‑1 (MQ=255)
                                                                                                   gcgcAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCAt                                          <  1:3251908/148‑1 (MQ=255)
                                                                                                   gcgcAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCAt                                          <  1:2878919/148‑1 (MQ=255)
                                                                                                   gcgcAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCAt                                          <  1:2223236/148‑1 (MQ=255)
                                                                                                        gACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCAtttttt                                     >  2:384348/1‑148 (MQ=255)
                                                                                                         aCAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTc                                    >  2:2531859/1‑148 (MQ=255)
                                                                                                          cAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTcc                                   >  1:1242354/1‑148 (MQ=255)
                                                                                                            gTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGt                                 <  1:518438/148‑1 (MQ=255)
                                                                                                            gTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGt                                 <  2:1591135/148‑1 (MQ=255)
                                                                                                            gTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGt                                 <  1:29073/148‑1 (MQ=255)
                                                                                                             tAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTc                                <  1:2047559/148‑1 (MQ=255)
                                                                                                              aaaCAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCt                               <  2:514582/148‑1 (MQ=255)
                                                                                                               aaCAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCta                              <  1:746992/148‑1 (MQ=255)
                                                                                                               aaCAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCta                              <  1:2405783/148‑1 (MQ=255)
                                                                                                                aCAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCtat                             >  2:2503759/1‑148 (MQ=255)
                                                                                                                aCAGCGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCtat                             >  2:240595/1‑148 (MQ=255)
                                                                                                                    cGACAGCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGc                         <  2:422839/148‑1 (MQ=255)
                                                                                                                         gCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGAc                    >  2:2370478/1‑148 (MQ=255)
                                                                                                                         gCCTGCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGAc                    <  1:1154052/148‑1 (MQ=255)
                                                                                                                             gCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGc                <  1:2026118/148‑1 (MQ=255)
                                                                                                                             gCAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGc                >  1:2177586/1‑148 (MQ=255)
                                                                                                                              cAAAGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCg               <  1:2757864/148‑1 (MQ=255)
                                                                                                                                aaGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGa             <  2:2780928/148‑1 (MQ=255)
                                                                                                                                aaGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGa             >  1:1430407/1‑148 (MQ=255)
                                                                                                                                 aGCGTTAAATTACCGGAGGGTTTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGaa            <  1:489337/148‑1 (MQ=255)
                                                                                                                                 aGCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGaa            >  1:2892203/1‑148 (MQ=255)
                                                                                                                                  gCGTTAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGAat           >  2:2737019/1‑148 (MQ=255)
                                                                                                                                     ttAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGAATATc        >  1:1785573/1‑148 (MQ=255)
                                                                                                                                     ttAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGAATATc        >  1:113438/1‑148 (MQ=255)
                                                                                                                                      tAAATTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGAATATCg       <  1:2875355/148‑1 (MQ=255)
                                                                                                                                         aTTACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGACTATCGAcc    >  2:1820284/1‑148 (MQ=255)
                                                                                                                                           tACCGGAGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGAATATCGACCGc  >  2:2595532/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GAAAGTACCCGCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCGCAAATCTTATTGCAGCGCCCGGATATTATGGAAGCTGAACACGCGTTAATGGCGGCTAATGCCAATATTGGCGCTGCACGCGCGGCATTTTTTCCGTCTATAAGCCTGACCAGCGGAATATCGACCGC  >  CP009273/591719‑592005

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: