Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 834,794 C→A 5.1% D732Y (GAT→TAT)  fiu ← catecholate siderophore receptor Fiu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273834,7940CA 5.1% 835.8 / 15.6 213D732Y (GAT→TAT) fiucatecholate siderophore receptor Fiu
Reads supporting (aligned to +/- strand):  ref base C (103/99);  new base A (7/4);  total (110/103)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 218 total aligned reads displayed.

CGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGATATAGC  >  CP009273/834645‑834943
                                                                                                                                                     |                                                                                                                                                     
 gggAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCg                                                                                                                                                      <  2:3532347/150‑1 (MQ=255)
  ggAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCgg                                                                                                                                                     >  2:3296111/1‑150 (MQ=255)
    aaTGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTa                                                                                                                                                   <  1:3387501/150‑1 (MQ=255)
       gtgGTACATCATACAAGTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCa                                                                                                                                                <  2:3646344/150‑1 (MQ=255)
        tgGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCaa                                                                                                                                               >  1:2275147/1‑150 (MQ=255)
          gTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAAc                                                                                                                                             <  1:1137322/150‑1 (MQ=255)
          gTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAAc                                                                                                                                             <  1:1749678/150‑1 (MQ=255)
            aCATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAg                                                                                                                                           <  1:6391/150‑1 (MQ=255)
            aCATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAg                                                                                                                                           <  1:3024678/150‑1 (MQ=255)
             catcatACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAgg                                                                                                                                          >  1:2363226/1‑150 (MQ=255)
                catACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTg                                                                                                                                       >  2:1255869/1‑150 (MQ=255)
                catACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTg                                                                                                                                       <  1:1030376/150‑1 (MQ=255)
                 atACACTGTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGt                                                                                                                                      <  2:3151475/150‑1 (MQ=255)
                 atACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGACATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCGCGCCCGGGTGATAACGGTAGCCGCTCTGGTTGTTTGAGGCGACGTAATCGGTATCAAACAGGTTGt                                                                                                                                      >  2:3411617/1‑150 (MQ=255)
                 atACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGt                                                                                                                                      >  2:2934470/1‑150 (MQ=255)
                   aCAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTaa                                                                                                                                    >  1:2740207/1‑150 (MQ=255)
                     aaTTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAAc                                                                                                                                  <  1:3190379/150‑1 (MQ=255)
                      aTTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACg                                                                                                                                 >  2:217270/1‑150 (MQ=255)
                        ttCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGtt                                                                                                                               <  2:740083/150‑1 (MQ=255)
                         tctcCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGttt                                                                                                                              <  1:1395290/150‑1 (MQ=255)
                           tcCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAg                                                                                                                            <  1:2318056/150‑1 (MQ=255)
                           tcCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAg                                                                                                                            <  1:2780509/150‑1 (MQ=255)
                             cAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCt                                                                                                                          >  2:3352821/1‑150 (MQ=255)
                              aaaaaGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTg                                                                                                                         >  1:3802253/1‑150 (MQ=255)
                                aaaGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGa                                                                                                                       <  1:1991276/150‑1 (MQ=255)
                                  aGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATAGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAg                                                                                                                     <  1:2449783/150‑1 (MQ=255)
                                   gTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGt                                                                                                                    <  1:1362049/150‑1 (MQ=255)
                                   gTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGt                                                                                                                    >  2:3118368/1‑150 (MQ=255)
                                     ggggCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCg                                                                                                                  <  2:697549/150‑1 (MQ=255)
                                     ggggCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCg                                                                                                                  >  2:1581448/1‑150 (MQ=255)
                                     ggggCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGATCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCg                                                                                                                  >  1:1823860/1‑150 (MQ=255)
                                        gCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCgaga                                                                                                               <  1:510353/150‑1 (MQ=255)
                                          cTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGAtt                                                                                                             <  1:982387/150‑1 (MQ=255)
                                          cTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGAtt                                                                                                             <  2:128421/150‑1 (MQ=255)
                                           tGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTg                                                                                                            <  2:1319115/150‑1 (MQ=255)
                                             cgcCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCg                                                                                                          <  2:1520651/150‑1 (MQ=255)
                                              gcCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTTGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGa                                                                                                         >  2:1550298/1‑150 (MQ=255)
                                               ccccACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGAt                                                                                                        <  1:3557302/150‑1 (MQ=255)
                                               ccccACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGAt                                                                                                        >  1:92525/1‑150 (MQ=255)
                                                         gAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATa                                                                                              >  1:2958597/1‑150 (MQ=255)
                                                          aaTCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATAc                                                                                             >  2:1054543/1‑150 (MQ=255)
                                                           aTCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATAcc                                                                                            >  2:664840/1‑150 (MQ=255)
                                                           aTCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATAcc                                                                                            <  1:1315437/150‑1 (MQ=255)
                                                              aGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCa                                                                                         >  1:910468/1‑150 (MQ=255)
                                                              aGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCa                                                                                         <  1:2152107/150‑1 (MQ=255)
                                                                 aaTGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGtt                                                                                      >  1:1612960/1‑150 (MQ=255)
                                                                 aaTGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGtt                                                                                      >  2:2884085/1‑150 (MQ=255)
                                                                 aaTGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATAGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGtt                                                                                      >  2:2039361/1‑150 (MQ=255)
                                                                   tGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTg                                                                                    <  1:2529/150‑1 (MQ=255)
                                                                     cATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGc                                                                                  >  1:2238830/1‑150 (MQ=255)
                                                                      atatTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCa                                                                                 <  2:2142324/150‑1 (MQ=255)
                                                                      atatTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCa                                                                                 <  2:498522/150‑1 (MQ=255)
                                                                      atatTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCa                                                                                 <  2:1809924/150‑1 (MQ=255)
                                                                       tatTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCAt                                                                                >  2:3244209/1‑150 (MQ=255)
                                                                            gCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACTGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCg                                                                           >  1:3488072/1‑150 (MQ=255)
                                                                            gCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCg                                                                           >  1:3577616/1‑150 (MQ=255)
                                                                             cTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGa                                                                          <  2:3112608/150‑1 (MQ=255)
                                                                                 gAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCa                                                                      <  2:2643293/150‑1 (MQ=255)
                                                                                 gAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCa                                                                      >  1:2595428/1‑150 (MQ=255)
                                                                                     aagaagGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTaa                                                                  >  2:1105926/1‑150 (MQ=255)
                                                                                       gaagGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAAcc                                                                >  2:350447/1‑150 (MQ=255)
                                                                                        aagGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCt                                                               >  2:2604753/1‑150 (MQ=255)
                                                                                         agGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCtt                                                              <  1:1715414/150‑1 (MQ=255)
                                                                                           gTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCg                                                            <  2:1642208/150‑1 (MQ=255)
                                                                                           gTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCg                                                            >  2:189610/1‑150 (MQ=255)
                                                                                                tGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTaaa                                                       >  1:3215456/1‑150 (MQ=255)
                                                                                                  gCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACg                                                     <  1:894800/150‑1 (MQ=255)
                                                                                                  gCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACg                                                     >  2:2693960/1‑150 (MQ=255)
                                                                                                    tCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCt                                                   >  1:4147738/1‑150 (MQ=255)
                                                                                                      gCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTgg                                                 <  1:2192201/150‑1 (MQ=255)
                                                                                                      gCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTgg                                                 <  2:3209839/150‑1 (MQ=255)
                                                                                                      gCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTgg                                                 >  1:965999/1‑150 (MQ=255)
                                                                                                             tGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTccc                                          <  1:1879921/150‑1 (MQ=255)
                                                                                                             tGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTccc                                          <  1:945935/150‑1 (MQ=255)
                                                                                                                tAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCAcc                                       >  1:1845169/1‑150 (MQ=255)
                                                                                                                tAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCAcc                                       <  2:3901065/150‑1 (MQ=255)
                                                                                                                tAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCAcc                                       <  1:57752/150‑1 (MQ=255)
                                                                                                                 aaCGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACcg                                      >  1:1499731/1‑150 (MQ=255)
                                                                                                                  aCGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACcgc                                     >  1:4170435/1‑150 (MQ=255)
                                                                                                                   cGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACcgcg                                    <  2:4126029/150‑1 (MQ=255)
                                                                                                                         ccGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCt                              <  2:373747/150‑1 (MQ=255)
                                                                                                                         ccGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCt                              <  1:3100929/150‑1 (MQ=255)
                                                                                                                            ctctTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGaa                           <  2:2116340/150‑1 (MQ=255)
                                                                                                                              ctTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAAcc                         >  1:1365667/1‑150 (MQ=255)
                                                                                                                                tgttgATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCtt                       <  1:522365/150‑1 (MQ=255)
                                                                                                                                   tgattgaGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTAt                    >  2:229702/1‑150 (MQ=255)
                                                                                                                                   tgattgaGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTAt                    <  1:1654042/150‑1 (MQ=255)
                                                                                                                                    gattgaGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATg                   <  2:1728391/150‑1 (MQ=255)
                                                                                                                                      ttgaGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCa                 <  2:2532414/150‑1 (MQ=255)
                                                                                                                                        gaGGCGACGTAATCTGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTAGAGATTGCGATTAACTCGATACCCCAGTTTGGAATCGGCGACCCAGTAAACTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACATTTATGCATa               >  1:180414/1‑150 (MQ=255)
                                                                                                                                          ggCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACt             <  2:143379/150‑1 (MQ=255)
                                                                                                                                            cGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGc           >  1:3882656/1‑150 (MQ=255)
                                                                                                                                               cGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGa        >  2:3108125/1‑150 (MQ=255)
                                                                                                                                                gTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGat       >  1:3480148/1‑150 (MQ=255)
                                                                                                                                                gTAATAGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGat       <  1:1190937/150‑1 (MQ=255)
                                                                                                                                                   aTCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGAtata    >  1:1498452/1‑150 (MQ=255)
                                                                                                                                                    tCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGATATAg   >  2:3250440/1‑150 (MQ=255)
                                                                                                                                                    tCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGATATAg   >  1:380466/1‑150 (MQ=255)
                                                                                                                                                     cGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGATATAGc  <  2:1597187/150‑1 (MQ=255)
                                                                                                                                                     cGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGATATAGc  <  2:1419135/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGAATGTGGTACATCATACAATTTCTCCAAAAAGTGGGGCCTGCGCCCCACATCTGAATCAGAAATGCATATTGGCTGTGAGCAAGAAGGTTCTTGGCTCGCCCGGGTGATAACGGTAGCCGCTCTTGTTGATTGAGGCGACGTAATCGGTATCAAACAGGTTGTAAACGTTTAGCTGGAAGTCGAGATTGCGATTAACTCGATACCCCAGTTTGGCATCGGCGACCCAGTAACCTTCGGTAAACGCTGGCGTTCCCACCGCGCCGTCTGAACCTTTATGCATACTGCCGATATAGC  >  CP009273/834645‑834943

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: ‑