Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 592,124 A→T 6.6% I357F (ATT→TTT)  cusC → copper/silver efflux system, outer membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273592,1240AT 6.6% 1326.1 / 62.3 362I357F (ATT→TTT) cusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base A (178/160);  new base T (10/14);  total (188/174)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 367 total aligned reads displayed.

TATAAGCCTGACCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGC  >  CP009273/591976‑592271
                                                                                                                                                    |                                                                                                                                                   
 ataAGCCTGACCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAAt                                                                                                                                                    >  1:698353/1‑149 (MQ=255)
  taAGCCTGACCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAAtt                                                                                                                                                   <  1:3049089/149‑1 (MQ=255)
    aGCCTGACCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCg                                                                                                                                                 >  2:1481995/1‑149 (MQ=255)
       cTGACCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCca                                                                                                                                              <  2:302413/149‑1 (MQ=255)
         gACCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCcagc                                                                                                                                            <  2:3330716/149‑1 (MQ=255)
          aCCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCcagca                                                                                                                                           >  2:4471110/1‑149 (MQ=255)
           ccAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATCTCGCCGAAATTCGCcagcag                                                                                                                                          >  1:4373977/1‑149 (MQ=255)
           ccAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCcagcag                                                                                                                                          <  2:145261/149‑1 (MQ=255)
            cAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCcagcagc                                                                                                                                         >  2:1675008/1‑149 (MQ=255)
             aGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCcagcagca                                                                                                                                        >  1:1572626/1‑149 (MQ=255)
              gCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCcagcagcag                                                                                                                                       <  1:1340435/149‑1 (MQ=255)
                  aaTATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCgg                                                                                                                                   <  2:3057765/149‑1 (MQ=255)
                   atatCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCggt                                                                                                                                  >  1:1913860/1‑149 (MQ=255)
                    tatCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCggtg                                                                                                                                 <  2:3987171/149‑1 (MQ=255)
                     atCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCggtgg                                                                                                                                <  2:2961679/149‑1 (MQ=255)
                       cGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGgtggtg                                                                                                                              <  2:1005908/149‑1 (MQ=255)
                        gACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGa                                                                                                                             >  2:2349526/1‑149 (MQ=255)
                         aCCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGaa                                                                                                                            <  1:2408921/149‑1 (MQ=255)
                          ccGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAAt                                                                                                                           >  1:4106562/1‑149 (MQ=255)
                          ccGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCAGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCTCCAGCAGCAGTCGGTGGTGAAt                                                                                                                           >  1:3113162/1‑149 (MQ=255)
                              cAGTAGCGATCTCTCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCGGGATATCGCCGAATTTCGCCAGCAGCAGTCGGTGGTGAATTATg                                                                                                                       <  2:2054970/149‑1 (MQ=255)
                              cAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATg                                                                                                                       <  1:2624064/149‑1 (MQ=255)
                                  aGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACa                                                                                                                   <  2:2519072/149‑1 (MQ=255)
                                   gCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATTCAGGACGAAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAg                                                                                                                  >  1:3687473/1‑149 (MQ=255)
                                   gCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAg                                                                                                                  >  1:4630815/1‑149 (MQ=255)
                                         tATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATc                                                                                                            >  1:3271218/1‑149 (MQ=255)
                                          aTCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATcc                                                                                                           <  1:2019963/149‑1 (MQ=255)
                                           tCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCa                                                                                                          <  2:1089995/149‑1 (MQ=255)
                                           tCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCa                                                                                                          >  2:4610344/1‑149 (MQ=255)
                                           tCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCa                                                                                                          >  1:243525/1‑149 (MQ=255)
                                           tCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCa                                                                                                          >  1:265931/1‑149 (MQ=255)
                                             gTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGa                                                                                                        >  2:2904383/1‑149 (MQ=255)
                                             gTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGa                                                                                                        <  1:984534/149‑1 (MQ=255)
                                              tCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGaa                                                                                                       <  1:4500429/149‑1 (MQ=255)
                                              tCATTATATAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGaa                                                                                                       <  1:338150/149‑1 (MQ=255)
                                               cATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAAc                                                                                                      >  2:5014653/1‑149 (MQ=255)
                                               cATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAAc                                                                                                      <  1:4901685/149‑1 (MQ=255)
                                                 ttattTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGc                                                                                                    >  1:1722091/1‑149 (MQ=255)
                                                  tattTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGcc                                                                                                   >  2:3709969/1‑149 (MQ=255)
                                                   attTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCt                                                                                                  <  2:2758973/149‑1 (MQ=255)
                                                    tttAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCtt                                                                                                 <  1:4031853/149‑1 (MQ=255)
                                                           ccAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAATTTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAaagaa                                                                                          >  2:1618127/1‑149 (MQ=255)
                                                           ccAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAaagaa                                                                                          >  2:1515166/1‑149 (MQ=255)
                                                           ccAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAaagaa                                                                                          >  2:1779772/1‑149 (MQ=255)
                                                            cagcagCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAaagaag                                                                                         >  2:313263/1‑149 (MQ=255)
                                                             agcagcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGt                                                                                        <  2:4824997/149‑1 (MQ=255)
                                                              gcagcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTg                                                                                       >  1:1249743/1‑149 (MQ=255)
                                                              gcagcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTg                                                                                       >  1:5088029/1‑149 (MQ=255)
                                                               cagcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTgg                                                                                      >  1:2848042/1‑149 (MQ=255)
                                                                 gcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCa                                                                                    >  2:3865933/1‑149 (MQ=255)
                                                                 gcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCa                                                                                    >  1:4234617/1‑149 (MQ=255)
                                                                   gggATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGa                                                                                  <  2:4225588/149‑1 (MQ=255)
                                                                     gATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATg                                                                                <  2:2632270/149‑1 (MQ=255)
                                                                      aTGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATgc                                                                               >  2:796804/1‑149 (MQ=255)
                                                                           gAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTgg                                                                          <  1:1078024/149‑1 (MQ=255)
                                                                             aTTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCa                                                                        >  2:4973459/1‑149 (MQ=255)
                                                                              ttttATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCAt                                                                       >  2:3953681/1‑149 (MQ=255)
                                                                              ttttATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCAt                                                                       <  1:2260372/149‑1 (MQ=255)
                                                                                ttATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTa                                                                     <  1:1585614/149‑1 (MQ=255)
                                                                                  ataCCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACg                                                                   >  2:138990/1‑149 (MQ=255)
                                                                                  ataCCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACg                                                                   >  1:3943033/1‑149 (MQ=255)
                                                                                   taCCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGt                                                                  <  1:2154332/149‑1 (MQ=255)
                                                                                   taCCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGt                                                                  >  1:4721345/1‑149 (MQ=255)
                                                                                    aCCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTc                                                                 <  2:353638/149‑1 (MQ=255)
                                                                                          aaTTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGcc                                                           <  2:319370/149‑1 (MQ=255)
                                                                                          aaTTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGcc                                                           >  1:595871/1‑149 (MQ=255)
                                                                                               agaTCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAAc                                                      >  2:5184806/1‑149 (MQ=255)
                                                                                                gaTCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACg                                                     >  1:2719417/1‑149 (MQ=255)
                                                                                                  tCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGAt                                                   >  2:1566318/1‑149 (MQ=255)
                                                                                                   ccccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATc                                                  >  2:1165518/1‑149 (MQ=255)
                                                                                                     ccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCaa                                                >  1:2497040/1‑149 (MQ=255)
                                                                                                      cATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCaaa                                               <  2:2426139/149‑1 (MQ=255)
                                                                                                       aTTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAATTTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAAt                                              <  2:2089803/149‑1 (MQ=255)
                                                                                                       aTTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAAt                                              <  1:1146763/149‑1 (MQ=255)
                                                                                                       aTTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAAt                                              <  2:441105/149‑1 (MQ=255)
                                                                                                        tttttAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAAtt                                             >  2:520035/1‑149 (MQ=255)
                                                                                                               tGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGccc                                      >  1:1529369/1‑149 (MQ=255)
                                                                                                                 ccGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCcag                                    <  1:4813043/149‑1 (MQ=255)
                                                                                                                  cGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCcagc                                   >  2:4933849/1‑149 (MQ=255)
                                                                                                                     aCGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCagcagc                                <  2:3616656/149‑1 (MQ=255)
                                                                                                                     aCGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCagcagc                                >  1:2082145/1‑149 (MQ=255)
                                                                                                                     aCGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCagcagc                                >  1:974285/1‑149 (MQ=255)
                                                                                                                      cGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCg                               >  1:4981128/1‑149 (MQ=255)
                                                                                                                          aCCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTAt                           <  1:5161442/149‑1 (MQ=255)
                                                                                                                           ccAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGTGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATc                          <  1:3745131/149‑1 (MQ=255)
                                                                                                                             aGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTg                        <  2:282409/149‑1 (MQ=255)
                                                                                                                              ggCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTgg                       >  1:3848896/1‑149 (MQ=255)
                                                                                                                                 cAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCAt                    >  1:1612333/1‑149 (MQ=255)
                                                                                                                                  aaTCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATc                   <  2:668888/149‑1 (MQ=255)
                                                                                                                                  aaTCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATc                   <  1:288670/149‑1 (MQ=255)
                                                                                                                                  aaTCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATc                   <  2:3934738/149‑1 (MQ=255)
                                                                                                                                  aaTCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATc                   <  1:3600148/149‑1 (MQ=255)
                                                                                                                                   aTCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCg                  >  1:526547/1‑149 (MQ=255)
                                                                                                                                     cTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCt                <  1:78201/149‑1 (MQ=255)
                                                                                                                                     cTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCt                <  1:1031315/149‑1 (MQ=255)
                                                                                                                                        gATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCa             >  2:4312467/1‑149 (MQ=255)
                                                                                                                                        gATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCa             >  1:276964/1‑149 (MQ=255)
                                                                                                                                         atatCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCaa            >  1:3328950/1‑149 (MQ=255)
                                                                                                                                          tatCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCaaa           <  1:910223/149‑1 (MQ=255)
                                                                                                                                             cGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTa        <  2:4120940/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATAAGCCTGACCAGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGC  >  CP009273/591976‑592271

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: