Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,458,948 G→T 100% D503Y (GAT→TAT)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,9480GT98.4% 606.2 / 0.8 189D503Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (109/77);  total (110/79)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 189 total aligned reads displayed.

TGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAA  >  CP009273/4458802‑4459092
                                                                                                                                                  |                                                                                                                                                
 ggtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTAt                                                                                                                                                <  1:26422/148‑1 (MQ=255)
      tGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCt                                                                                                                                           >  2:3227139/1‑148 (MQ=255)
        gCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGa                                                                                                                                         >  1:1068743/1‑148 (MQ=255)
         cAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGaa                                                                                                                                        >  2:1354329/1‑148 (MQ=255)
          aGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAAt                                                                                                                                       <  1:2917468/148‑1 (MQ=255)
           gAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATc                                                                                                                                      >  1:2067106/1‑148 (MQ=255)
            aaaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTCTACGCTGAATCg                                                                                                                                     >  2:958378/1‑148 (MQ=255)
            aaaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCg                                                                                                                                     >  1:886834/1‑148 (MQ=255)
            aaaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCg                                                                                                                                     <  2:3491593/148‑1 (MQ=255)
            aaaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCg                                                                                                                                     >  2:307304/1‑148 (MQ=255)
               aTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCa                                                                                                                                  <  2:1614350/148‑1 (MQ=255)
                  ccGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAggg                                                                                                                               >  1:3255891/1‑148 (MQ=255)
                      gCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTg                                                                                                                           <  1:685009/148‑1 (MQ=255)
                      gCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTg                                                                                                                           >  2:3429062/1‑148 (MQ=255)
                          cATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGcgcg                                                                                                                       <  1:1537716/148‑1 (MQ=255)
                           aTCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGt                                                                                                                      <  1:2809861/148‑1 (MQ=255)
                            tCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTg                                                                                                                     >  1:1914521/1‑148 (MQ=255)
                              aGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGt                                                                                                                   <  1:3073937/148‑1 (MQ=255)
                              aGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGt                                                                                                                   <  2:1944546/148‑1 (MQ=255)
                              aGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGt                                                                                                                   <  2:3113082/148‑1 (MQ=255)
                               gCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGtt                                                                                                                  >  1:969542/1‑148 (MQ=255)
                                  ggTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCg                                                                                                               <  2:169110/148‑1 (MQ=255)
                                   cattGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATAGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTCGTTGCgg                                                                                                              >  2:2003670/3‑148 (MQ=255)
                                    tttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGt                                                                                                             >  2:1323604/1‑148 (MQ=255)
                                     ttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                            >  2:62895/1‑148 (MQ=255)
                                      tGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTgg                                                                                                           <  1:2418690/148‑1 (MQ=255)
                                        cAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcg                                                                                                         >  1:3263745/1‑148 (MQ=255)
                                          aaGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcgac                                                                                                       <  1:3289203/148‑1 (MQ=255)
                                          aaGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcgac                                                                                                       >  1:3487223/1‑148 (MQ=255)
                                             gTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaa                                                                                                    >  2:948071/1‑148 (MQ=255)
                                             gTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaa                                                                                                    >  2:947909/1‑148 (MQ=255)
                                               gCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAAt                                                                                                  <  2:2901859/148‑1 (MQ=255)
                                                    gCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTg                                                                                             >  2:3321586/1‑148 (MQ=255)
                                                        gAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAg                                                                                         >  1:179176/1‑148 (MQ=255)
                                                        gAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAg                                                                                         >  1:132715/1‑148 (MQ=255)
                                                        gAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAg                                                                                         >  1:1946334/1‑148 (MQ=255)
                                                          aaTCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAgcg                                                                                       >  1:2797013/1‑148 (MQ=255)
                                                             ccTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                    >  1:1665387/1‑148 (MQ=255)
                                                               tCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGa                                                                                  >  2:3142441/1‑148 (MQ=255)
                                                               tCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGa                                                                                  >  1:1023449/1‑148 (MQ=255)
                                                                  aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAgg                                                                               >  2:443141/1‑148 (MQ=255)
                                                                    gtgtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGgag                                                                             >  2:1950367/1‑148 (MQ=255)
                                                                    gtgtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGgag                                                                             <  1:947909/148‑1 (MQ=255)
                                                                    gtgtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGgag                                                                             <  1:948071/148‑1 (MQ=255)
                                                                      gtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAt                                                                           >  1:1028353/1‑148 (MQ=255)
                                                                         ttCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTAc                                                                        <  2:1276113/148‑1 (MQ=255)
                                                                          tCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTAcc                                                                       <  1:3121190/148‑1 (MQ=255)
                                                                          tCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTAcc                                                                       >  2:1564758/1‑148 (MQ=255)
                                                                            gCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAg                                                                     >  1:1210414/1‑148 (MQ=255)
                                                                             cAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgc                                                                    >  1:775124/1‑148 (MQ=255)
                                                                             cAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgc                                                                    >  1:2963608/1‑148 (MQ=255)
                                                                                          aTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATc                                                       >  1:2279873/1‑148 (MQ=255)
                                                                                           tGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATcc                                                      >  1:1436768/1‑148 (MQ=255)
                                                                                             gCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGa                                                    <  1:451465/148‑1 (MQ=255)
                                                                                              cGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCTGCGGGCGGATGAATCCGAc                                                   <  2:1898806/148‑1 (MQ=255)
                                                                                              cGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAc                                                   <  2:2279770/148‑1 (MQ=255)
                                                                                               gagaTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAcc                                                  <  2:3255891/148‑1 (MQ=255)
                                                                                                  aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGa                                               <  1:2403204/148‑1 (MQ=255)
                                                                                                  aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGa                                               <  1:3333665/148‑1 (MQ=255)
                                                                                                    tgtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATc                                             >  2:1992450/1‑148 (MQ=255)
                                                                                                         cGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGa                                        <  1:1819203/148‑1 (MQ=255)
                                                                                                            tCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAgg                                     <  2:179176/148‑1 (MQ=255)
                                                                                                              gATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGat                                   <  1:424014/148‑1 (MQ=255)
                                                                                                               aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  <  1:3450792/148‑1 (MQ=255)
                                                                                                               aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  <  1:2491114/148‑1 (MQ=255)
                                                                                                                tGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatat                                 >  2:1230684/1‑148 (MQ=255)
                                                                                                                 gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatata                                >  2:2704673/1‑148 (MQ=255)
                                                                                                                  aCATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatatat                               >  1:737615/1‑148 (MQ=255)
                                                                                                                   catcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATAtt                              >  1:1479551/1‑148 (MQ=255)
                                                                                                                    atcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTg                             <  2:2755208/148‑1 (MQ=255)
                                                                                                                    atcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTg                             >  2:2028172/1‑148 (MQ=255)
                                                                                                                    atcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTg                             >  2:3368292/1‑148 (MQ=255)
                                                                                                                      catGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCt                           >  1:1376049/1‑148 (MQ=255)
                                                                                                                        tgctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCttt                         >  2:2493193/1‑148 (MQ=255)
                                                                                                                        tgctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCttt                         >  2:2867406/1‑148 (MQ=255)
                                                                                                                        tgctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCttt                         >  1:2875104/1‑148 (MQ=255)
                                                                                                                         gctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCtttt                        >  2:871386/1‑148 (MQ=255)
                                                                                                                         gctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCtttt                        <  2:2643176/148‑1 (MQ=255)
                                                                                                                          ctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTc                       <  1:1498124/148‑1 (MQ=255)
                                                                                                                           tgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCt                      <  1:3321586/148‑1 (MQ=255)
                                                                                                                            gcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCtt                     >  1:1624234/1‑148 (MQ=255)
                                                                                                                            gcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACATGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGATTTTCtt                     >  1:339243/1‑148 (MQ=255)
                                                                                                                               tAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                  >  1:1511023/1‑148 (MQ=255)
                                                                                                                               tAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                  >  2:1631123/1‑148 (MQ=255)
                                                                                                                                aaGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATc                 >  1:1703723/1‑148 (MQ=255)
                                                                                                                                aaGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATc                 <  2:2113408/148‑1 (MQ=255)
                                                                                                                                 aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                >  1:1242957/1‑148 (MQ=255)
                                                                                                                                 aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                >  2:2627146/1‑148 (MQ=255)
                                                                                                                                 aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                >  1:3203622/1‑148 (MQ=255)
                                                                                                                                  gATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCg               <  2:1068743/148‑1 (MQ=255)
                                                                                                                                     tAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCa            >  1:3416939/1‑148 (MQ=255)
                                                                                                                                      aaGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCaa           <  2:1771057/148‑1 (MQ=255)
                                                                                                                                         cGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAag        >  2:1012834/1‑148 (MQ=255)
                                                                                                                                         cGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAag        >  2:1013649/1‑148 (MQ=255)
                                                                                                                                           ggTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAagag      >  2:2970977/1‑148 (MQ=255)
                                                                                                                                            gTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAgaga     >  1:435996/1‑148 (MQ=255)
                                                                                                                                             ttACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGAc    >  2:2909815/1‑148 (MQ=255)
                                                                                                                                             ttACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGAc    >  2:1160292/1‑148 (MQ=255)
                                                                                                                                               aCTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACaa  <  1:1872031/148‑1 (MQ=255)
                                                                                                                                               aCTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACaa  <  1:3227137/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAA  >  CP009273/4458802‑4459092

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