Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,504,603 G→T 6.3% T631K (ACG→AAG)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,504,6030GT 6.3% 1421.1 / 43.1 363T631K (ACG→AAG) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (167/173);  new base T (10/13);  total (177/186)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 369 total aligned reads displayed.

TCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCC  >  CP009273/4504454‑4504751
                                                                                                                                                     |                                                                                                                                                    
tCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCg                                                                                                                                                      >  2:935831/1‑150 (MQ=255)
  cAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGtt                                                                                                                                                    <  2:2800158/150‑1 (MQ=255)
  cAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGtt                                                                                                                                                    <  2:1577779/150‑1 (MQ=255)
   aGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTa                                                                                                                                                   >  2:2054945/1‑150 (MQ=255)
     tttCCTGGCTTGTAGTCCACGCCGAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCTCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGc                                                                                                                                                 <  2:3986600/150‑1 (MQ=255)
     tttCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGc                                                                                                                                                 >  1:3475544/1‑150 (MQ=255)
          tGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTAcc                                                                                                                                            <  2:2742860/150‑1 (MQ=255)
          tGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTAcc                                                                                                                                            <  2:4679519/150‑1 (MQ=255)
           ggCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTAccc                                                                                                                                           <  2:4426866/150‑1 (MQ=255)
            gCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCa                                                                                                                                          >  1:2381749/1‑150 (MQ=255)
             cTTGTAGTCCACGCCCAGCGTGCCTGTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAg                                                                                                                                         >  2:3117837/1‑150 (MQ=255)
              ttGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGa                                                                                                                                        <  2:2558166/150‑1 (MQ=255)
               tGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGAt                                                                                                                                       >  1:693574/1‑150 (MQ=255)
                gTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATc                                                                                                                                      <  2:4702489/150‑1 (MQ=255)
                  aGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGt                                                                                                                                    >  2:2615382/1‑150 (MQ=255)
                    tCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTaa                                                                                                                                  <  2:692668/150‑1 (MQ=255)
                     ccACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAAc                                                                                                                                 >  1:1330374/1‑150 (MQ=255)
                       aCGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGt                                                                                                                               <  2:2176206/150‑1 (MQ=255)
                           ccAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCt                                                                                                                           <  2:2164105/150‑1 (MQ=255)
                           ccAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCt                                                                                                                           <  2:4934165/150‑1 (MQ=255)
                           ccAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCt                                                                                                                           <  2:1370739/150‑1 (MQ=255)
                           ccAGCGTGCCTTTATGGTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCt                                                                                                                           <  2:2163808/150‑1 (MQ=255)
                             aGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGc                                                                                                                         >  1:684826/1‑150 (MQ=255)
                             aGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGc                                                                                                                         >  2:828456/1‑150 (MQ=255)
                             aGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGc                                                                                                                         >  1:1489075/1‑150 (MQ=255)
                               cGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGCGt                                                                                                                       >  1:62155/1‑150 (MQ=255)
                                gTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGtt                                                                                                                      <  1:3803976/150‑1 (MQ=255)
                                 tGCCTTTCTGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGttt                                                                                                                     <  2:3408499/150‑1 (MQ=255)
                                 tGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGttt                                                                                                                     <  1:4489219/150‑1 (MQ=255)
                                   ccTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTcc                                                                                                                   <  1:2159275/150‑1 (MQ=255)
                                   ccTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTcc                                                                                                                   >  1:1632795/1‑150 (MQ=255)
                                     tttATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAg                                                                                                                 >  1:4804150/1‑150 (MQ=255)
                                            ttCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTa                                                                                                          >  2:483754/1‑150 (MQ=255)
                                              cGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATg                                                                                                        <  1:2702765/150‑1 (MQ=255)
                                              cGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATg                                                                                                        <  2:4306065/150‑1 (MQ=255)
                                              cGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATg                                                                                                        >  1:3358778/1‑150 (MQ=255)
                                               ggggAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATgg                                                                                                       >  1:4366700/1‑150 (MQ=255)
                                               ggggAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATgg                                                                                                       <  2:4765841/150‑1 (MQ=255)
                                                 ggAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGgcg                                                                                                     >  2:1912656/1‑150 (MQ=255)
                                                  gagaATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGgcgc                                                                                                    >  2:1444554/1‑150 (MQ=255)
                                                       tGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGtttt                                                                                               >  1:4113132/1‑150 (MQ=255)
                                                       tGGTACCAGATTGCCGTAGGTGTCGACTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGAGtttt                                                                                               >  2:4872891/1‑150 (MQ=255)
                                                        ggTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTg                                                                                              <  1:1810534/150‑1 (MQ=255)
                                                            ccAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCAc                                                                                          <  2:2945996/150‑1 (MQ=255)
                                                              aGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGt                                                                                        >  1:4867382/1‑150 (MQ=255)
                                                              aGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGt                                                                                        <  1:4330069/150‑1 (MQ=255)
                                                                 ttGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCa                                                                                     >  2:3669481/1‑150 (MQ=255)
                                                                 ttGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCa                                                                                     >  1:4105057/1‑150 (MQ=255)
                                                                   gCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGt                                                                                   <  1:3180480/150‑1 (MQ=255)
                                                                       tAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACg                                                                               <  1:3344794/150‑1 (MQ=255)
                                                                        aGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACgg                                                                              >  1:3276207/1‑150 (MQ=255)
                                                                        aGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACgg                                                                              <  1:4587753/150‑1 (MQ=255)
                                                                          gtgtCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGgtg                                                                            <  2:3631419/150‑1 (MQ=255)
                                                                           tgtCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGgtgt                                                                           <  2:3907919/150‑1 (MQ=255)
                                                                            gtCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTc                                                                          <  2:3212126/150‑1 (MQ=255)
                                                                            gtCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTc                                                                          <  2:3925255/150‑1 (MQ=255)
                                                                              cGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGt                                                                        >  2:3453086/1‑150 (MQ=255)
                                                                              cGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGt                                                                        >  2:178626/1‑150 (MQ=255)
                                                                              cGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTCCGGTTCCAGCCCGGTATGGCGCGTTTTGCCACCTGCAGAGACGGTGTCGt                                                                        >  1:3882883/1‑150 (MQ=255)
                                                                               gCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGtt                                                                       >  1:3945745/1‑150 (MQ=255)
                                                                                 cTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTgg                                                                     >  1:2279044/1‑150 (MQ=255)
                                                                                   ttCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTc                                                                   >  2:4404608/1‑150 (MQ=255)
                                                                                     ctcGCGGATTTCCGCGTTCACAGACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTg                                                                 <  2:450771/150‑1 (MQ=255)
                                                                                          ggATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTgg                                                            >  2:4564151/1‑150 (MQ=255)
                                                                                           gATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTTGCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGa                                                           <  1:3639192/150‑1 (MQ=255)
                                                                                                 cgcgTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATAAAGCGTAGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTa                                                     >  2:5114813/1‑150 (MQ=255)
                                                                                                     ttCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGa                                                 >  2:3044009/1‑150 (MQ=255)
                                                                                                     ttCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGa                                                 >  2:804661/1‑150 (MQ=255)
                                                                                                     ttCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGa                                                 >  2:4094942/1‑150 (MQ=255)
                                                                                                        acaTACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTg                                              >  1:104042/1‑150 (MQ=255)
                                                                                                         caTACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGt                                             >  2:1195765/1‑150 (MQ=255)
                                                                                                         caTACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGt                                             >  2:3767514/1‑150 (MQ=255)
                                                                                                              gCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAg                                        >  2:3362637/1‑150 (MQ=255)
                                                                                                               cATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGt                                       <  1:854628/150‑1 (MQ=255)
                                                                                                                aTAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGtt                                      <  1:4983809/150‑1 (MQ=255)
                                                                                                                  aGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTaa                                    <  2:1328014/150‑1 (MQ=255)
                                                                                                                   gCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTaac                                   >  1:1899805/1‑149 (MQ=255)
                                                                                                                    cTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATc                                  >  1:1989987/1‑150 (MQ=255)
                                                                                                                      ggCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAg                                <  2:4953130/150‑1 (MQ=255)
                                                                                                                       gCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAgg                               >  2:2628632/1‑150 (MQ=255)
                                                                                                                       gCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAgg                               >  2:25411/1‑150 (MQ=255)
                                                                                                                       gCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAgg                               <  1:123555/150‑1 (MQ=255)
                                                                                                                         gTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGaa                             >  1:2676737/1‑150 (MQ=255)
                                                                                                                         gTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGaa                             <  2:1632795/150‑1 (MQ=255)
                                                                                                                           aGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACa                           >  1:847862/1‑150 (MQ=255)
                                                                                                                             aTGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGc                         >  2:4060247/1‑150 (MQ=255)
                                                                                                                             aTGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGc                         <  1:1214461/150‑1 (MQ=255)
                                                                                                                                gAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGcccc                      <  1:946336/150‑1 (MQ=255)
                                                                                                                                 aaaCGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCa                     >  2:2949648/1‑150 (MQ=255)
                                                                                                                                   aCGTTATCAAGCGTTGGCTTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCAtt                   <  1:1206386/150‑1 (MQ=255)
                                                                                                                                       tATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTccg               <  1:5090640/150‑1 (MQ=255)
                                                                                                                                        aTCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTccgc              >  2:2574464/1‑150 (MQ=255)
                                                                                                                                          cAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTccgccg            >  2:2856515/1‑150 (MQ=255)
                                                                                                                                           aaGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGt           >  2:4892364/1‑150 (MQ=255)
                                                                                                                                             gCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCa         >  2:1684446/1‑150 (MQ=255)
                                                                                                                                             gCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCa         <  1:1058493/150‑1 (MQ=255)
                                                                                                                                              cGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAg        >  1:2986752/1‑150 (MQ=255)
                                                                                                                                              cGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAg        >  1:2570702/1‑150 (MQ=255)
                                                                                                                                                    cGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGcc  <  2:3271767/150‑1 (MQ=255)
                                                                                                                                                    cGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGcc  >  1:1745333/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCC  >  CP009273/4504454‑4504751

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: