Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,902,413 G→A 22.1% P205S (CCT→TCT)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,902,4130GA22.1% 702.4 / 194.2 312P205S (CCT→TCT) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base G (140/103);  new base A (31/38);  total (171/141)
Fisher's exact test for biased strand distribution p-value = 7.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

CTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTT  >  CP009273/2902264‑2902562
                                                                                                                                                     |                                                                                                                                                     
cTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAgg                                                                                                                                                       <  2:2881337/150‑1 (MQ=255)
cTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAgg                                                                                                                                                       >  1:4150223/1‑150 (MQ=255)
  gagaTTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCt                                                                                                                                                     >  1:2088710/1‑150 (MQ=255)
     aTTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGa                                                                                                                                                  <  2:2754868/150‑1 (MQ=255)
      tttCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGAt                                                                                                                                                 <  1:62225/150‑1 (MQ=255)
         cAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATAcc                                                                                                                                              >  2:4115585/1‑150 (MQ=255)
            cAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGAt                                                                                                                                           <  1:4856386/150‑1 (MQ=255)
             aGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATg                                                                                                                                          <  2:3838336/150‑1 (MQ=255)
             aGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATg                                                                                                                                          <  2:3304853/150‑1 (MQ=255)
             aGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATg                                                                                                                                          >  1:986847/1‑150 (MQ=255)
             aGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATg                                                                                                                                          >  1:3910869/1‑150 (MQ=255)
               gCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGga                                                                                                                                        <  2:643057/150‑1 (MQ=255)
               gCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGga                                                                                                                                        <  1:2725301/150‑1 (MQ=255)
               gCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGga                                                                                                                                        <  2:3074349/150‑1 (MQ=255)
                cccGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGgag                                                                                                                                       <  1:1103635/150‑1 (MQ=255)
                cccGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGgag                                                                                                                                       <  1:865472/150‑1 (MQ=255)
                 ccGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGgaga                                                                                                                                      <  1:3254802/150‑1 (MQ=255)
                  cGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGagag                                                                                                                                     <  1:3632841/150‑1 (MQ=255)
                   ggAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGagc                                                                                                                                    >  1:4708130/1‑150 (MQ=255)
                   ggAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGagc                                                                                                                                    <  1:4227585/150‑1 (MQ=255)
                          tataGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCtt                                                                                                                             <  2:1099588/150‑1 (MQ=255)
                          tataGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCtt                                                                                                                             <  2:4765479/150‑1 (MQ=255)
                           ataGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCttt                                                                                                                            <  2:763482/150‑1 (MQ=255)
                             aGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTa                                                                                                                          <  2:4569267/150‑1 (MQ=255)
                              gATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTAc                                                                                                                         <  2:2257681/150‑1 (MQ=255)
                                 ggAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACaga                                                                                                                      >  2:1411689/1‑150 (MQ=255)
                                   aaTCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGt                                                                                                                    >  2:4018521/1‑150 (MQ=255)
                                    aTCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTg                                                                                                                   >  2:2178640/1‑150 (MQ=255)
                                    aTCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTg                                                                                                                   >  2:3356500/1‑150 (MQ=255)
                                        aCGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGa                                                                                                               <  1:4592620/150‑1 (MQ=255)
                                        aCGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGa                                                                                                               <  1:2140462/150‑1 (MQ=255)
                                        aCGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGa                                                                                                               >  2:4099518/1‑150 (MQ=255)
                                           tCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTc                                                                                                            >  1:4403136/1‑150 (MQ=255)
                                            cTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCg                                                                                                           >  1:3590080/1‑150 (MQ=255)
                                            cTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCg                                                                                                           >  1:1002245/1‑150 (MQ=255)
                                            cTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCg                                                                                                           <  1:492555/150‑1 (MQ=255)
                                                cAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCggttt                                                                                                       <  1:76575/150‑1 (MQ=255)
                                                cAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCggttt                                                                                                       >  2:127872/1‑150 (MQ=255)
                                                     aaaTAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCggtttggttt                                                                                                  >  2:3637796/1‑150 (MQ=255)
                                                      aaTAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGtttt                                                                                                 >  2:2481427/1‑150 (MQ=255)
                                                         aaCCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGAc                                                                                              >  1:1570358/1‑150 (MQ=255)
                                                              cTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCAc                                                                                         >  1:2342872/1‑150 (MQ=255)
                                                              cTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCAc                                                                                         >  2:158160/1‑150 (MQ=255)
                                                              cTTTTACCGGAACATTACAGAACAATGCAATCTTCGCACGTACGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGGTTGGTTTTGACTTCAc                                                                                         >  1:397316/1‑150 (MQ=255)
                                                                    ccGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTTGTTTTGACTTCACCAGACg                                                                                   >  2:2787279/1‑150 (MQ=255)
                                                                    ccGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACg                                                                                   <  1:2910190/150‑1 (MQ=255)
                                                                    ccGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACg                                                                                   >  2:4388044/1‑150 (MQ=255)
                                                                         aCATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATCCCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGGTTGGTTTTGACTTCACCAGACGCTGcc                                                                              >  2:577110/1‑150 (MQ=255)
                                                                          cATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCa                                                                             <  1:1215623/150‑1 (MQ=255)
                                                                          cATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCa                                                                             <  1:2504698/150‑1 (MQ=255)
                                                                           aTTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCAt                                                                            >  2:1881405/1‑150 (MQ=255)
                                                                           aTTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCAt                                                                            <  1:4409226/150‑1 (MQ=255)
                                                                             tACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGt                                                                          <  2:2124794/150‑1 (MQ=255)
                                                                              aCAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTa                                                                         <  1:3288475/150‑1 (MQ=255)
                                                                              aCAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTa                                                                         <  2:4758274/150‑1 (MQ=255)
                                                                               cAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTAc                                                                        >  2:4252593/1‑150 (MQ=255)
                                                                                aGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGGTTGGTTTTGACTTCACCAGACGCTGCCATGTACg                                                                       >  1:1376893/1‑150 (MQ=255)
                                                                                  aaCAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCCGACGCTGCCATGTACGGc                                                                     >  1:1358806/1‑150 (MQ=255)
                                                                                  aaCAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGc                                                                     >  2:2425603/1‑150 (MQ=255)
                                                                                   aCAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCa                                                                    >  2:2578656/1‑150 (MQ=255)
                                                                                       tGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAg                                                                <  1:1879268/150‑1 (MQ=255)
                                                                                        gCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGc                                                               >  2:348044/1‑150 (MQ=255)
                                                                                         cAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCg                                                              <  1:4360373/150‑1 (MQ=255)
                                                                                         cAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCg                                                              >  1:4277205/1‑150 (MQ=255)
                                                                                          aaTCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGt                                                             >  1:2215170/1‑150 (MQ=255)
                                                                                           aTCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTc                                                            >  1:4774599/1‑150 (MQ=255)
                                                                                            tCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCa                                                           >  1:1705064/1‑150 (MQ=255)
                                                                                              ttCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAgg                                                         <  1:3394254/150‑1 (MQ=255)
                                                                                               tCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGt                                                        <  2:23346/150‑1 (MQ=255)
                                                                                                cGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTg                                                       >  1:4145652/1‑150 (MQ=255)
                                                                                                cGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTg                                                       <  2:318126/150‑1 (MQ=255)
                                                                                                  cACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGTCTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCa                                                     >  2:4430412/1‑150 (MQ=255)
                                                                                                   aCGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCAt                                                    >  2:3592780/1‑150 (MQ=255)
                                                                                                    cgttcgttCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATa                                                   >  1:1480134/1‑150 (MQ=255)
                                                                                                       tcgttcgCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAAc                                                <  2:1552356/150‑1 (MQ=255)
                                                                                                        cgttcgCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCCGACGCTGCCATGTTCGGCCCCAGCGTCAGGTGCATAAACa                                               >  1:2868268/1‑150 (MQ=255)
                                                                                                         gttcgCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACag                                              >  1:3165618/1‑150 (MQ=255)
                                                                                                         gttcgCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACag                                              >  2:1802018/1‑150 (MQ=255)
                                                                                                              ccGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAgtgt                                         >  1:4101300/1‑150 (MQ=255)
                                                                                                               cGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGtgtg                                        >  2:1343613/1‑150 (MQ=255)
                                                                                                                 gAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCt                                      >  1:511291/1‑150 (MQ=255)
                                                                                                                 gAACAGCGCGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCt                                      <  2:2597224/150‑1 (MQ=255)
                                                                                                                   aCAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCa                                    <  1:1835058/150‑1 (MQ=255)
                                                                                                                    cAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCAc                                   <  2:1703347/150‑1 (MQ=255)
                                                                                                                     aGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACg                                  <  1:1610671/150‑1 (MQ=255)
                                                                                                                      gcgcGATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACgg                                 <  1:2623177/150‑1 (MQ=255)
                                                                                                                       cgcgATCTGAACGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGc                                <  1:2864032/150‑1 (MQ=255)
                                                                                                                        gcgATCTGAACGACAAATCAGGATGTCAGGCGGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTTACTTCAACAGACGCTGCCCTGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGcc                               >  1:929381/1‑150 (MQ=255)
                                                                                                                            tCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAAtt                           <  1:2986310/150‑1 (MQ=255)
                                                                                                                             cTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAAttt                          >  1:701960/1‑150 (MQ=255)
                                                                                                                              tGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTTACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAAAAGAGTGTGCTCCCGGCCAAatta                         >  1:4471004/1‑149 (MQ=255)
                                                                                                                               gAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCa                        <  1:1689543/150‑1 (MQ=255)
                                                                                                                                 aCGACAAATCAGGATGTCAGACTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAAc                      >  2:3425813/1‑150 (MQ=255)
                                                                                                                                  cGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCaaaa                     >  1:656453/1‑148 (MQ=255)
                                                                                                                                         tCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCt              >  1:1793490/1‑150 (MQ=255)
                                                                                                                                          cAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTg             <  2:4537962/150‑1 (MQ=255)
                                                                                                                                            ggATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTGTGGTTTTGACTTAACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCCATTTCAACAGCCATCTGGc           >  1:4079321/1‑150 (MQ=255)
                                                                                                                                             gATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCg          <  2:1177424/150‑1 (MQ=255)
                                                                                                                                             gATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCg          >  2:2342824/1‑150 (MQ=255)
                                                                                                                                             gATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCg          >  2:3524849/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGAGTGCTGAGTCGGTTTGGTTTTGACTTCACCAGACGCTGCCATGTACGGCACCAGCGTCAGGTGCATAAACAGAGTGTGCTCACGGCCAATTTCAACAGCCATCTGGCGAATCGCTT  >  CP009273/2902264‑2902562

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: