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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,6090GT98.9% 607.7 / ‑1.9 188D390Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (114/72);  total (116/72)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-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.

GGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC  >  CP009273/4458464‑4458753
                                                                                                                                                 |                                                                                                                                                
    cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATAtct                                                                                                                                            <  2:2962020/148‑1 (MQ=255)
        gcgcAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCgg                                                                                                                                        <  2:267595/148‑1 (MQ=255)
           cAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaa                                                                                                                                     >  1:2476774/1‑148 (MQ=255)
           cAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaa                                                                                                                                     >  1:690111/1‑148 (MQ=255)
             gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaaaa                                                                                                                                   <  2:1094537/148‑1 (MQ=255)
             gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaaaa                                                                                                                                   >  2:2857216/1‑148 (MQ=255)
             gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaaaa                                                                                                                                   >  2:1718503/1‑148 (MQ=255)
               aaaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAcc                                                                                                                                 >  1:3290541/1‑148 (MQ=255)
                aaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCa                                                                                                                                >  1:2935701/1‑148 (MQ=255)
                 aGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAg                                                                                                                               <  2:2585071/148‑1 (MQ=255)
                  gCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGc                                                                                                                              <  1:2328154/148‑1 (MQ=255)
                    tGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGa                                                                                                                            <  1:3096952/148‑1 (MQ=255)
                         aaaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAAcgcg                                                                                                                       >  2:3134034/1‑148 (MQ=255)
                         aaaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAAcgcg                                                                                                                       >  1:994647/1‑148 (MQ=255)
                          aaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATCAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGt                                                                                                                      >  2:3111356/1‑148 (MQ=255)
                          aaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGt                                                                                                                      >  1:2601057/1‑148 (MQ=255)
                           aCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAACCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtg                                                                                                                     <  1:5425/148‑1 (MQ=255)
                            cAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgc                                                                                                                    >  1:3371870/1‑148 (MQ=255)
                            cAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgc                                                                                                                    >  2:2819898/1‑148 (MQ=255)
                              gAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgctg                                                                                                                  >  1:3228546/1‑148 (MQ=255)
                              gAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgctg                                                                                                                  >  2:1267578/1‑148 (MQ=255)
                                aaGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCa                                                                                                                >  1:2933740/1‑148 (MQ=255)
                                   tGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAg                                                                                                             <  1:2158391/148‑1 (MQ=255)
                                    gATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgc                                                                                                            >  2:1188289/1‑148 (MQ=255)
                                    gATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgc                                                                                                            >  2:2982399/1‑148 (MQ=255)
                                        gTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGt                                                                                                        >  2:2536382/1‑148 (MQ=255)
                                           aaaCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGt                                                                                                     <  1:990435/148‑1 (MQ=255)
                                              cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGa                                                                                                  >  2:1126623/1‑148 (MQ=255)
                                                tCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                <  2:1308006/148‑1 (MQ=255)
                                                tCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                >  1:1939430/1‑148 (MQ=255)
                                                 cTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACa                                                                                               >  1:2435848/1‑148 (MQ=255)
                                                  tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAg                                                                                              <  2:2251478/148‑1 (MQ=255)
                                                   ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                             <  2:1860850/148‑1 (MQ=255)
                                                   ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                             <  1:3520787/148‑1 (MQ=255)
                                                   ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGAGCGTGGTTGAACAGt                                                                                             >  2:1907276/1‑148 (MQ=255)
                                                        cTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTa                                                                                        <  1:1585550/148‑1 (MQ=255)
                                                            tCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAg                                                                                    <  1:1206422/148‑1 (MQ=255)
                                                            tCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAg                                                                                    >  1:819947/1‑148 (MQ=255)
                                                                  cTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGa                                                                              >  1:3405004/1‑148 (MQ=255)
                                                                   tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGAc                                                                             <  2:833749/148‑1 (MQ=255)
                                                                      ggCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTa                                                                          >  1:1153083/1‑148 (MQ=255)
                                                                      ggCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTa                                                                          >  1:2930031/1‑148 (MQ=255)
                                                                      ggCGCAATGGATATTCTGTGAACTGATAAAACCGGCACACTGACTAAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCAGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTa                                                                          >  1:16271/1‑148 (MQ=255)
                                                                        cgcAATGGATATTCTGTTCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaa                                                                        >  2:287640/1‑148 (MQ=255)
                                                                          cAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaaaa                                                                      <  2:2309247/148‑1 (MQ=255)
                                                                            aTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAcc                                                                    >  2:1244725/1‑148 (MQ=255)
                                                                            aTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAcc                                                                    >  2:169909/1‑148 (MQ=255)
                                                                            aTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAcc                                                                    >  2:3019886/1‑148 (MQ=255)
                                                                             tGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCt                                                                   >  2:1963491/1‑148 (MQ=255)
                                                                             tGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCt                                                                   >  2:162530/1‑148 (MQ=255)
                                                                              ggATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTg                                                                  >  1:1852347/1‑148 (MQ=255)
                                                                               gATATTCTGTGCACTGATAAAACTGGCACCCTGACACAGGACAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGc                                                                 >  1:2128557/1‑148 (MQ=255)
                                                                                atatTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCt                                                                >  2:806959/1‑148 (MQ=255)
                                                                                   ttCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGa                                                             >  1:1265192/1‑148 (MQ=255)
                                                                                        tgCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAg                                                        >  2:410970/1‑148 (MQ=255)
                                                                                        tgCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAg                                                        >  2:410885/1‑148 (MQ=255)
                                                                                          cACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCg                                                      <  2:37269/148‑1 (MQ=255)
                                                                                           aCTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCgg                                                     >  1:213512/1‑148 (MQ=255)
                                                                                              gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGANCGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                  <  1:1807377/148‑1 (MQ=255)
                                                                                              gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaaa                                                                       <  1:925426/127‑1 (MQ=255)
                                                                                              gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaaa                                                                       >  2:925365/1‑127 (MQ=255)
                                                                                              gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaaa                                                                       <  1:925365/127‑1 (MQ=255)
                                                                                              gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                  >  1:1332361/1‑148 (MQ=255)
                                                                                               aTAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCt                                                 >  2:1575960/1‑148 (MQ=255)
                                                                                                 aaaaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCg                                               >  1:177444/1‑148 (MQ=255)
                                                                                                   aaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGaa                                             >  1:2238667/1‑148 (MQ=255)
                                                                                                   aaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGaa                                             >  1:1785247/1‑148 (MQ=255)
                                                                                                    aCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                            <  1:1740713/148‑1 (MQ=255)
                                                                                                    aCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                            <  1:1548430/148‑1 (MQ=255)
                                                                                                    aCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                            <  1:3319586/148‑1 (MQ=255)
                                                                                                        gCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTAc                                        >  2:2322647/1‑148 (MQ=255)
                                                                                                         cACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACg                                       >  1:2244149/1‑148 (MQ=255)
                                                                                                          aCCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACgg                                      <  1:3134034/148‑1 (MQ=255)
                                                                                                           cccTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGa                                     <  2:2476774/148‑1 (MQ=255)
                                                                                                           cccTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGa                                     <  2:690111/148‑1 (MQ=255)
                                                                                                              tGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGa                                  >  2:926258/1‑148 (MQ=255)
                                                                                                              tGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGa                                  >  1:571545/1‑148 (MQ=255)
                                                                                                              tGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGa                                  >  1:571800/1‑148 (MQ=255)
                                                                                                               gACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGaa                                 <  2:1456912/148‑1 (MQ=255)
                                                                                                               gACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGaa                                 <  2:3228546/148‑1 (MQ=255)
                                                                                                                   cAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGt                             <  2:2934551/148‑1 (MQ=255)
                                                                                                                   cAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGt                             <  2:2935701/148‑1 (MQ=255)
                                                                                                                   cAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGt                             <  1:1350555/148‑1 (MQ=255)
                                                                                                                     ggATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCa                           >  1:228633/1‑148 (MQ=255)
                                                                                                                        tAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGATCGAAGGTACGGATGAAGAGTCAgcg                        >  1:3397607/1‑148 (MQ=255)
                                                                                                                         aaaaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATAAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGCATGAAGAGACAgcgc                       >  1:2832627/1‑148 (MQ=255)
                                                                                                                             ttGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctc                   >  1:1406666/1‑148 (MQ=255)
                                                                                                                              tgtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctcg                  >  2:2764621/1‑148 (MQ=255)
                                                                                                                               gtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctcgc                 <  2:824029/148‑1 (MQ=255)
                                                                                                                                    ggAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGcc            <  2:3371895/148‑1 (MQ=255)
                                                                                                                                     gagaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCa           >  1:3555315/1‑148 (MQ=255)
                                                                                                                                      agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAg          <  2:1406666/148‑1 (MQ=255)
                                                                                                                                      agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAg          >  2:1157015/1‑148 (MQ=255)
                                                                                                                                      agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAg          <  2:2631604/148‑1 (MQ=255)
                                                                                                                                       gaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGt         >  2:2230492/1‑148 (MQ=255)
                                                                                                                                           cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt     >  2:1431270/1‑148 (MQ=255)
                                                                                                                                           cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt     <  1:2638606/148‑1 (MQ=255)
                                                                                                                                           cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt     >  2:1760884/1‑148 (MQ=255)
                                                                                                                                           cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt     >  2:1760897/1‑148 (MQ=255)
                                                                                                                                              aCCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc  >  2:397556/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC  >  CP009273/4458464‑4458753

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: