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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,6090GT71.9% 415.4 / 274.2 277D390Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (48/30);  new base T (99/99);  total (148/129)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 278 total aligned reads displayed.

CTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAA  >  CP009273/4458462‑4458758
                                                                                                                                                   |                                                                                                                                                     
 tGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGata                                                                                                                                                    >  1:3157112/1‑150 (MQ=255)
   gcgcgcGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATAtc                                                                                                                                                  <  2:377263/150‑1 (MQ=255)
        cggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCgg                                                                                                                                             >  1:1301905/1‑150 (MQ=255)
              aGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAc                                                                                                                                       >  1:1007359/1‑150 (MQ=255)
              aGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAc                                                                                                                                       >  1:4203297/1‑150 (MQ=255)
               gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAcc                                                                                                                                      <  1:2421683/150‑1 (MQ=255)
               gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAAcc                                                                                                                                      >  2:4201866/1‑150 (MQ=255)
                 aaaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAg                                                                                                                                    >  2:413352/1‑150 (MQ=255)
                   aGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCg                                                                                                                                  >  2:3374026/1‑150 (MQ=255)
                         cGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAAcgcg                                                                                                                            >  1:3010476/1‑150 (MQ=255)
                          gAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGt                                                                                                                           <  1:3151378/150‑1 (MQ=255)
                           aaaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtg                                                                                                                          >  2:145507/1‑150 (MQ=255)
                           aaaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGtg                                                                                                                          <  1:2745278/150‑1 (MQ=255)
                            aaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGtgc                                                                                                                         <  1:3178493/150‑1 (MQ=255)
                             aCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgct                                                                                                                        <  2:1746953/150‑1 (MQ=255)
                              cAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgctg                                                                                                                       >  1:1467716/1‑150 (MQ=255)
                               aGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGtgctgc                                                                                                                      >  2:596581/1‑150 (MQ=255)
                               aGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGtgctgc                                                                                                                      >  1:1801970/1‑150 (MQ=255)
                                gAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCa                                                                                                                     <  1:3941310/150‑1 (MQ=255)
                                gAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCa                                                                                                                     >  2:1396074/1‑150 (MQ=255)
                                 aaaGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCAt                                                                                                                    >  1:350171/1‑150 (MQ=255)
                                    gTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgc                                                                                                                 <  1:188627/150‑1 (MQ=255)
                                    gTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgc                                                                                                                 <  1:3872322/150‑1 (MQ=255)
                                    gTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgc                                                                                                                 >  1:697804/1‑150 (MQ=255)
                                     tGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgcg                                                                                                                <  2:3095317/150‑1 (MQ=255)
                                       aTCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGcgcg                                                                                                              <  1:2378724/150‑1 (MQ=255)
                                       aTCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGcgcg                                                                                                              >  2:2612316/1‑150 (MQ=255)
                                        tCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGt                                                                                                             >  1:4358236/1‑150 (MQ=255)
                                         cGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTg                                                                                                            <  2:2067219/150‑1 (MQ=255)
                                         cGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTg                                                                                                            <  1:4086953/150‑1 (MQ=255)
                                            cAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGtt                                                                                                         >  1:847628/1‑150 (MQ=255)
                                            cAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGtt                                                                                                         <  1:3907174/150‑1 (MQ=255)
                                            cAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGtt                                                                                                         <  2:4014134/150‑1 (MQ=255)
                                               aCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGaa                                                                                                      >  1:3035375/1‑150 (MQ=255)
                                                cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                     >  2:1572513/1‑150 (MQ=255)
                                                cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                     >  2:1732666/1‑150 (MQ=255)
                                                cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                     <  1:4866920/150‑1 (MQ=255)
                                                cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                     >  2:3224685/1‑150 (MQ=255)
                                                   cTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                                  >  2:1534281/1‑150 (MQ=255)
                                                    tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                                 <  1:789066/150‑1 (MQ=255)
                                                    tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                                 <  2:2923989/150‑1 (MQ=255)
                                                    tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                                 <  2:951507/150‑1 (MQ=255)
                                                    tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                                 >  2:3845019/1‑150 (MQ=255)
                                                        tGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTa                                                                                             <  2:4418534/150‑1 (MQ=255)
                                                        tGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCAATATCTCCGGTAAAACCAGCGAACCCGTGCTGCATAGCGCGTGGTTGAACAGTCATTa                                                                                             >  1:2660081/1‑150 (MQ=255)
                                                          cTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATc                                                                                           >  2:458602/1‑150 (MQ=255)
                                                           tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCa                                                                                          >  1:1147983/1‑150 (MQ=255)
                                                           tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCa                                                                                          <  2:3270793/150‑1 (MQ=255)
                                                           tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCa                                                                                          <  2:2372266/150‑1 (MQ=255)
                                                            aTTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAg                                                                                         <  2:471512/150‑1 (MQ=255)
                                                              tCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGAc                                                                                       <  1:2570025/150‑1 (MQ=255)
                                                                  aaCTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGa                                                                                   >  2:2573255/1‑150 (MQ=255)
                                                                     tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACtt                                                                                <  2:868656/150‑1 (MQ=255)
                                                                      ttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTa                                                                               >  1:3086960/1‑150 (MQ=255)
                                                                      ttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTa                                                                               <  2:1408810/150‑1 (MQ=255)
                                                                       tGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaa                                                                              >  1:1722179/1‑150 (MQ=255)
                                                                        ggCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaa                                                                             >  1:1842358/1‑150 (MQ=255)
                                                                           gcAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAc                                                                          >  1:1870051/1‑150 (MQ=255)
                                                                           gcAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAc                                                                          >  1:4474894/1‑150 (MQ=255)
                                                                            cAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAcc                                                                         >  1:2007448/1‑150 (MQ=255)
                                                                             aaTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAcca                                                                        >  1:2256380/1‑149 (MQ=255)
                                                                               tGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGc                                                                      >  1:1082713/1‑150 (MQ=255)
                                                                                ggATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCt                                                                     <  1:1389096/150‑1 (MQ=255)
                                                                                   tatTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGa                                                                  >  2:3540612/1‑150 (MQ=255)
                                                                                        tgtgCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAg                                                             >  2:29139/1‑150 (MQ=255)
                                                                                             aCTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTg                                                        <  2:34685/150‑1 (MQ=255)
                                                                                              cTGATAATACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                       >  2:2194107/1‑150 (MQ=255)
                                                                                              cTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                       <  1:3455078/150‑1 (MQ=255)
                                                                                               tGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCt                                                      <  2:3975248/150‑1 (MQ=255)
                                                                                                 aTAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCg                                                    >  2:4272673/1‑150 (MQ=255)
                                                                                                    aagccGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCAGCTCGATACAGCGGTGCTCGAAg                                                 >  2:2055090/4‑150 (MQ=255)
                                                                                                    aaaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCAGCTCGATACAGCGGTGCTCGAAg                                                 >  2:4111844/1‑150 (MQ=255)
                                                                                                       ccGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTa                                              <  1:3720154/150‑1 (MQ=255)
                                                                                                        cGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCAGCTCGATACAGCGGTGCTCGAAGGTAc                                             <  2:2267783/150‑1 (MQ=255)
                                                                                                         ggCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACg                                            >  1:4190180/1‑150 (MQ=255)
                                                                                                          gCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACgg                                           >  1:3361333/1‑150 (MQ=255)
                                                                                                          gCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACgg                                           <  1:255283/150‑1 (MQ=255)
                                                                                                          gCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACgg                                           >  1:919824/1‑150 (MQ=255)
                                                                                                             cccTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATg                                        >  1:2143757/1‑150 (MQ=255)
                                                                                                                    gCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCAGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTc                                 <  2:4368844/150‑1 (MQ=255)
                                                                                                                       ggATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgc                              <  1:3441644/150‑1 (MQ=255)
                                                                                                                       ggATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgc                              <  1:995483/150‑1 (MQ=255)
                                                                                                                          tAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTCCGGATGAAGAGTCAgcgcg                           >  2:808335/1‑150 (MQ=255)
                                                                                                                          tAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgcgcg                           >  2:2090613/1‑150 (MQ=255)
                                                                                                                            aaaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgct                         >  2:3027278/1‑150 (MQ=255)
                                                                                                                                 gtgCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTg                    <  2:1109530/150‑1 (MQ=255)
                                                                                                                                    cTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGcc                 >  1:3318564/1‑150 (MQ=255)
                                                                                                                                      ggAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAg               <  2:2793207/150‑1 (MQ=255)
                                                                                                                                        agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTc             <  1:1500774/150‑1 (MQ=255)
                                                                                                                                        agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTc             >  2:3279593/1‑150 (MQ=255)
                                                                                                                                        agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTc             <  1:4709946/150‑1 (MQ=255)
                                                                                                                                          aaTCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGt           >  2:1269878/1‑150 (MQ=255)
                                                                                                                                           aTCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt          <  2:488073/150‑1 (MQ=255)
                                                                                                                                             cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTgg        <  1:4158505/150‑1 (MQ=255)
                                                                                                                                             cATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCAGCTCGAGACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTgg        <  2:3157112/150‑1 (MQ=255)
                                                                                                                                              aTACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc       >  1:4292998/1‑150 (MQ=255)
                                                                                                                                              aTACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGACGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc       >  1:2795491/1‑150 (MQ=255)
                                                                                                                                              aTACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc       >  1:365388/1‑150 (MQ=255)
                                                                                                                                               tACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCa      <  1:582075/150‑1 (MQ=255)
                                                                                                                                                   gATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCAGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCaaaaa  >  2:3653201/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAAAA  >  CP009273/4458462‑4458758

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: