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,6090GT96.0% 1123.0 / 29.2 328D390Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (3/10);  new base T (167/148);  total (170/158)
Fisher's exact test for biased strand distribution p-value = 4.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 332 total aligned reads displayed.

CTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAA  >  CP009273/4458462‑4458756
                                                                                                                                                   |                                                                                                                                                   
cTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTat                                                                                                                                                   <  2:944677/150‑1 (MQ=255)
cTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTat                                                                                                                                                   >  1:4158710/1‑150 (MQ=255)
 tGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCtata                                                                                                                                                  <  2:2651296/150‑1 (MQ=255)
 tGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCtata                                                                                                                                                  <  2:1104215/150‑1 (MQ=255)
    cgcgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGTACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATAtct                                                                                                                                               <  2:4682825/150‑1 (MQ=255)
     gcgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATAtctc                                                                                                                                              <  1:3820502/150‑1 (MQ=255)
       gcggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCg                                                                                                                                            <  1:2016413/150‑1 (MQ=255)
       gcggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCg                                                                                                                                            <  2:1937447/150‑1 (MQ=255)
            gcAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaaa                                                                                                                                       >  1:1992801/1‑150 (MQ=255)
            gcAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaaa                                                                                                                                       <  2:2277633/150‑1 (MQ=255)
            gcAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaaa                                                                                                                                       >  1:657575/1‑150 (MQ=255)
             cAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTaaaa                                                                                                                                      >  1:2482021/1‑150 (MQ=255)
              aGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAc                                                                                                                                     <  1:4071883/150‑1 (MQ=255)
              aGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAc                                                                                                                                     <  1:3178650/150‑1 (MQ=255)
               gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAcc                                                                                                                                    >  1:975579/1‑150 (MQ=255)
               gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAcc                                                                                                                                    >  2:1747045/1‑150 (MQ=255)
               gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAAcc                                                                                                                                    <  2:3588635/150‑1 (MQ=255)
               gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAAcc                                                                                                                                    <  2:2255996/150‑1 (MQ=255)
                 aaaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAg                                                                                                                                  <  1:4648303/150‑1 (MQ=255)
                  aaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGc                                                                                                                                 >  1:4316893/1‑150 (MQ=255)
                   aGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCg                                                                                                                                >  2:594781/1‑150 (MQ=255)
                          gAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGt                                                                                                                         >  2:1219325/1‑150 (MQ=255)
                                gAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCa                                                                                                                   >  2:2161611/1‑150 (MQ=255)
                                gAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCa                                                                                                                   >  1:4617633/1‑150 (MQ=255)
                                 aaaGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCAt                                                                                                                  >  2:262820/1‑150 (MQ=255)
                                 aaaGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCAGGTAAAACCAGCGAACGCGTGATGCAt                                                                                                                  >  2:259860/1‑150 (MQ=255)
                                       aTCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGcgcg                                                                                                            <  1:2246509/150‑1 (MQ=255)
                                           tCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGt                                                                                                        <  1:3058765/150‑1 (MQ=255)
                                            cAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGtt                                                                                                       <  1:2134766/150‑1 (MQ=255)
                                              aaCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGa                                                                                                     <  1:2616682/150‑1 (MQ=255)
                                                cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                   <  1:4650899/150‑1 (MQ=255)
                                                 aTCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACa                                                                                                  >  1:2785348/1‑150 (MQ=255)
                                                   cTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                                >  2:3635479/1‑150 (MQ=255)
                                                   cTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                                >  1:1185452/1‑150 (MQ=255)
                                                    tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                               <  1:2285080/150‑1 (MQ=255)
                                                    tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                               >  2:507650/1‑150 (MQ=255)
                                                       aTGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCAtt                                                                                            <  2:3141258/150‑1 (MQ=255)
                                                        tGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTa                                                                                           <  2:2896450/150‑1 (MQ=255)
                                                          cTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATc                                                                                         >  1:2712890/1‑150 (MQ=255)
                                                           tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCa                                                                                        >  2:1010663/1‑150 (MQ=255)
                                                            aTTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAg                                                                                       <  1:2283794/150‑1 (MQ=255)
                                                            aTTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAg                                                                                       >  1:3367488/1‑150 (MQ=255)
                                                            aTTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAg                                                                                       >  2:2998255/1‑150 (MQ=255)
                                                            aTTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAg                                                                                       >  1:2421732/1‑150 (MQ=255)
                                                               cAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGAcc                                                                                    >  2:2934605/1‑150 (MQ=255)
                                                                     tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACtt                                                                              <  2:2031986/150‑1 (MQ=255)
                                                                     tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACtt                                                                              <  1:3266849/150‑1 (MQ=255)
                                                                      ttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTa                                                                             >  1:4708443/1‑150 (MQ=255)
                                                                      ttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTa                                                                             >  1:354923/1‑150 (MQ=255)
                                                                        ggCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaa                                                                           >  1:2539988/1‑150 (MQ=255)
                                                                          cgcAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaaaa                                                                         >  1:4652291/1‑150 (MQ=255)
                                                                           gcAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAc                                                                        >  2:4134189/1‑150 (MQ=255)
                                                                           gcAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAc                                                                        <  1:4780429/150‑1 (MQ=255)
                                                                           gcAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAc                                                                        >  1:3454738/1‑150 (MQ=255)
                                                                            cAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAcc                                                                       >  2:2456995/1‑150 (MQ=255)
                                                                             aaTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCt                                                                      >  2:813673/1‑150 (MQ=255)
                                                                              aTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTg                                                                     >  1:233705/1‑150 (MQ=255)
                                                                                 gATATTCTGTGCACTGATAAAACCGGCACCCTGACTCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTc                                                                  >  1:1645956/1‑150 (MQ=255)
                                                                                 gATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTc                                                                  <  1:1837832/150‑1 (MQ=255)
                                                                                  atatTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCg                                                                 >  2:1819575/1‑150 (MQ=255)
                                                                                    atTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGAt                                                               <  1:3711144/150‑1 (MQ=255)
                                                                                          tgCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCg                                                         >  1:1064552/1‑150 (MQ=255)
                                                                                            cACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGt                                                       <  2:393058/150‑1 (MQ=255)
                                                                                                gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTc                                                   <  1:4251106/150‑1 (MQ=255)
                                                                                                gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTc                                                   >  1:1330365/1‑150 (MQ=255)
                                                                                                   aaaaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGaa                                                <  1:3978889/150‑1 (MQ=255)
                                                                                                    aaaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                               >  2:927697/1‑150 (MQ=255)
                                                                                                     aaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAgg                                              <  2:1627939/150‑1 (MQ=255)
                                                                                                     aaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAgg                                              >  2:3119554/1‑150 (MQ=255)
                                                                                                           cACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGa                                        <  2:1391506/150‑1 (MQ=255)
                                                                                                               cTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGaa                                    >  2:438664/1‑150 (MQ=255)
                                                                                                                tGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAag                                   >  1:3776793/1‑150 (MQ=255)
                                                                                                                     cAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCa                              >  1:3903237/1‑150 (MQ=255)
                                                                                                                     cAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCa                              >  1:136071/1‑150 (MQ=255)
                                                                                                                       ggATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGGAGAGTCAgc                            >  2:1300204/1‑150 (MQ=255)
                                                                                                                       ggATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgc                            >  2:1071419/1‑150 (MQ=255)
                                                                                                                           aaaaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgcgcgc                        >  1:2929964/1‑150 (MQ=255)
                                                                                                                             aaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctc                      <  1:885248/150‑1 (MQ=255)
                                                                                                                             aaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctc                      <  1:1079850/150‑1 (MQ=255)
                                                                                                                             aaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctc                      >  1:3286124/1‑150 (MQ=255)
                                                                                                                             aaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctc                      <  1:2934605/150‑1 (MQ=255)
                                                                                                                             aaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGAGACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctc                      <  2:4287732/150‑1 (MQ=255)
                                                                                                                                tgtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctcgct                   <  2:2368800/150‑1 (MQ=255)
                                                                                                                                 gtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTg                  <  2:4158710/150‑1 (MQ=255)
                                                                                                                                 gtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTg                  <  1:1010632/150‑1 (MQ=255)
                                                                                                                                 gtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTg                  >  2:1266042/1‑150 (MQ=255)
                                                                                                                                   gCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGc                <  1:4490561/150‑1 (MQ=255)
                                                                                                                                    cTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGcc               >  1:1478672/1‑150 (MQ=255)
                                                                                                                                     tGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCa              <  2:483366/150‑1 (MQ=255)
                                                                                                                                     tGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCa              <  1:4399237/150‑1 (MQ=255)
                                                                                                                                     tGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCa              <  2:279388/150‑1 (MQ=255)
                                                                                                                                       gagaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGt            >  2:15051/1‑150 (MQ=255)
                                                                                                                                        agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTCCGGATGTAGTGTCAGCGCGCTCGCTGGCCAGTc           >  2:3194984/1‑150 (MQ=255)
                                                                                                                                        agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTc           >  1:3461016/1‑150 (MQ=255)
                                                                                                                                        agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCAGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTc           >  1:1744210/1‑150 (MQ=255)
                                                                                                                                         gaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCg          >  1:282891/1‑150 (MQ=255)
                                                                                                                                          aaTCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGt         >  2:305762/1‑150 (MQ=255)
                                                                                                                                              aTACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc     >  1:530428/1‑150 (MQ=255)
                                                                                                                                               tACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCa    >  1:2655833/1‑150 (MQ=255)
                                                                                                                                                 ccGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCaaa  <  2:187643/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAAA  >  CP009273/4458462‑4458756

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: