Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,535,723 G→T 10.9% D271Y (GAT→TAT)  fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,535,7230GT10.9% 880.8 / 53.3 246D271Y (GAT→TAT) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base G (116/103);  new base T (6/21);  total (122/124)
Fisher's exact test for biased strand distribution p-value = 3.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 250 total aligned reads displayed.

GTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATT  >  CP009273/4535574‑4535872
                                                                                                                                                     |                                                                                                                                                     
gTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCg                                                                                                                                                       >  1:1269679/1‑150 (MQ=255)
 tttACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGa                                                                                                                                                      <  2:4494772/150‑1 (MQ=255)
     cgcgACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATAtt                                                                                                                                                  <  2:4012274/150‑1 (MQ=255)
      gcgACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCTATAttt                                                                                                                                                 <  2:2882663/150‑1 (MQ=255)
       cgACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATAtttt                                                                                                                                                >  1:1171664/1‑150 (MQ=255)
         aCAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCg                                                                                                                                              <  1:3932110/150‑1 (MQ=255)
           aaTACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGAt                                                                                                                                            >  2:2949449/1‑150 (MQ=255)
            aTACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATg                                                                                                                                           <  2:299914/150‑1 (MQ=255)
                  ccTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTa                                                                                                                                     <  1:1581252/150‑1 (MQ=255)
                   cTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGATCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTaa                                                                                                                                    <  2:2467108/150‑1 (MQ=255)
                   cTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTaa                                                                                                                                    <  2:902978/150‑1 (MQ=255)
                     ggAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACt                                                                                                                                  >  2:132633/1‑150 (MQ=255)
                     ggAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACt                                                                                                                                  >  1:1922778/1‑150 (MQ=255)
                     ggAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACt                                                                                                                                  >  1:2735111/1‑150 (MQ=255)
                     ggAGTTATAACAGAAGCGACATATCATCAGGTAGCAAAAAAAAATGCCAGCAAATCAATACCTCGCTTCAGAGAGACAAAAAACCGATACGTTCCCGGCTGACGCTGGGAGATGGTAATACTCAGGGCGAAATTTTCGATGGTATTAAca                                                                                                                                  >  2:3619113/1‑149 (MQ=255)
                       aGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACttt                                                                                                                                <  1:4210851/150‑1 (MQ=255)
                         ttATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTcg                                                                                                                              >  1:3929090/1‑150 (MQ=255)
                           ataACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTcgcg                                                                                                                            <  2:3893166/150‑1 (MQ=255)
                                   aGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACaa                                                                                                                    >  2:2107667/1‑150 (MQ=255)
                                   aGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACaa                                                                                                                    >  1:1531540/1‑150 (MQ=255)
                                     cGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAAtt                                                                                                                  >  2:4404884/1‑150 (MQ=255)
                                     cGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAAtt                                                                                                                  >  2:2077920/1‑150 (MQ=255)
                                      gACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTg                                                                                                                 >  2:4360533/1‑150 (MQ=255)
                                      gACAGATCATCAGGTAGAAAAAATAAATGGCATCATATCAATTCCTGGCTTGAGCGAGACATAATACCGTTACGTTCCAGGCTGACGCTGGTTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTg                                                                                                                 <  2:1443899/150‑1 (MQ=255)
                                       aCAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTgg                                                                                                                <  2:556272/150‑1 (MQ=255)
                                       aCAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTgg                                                                                                                >  2:3414824/1‑150 (MQ=255)
                                        cAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGc                                                                                                               >  2:1430570/1‑150 (MQ=255)
                                          gATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCt                                                                                                             >  2:2493654/1‑150 (MQ=255)
                                           atcatcAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTc                                                                                                            <  2:784984/150‑1 (MQ=255)
                                              atcaGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGa                                                                                                         >  1:3030732/1‑150 (MQ=255)
                                               tcaGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCATCAGAt                                                                                                        >  1:3343491/1‑150 (MQ=255)
                                                 aGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGa                                                                                                      >  2:3213732/1‑150 (MQ=255)
                                                 aGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGa                                                                                                      >  2:21507/1‑150 (MQ=255)
                                                  ggTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGAc                                                                                                     <  2:4092379/150‑1 (MQ=255)
                                                   gTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACa                                                                                                    <  2:470281/150‑1 (MQ=255)
                                                     aGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAat                                                                                                  >  2:1279633/1‑150 (MQ=255)
                                                       cAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAatat                                                                                                >  1:3525087/1‑150 (MQ=255)
                                                       cAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAatat                                                                                                <  1:2686308/150‑1 (MQ=255)
                                                        aaaaaTAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATg                                                                                               >  2:2624247/1‑150 (MQ=255)
                                                        aaaaaTAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATg                                                                                               >  2:4275395/1‑150 (MQ=255)
                                                        aaaaaTAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATg                                                                                               <  1:2165843/150‑1 (MQ=255)
                                                        aaaaaTAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATg                                                                                               >  1:2078388/1‑150 (MQ=255)
                                                         aaaaTAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGt                                                                                              <  1:4502484/150‑1 (MQ=255)
                                                              aaatGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTAccc                                                                                         >  1:1397514/1‑150 (MQ=255)
                                                              aaatGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTAccc                                                                                         <  2:4500887/150‑1 (MQ=255)
                                                               aatGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCg                                                                                        <  2:2724235/150‑1 (MQ=255)
                                                               aatGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCg                                                                                        >  1:5670/1‑150 (MQ=255)
                                                               aatGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGACTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCg                                                                                        <  2:3756171/150‑1 (MQ=255)
                                                                 tGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGAt                                                                                      >  1:3164930/1‑150 (MQ=255)
                                                                   gcagcaTATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAg                                                                                    <  2:2514813/150‑1 (MQ=255)
                                                                    cagcaTATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGt                                                                                   >  1:4317686/1‑150 (MQ=255)
                                                                     agcaTATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTc                                                                                  >  1:2549875/1‑150 (MQ=255)
                                                                      gcaTATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGGTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCa                                                                                 <  1:1031567/150‑1 (MQ=255)
                                                                       caTATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGAGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCaa                                                                                >  2:105494/1‑150 (MQ=255)
                                                                          atCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAaga                                                                             >  1:2003910/1‑150 (MQ=255)
                                                                          atCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAaga                                                                             >  1:4037543/1‑150 (MQ=255)
                                                                            cAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAgg                                                                           >  2:2836170/1‑150 (MQ=255)
                                                                              aTACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGAt                                                                         >  1:3792613/1‑150 (MQ=255)
                                                                                aCCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGAttt                                                                       >  1:4605198/1‑150 (MQ=255)
                                                                                 ccTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTg                                                                      <  1:2480728/150‑1 (MQ=255)
                                                                                 ccTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTg                                                                      <  1:3267597/150‑1 (MQ=255)
                                                                                  cTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGc                                                                     >  1:940250/1‑150 (MQ=255)
                                                                                   tGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGcc                                                                    >  2:3586328/1‑150 (MQ=255)
                                                                                           gCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGAt                                                            >  1:555834/1‑150 (MQ=255)
                                                                                           gCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGAt                                                            >  1:4384926/1‑150 (MQ=255)
                                                                                            cGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCACAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATc                                                           >  2:4189303/1‑150 (MQ=255)
                                                                                                  ataataCCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGt                                                     >  1:2286460/1‑150 (MQ=255)
                                                                                                  ataataCCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGt                                                     <  1:3165065/150‑1 (MQ=255)
                                                                                                    aataCCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTAt                                                   >  2:1140713/1‑150 (MQ=255)
                                                                                                     ataCCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTAtt                                                  <  1:4565381/150‑1 (MQ=255)
                                                                                                     ataCCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTAtt                                                  <  1:3118618/150‑1 (MQ=255)
                                                                                                          gTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCg                                             >  2:2340423/1‑150 (MQ=255)
                                                                                                               gTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTa                                        <  1:2866316/150‑1 (MQ=255)
                                                                                                                    cGGCTGACGCTGGGTGATGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGca                                   <  2:1092749/150‑1 (MQ=255)
                                                                                                                     ggCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGcac                                  >  2:3304485/1‑150 (MQ=255)
                                                                                                                      gCTGACGCTGGGTGATGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGcaca                                 >  1:852468/1‑150 (MQ=255)
                                                                                                                        tGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAgg                               >  1:3401848/1‑150 (MQ=255)
                                                                                                                         gACGCTGGGTGATGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGt                              <  2:3346819/150‑1 (MQ=255)
                                                                                                                           cGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCa                            >  1:1244937/1‑150 (MQ=255)
                                                                                                                               gggTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTAt                        <  2:3084760/150‑1 (MQ=255)
                                                                                                                                 gTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCAATATTa                      <  1:259979/150‑1 (MQ=255)
                                                                                                                                    aTGGTTATACTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAAc                   <  1:1827424/150‑1 (MQ=255)
                                                                                                                                    aTGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAAc                   <  1:2349387/150‑1 (MQ=255)
                                                                                                                                      ggTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACaa                 >  1:3810444/1‑150 (MQ=255)
                                                                                                                                      ggTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACaa                 >  2:2769414/1‑150 (MQ=255)
                                                                                                                                      ggTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACaa                 >  1:1297067/1‑150 (MQ=255)
                                                                                                                                       gTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACaaa                <  1:176396/150‑1 (MQ=255)
                                                                                                                                         tataCTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAAt              >  1:4634621/1‑150 (MQ=255)
                                                                                                                                           taCTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACACTATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATgg            <  2:1408170/150‑1 (MQ=255)
                                                                                                                                            aCTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGAACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACTGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATggg           >  2:1946186/1‑150 (MQ=255)
                                                                                                                                             cTCAGGGCTATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGt          <  2:3126885/150‑1 (MQ=255)
                                                                                                                                             cTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGt          >  1:3122108/1‑150 (MQ=255)
                                                                                                                                             cTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGt          <  2:666988/150‑1 (MQ=255)
                                                                                                                                              tCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTa         <  2:173614/150‑1 (MQ=255)
                                                                                                                                               cAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTAt        <  1:3772720/150‑1 (MQ=255)
                                                                                                                                               cAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTAt        >  2:319408/1‑150 (MQ=255)
                                                                                                                                                aGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATg       <  2:572593/150‑1 (MQ=255)
                                                                                                                                                   gCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACa    <  1:260642/150‑1 (MQ=255)
                                                                                                                                                   gCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACa    <  1:3308477/150‑1 (MQ=255)
                                                                                                                                                    cGATATTTTCGAGGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACAt   >  2:3165157/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTCGATGGTATTAACTTTCGCGGCGCACAATTGGCCTCAGATGACAATATGTTACCCGATAGTCAAAGAGGATTTGCCCCGGTGATCCACGGTATTGCTCGTGGTACTGCACAGGTCACTATTAAACAAAATGGGTATGACATT  >  CP009273/4535574‑4535872

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: