Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,530,002 C→A 18.1% intergenic (‑684/‑772) nanC ← / → fimB N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,530,0020CA18.1% 826.4 / 120.8 279intergenic (‑684/‑772)nanC/fimBN‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base C (127/101);  new base A (26/25);  total (153/126)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGA  >  CP009273/4529853‑4530149
                                                                                                                                                     |                                                                                                                                                   
  aGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAg                                                                                                                                                   <  2:4345071/150‑1 (MQ=255)
     cAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCt                                                                                                                                                >  1:3656768/1‑150 (MQ=255)
       gACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGa                                                                                                                                              >  1:1403236/1‑150 (MQ=255)
        aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAg                                                                                                                                             <  2:302647/150‑1 (MQ=255)
        aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAg                                                                                                                                             <  2:2413240/150‑1 (MQ=255)
        aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAg                                                                                                                                             <  1:149019/150‑1 (MQ=255)
         cGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAgt                                                                                                                                            >  1:4543260/1‑149 (MQ=255)
          gAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCa                                                                                                                                           >  2:111259/1‑150 (MQ=255)
          gAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCa                                                                                                                                           >  2:4440839/1‑150 (MQ=255)
          gAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCa                                                                                                                                           >  2:3399352/1‑150 (MQ=255)
            acacTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATa                                                                                                                                         >  1:1736331/1‑150 (MQ=255)
             cacTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATaa                                                                                                                                        >  1:4475252/1‑150 (MQ=255)
              acTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAt                                                                                                                                       >  1:2073973/1‑150 (MQ=255)
              acTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAt                                                                                                                                       >  2:2522792/1‑150 (MQ=255)
                 tGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATtct                                                                                                                                    <  2:1454286/150‑1 (MQ=255)
                 tGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATtct                                                                                                                                    <  1:2157075/150‑1 (MQ=255)
                 tGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATtct                                                                                                                                    <  1:5100305/150‑1 (MQ=255)
                  gCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATtctc                                                                                                                                   <  1:3416073/150‑1 (MQ=255)
                    ggtggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCtcat                                                                                                                                 >  1:1018904/1‑150 (MQ=255)
                     gtggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcatc                                                                                                                                <  1:3091592/150‑1 (MQ=255)
                     gtggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcatc                                                                                                                                >  1:1913007/1‑150 (MQ=255)
                      tggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcatca                                                                                                                               <  1:482189/150‑1 (MQ=255)
                       ggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTcatcat                                                                                                                              >  1:2501297/1‑150 (MQ=255)
                       ggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTcatcat                                                                                                                              >  2:3993144/1‑150 (MQ=255)
                        gtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATg                                                                                                                             >  1:4770984/1‑150 (MQ=255)
                              aaaaaCTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAatat                                                                                                                       <  1:2693058/150‑1 (MQ=255)
                                 aaCTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGtt                                                                                                                    >  1:2916637/1‑150 (MQ=255)
                                  aCTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGttt                                                                                                                   <  1:2610935/150‑1 (MQ=255)
                                  aCTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTAGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATACTCATCATGAAATATGttt                                                                                                                   >  1:1513293/1‑150 (MQ=255)
                                   cTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTc                                                                                                                  <  2:1733831/150‑1 (MQ=255)
                                   cTAAAGATCGTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAACTCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTc                                                                                                                  <  2:1250343/150‑1 (MQ=255)
                                    tAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTcc                                                                                                                 >  1:4068362/1‑150 (MQ=255)
                                     aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCt                                                                                                                >  1:5129740/1‑150 (MQ=255)
                                     aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCt                                                                                                                <  2:1373825/150‑1 (MQ=255)
                                      aaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTg                                                                                                               <  1:4298183/150‑1 (MQ=255)
                                      aaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGTATAATTCTCATCATGAAATATGTTTCCTg                                                                                                               <  1:4129616/150‑1 (MQ=255)
                                          tCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGttt                                                                                                           <  1:629012/150‑1 (MQ=255)
                                          tCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGttt                                                                                                           <  1:4348656/150‑1 (MQ=255)
                                           cTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTtg                                                                                                          >  1:2953761/1‑150 (MQ=255)
                                             tAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGTATAATTCTCATCATGAAATATGTTTCCTGGTTtgtg                                                                                                        >  2:1766556/1‑150 (MQ=255)
                                              aGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTgg                                                                                                       >  2:5041504/1‑150 (MQ=255)
                                               gTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGTATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGc                                                                                                      >  2:570414/1‑150 (MQ=255)
                                                 ttAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCtt                                                                                                    <  2:2679338/150‑1 (MQ=255)
                                                   aaCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGt                                                                                                  <  2:1551207/150‑1 (MQ=255)
                                                    aCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGATCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTa                                                                                                 >  1:499548/1‑150 (MQ=255)
                                                    aCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTa                                                                                                 >  2:678640/1‑150 (MQ=255)
                                                    aCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTa                                                                                                 >  1:574104/1‑150 (MQ=255)
                                                    aCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTa                                                                                                 >  2:4733981/1‑150 (MQ=255)
                                                     cTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTaa                                                                                                <  2:2916637/150‑1 (MQ=255)
                                                       aTTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACt                                                                                              <  2:3595021/150‑1 (MQ=255)
                                                           gTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTc                                                                                          >  2:2224524/1‑150 (MQ=255)
                                                            ttttATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCa                                                                                         >  1:2554295/1‑150 (MQ=255)
                                                                 taaataaTTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCt                                                                                    <  2:2132032/150‑1 (MQ=255)
                                                                   antaatTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGa                                                                                  >  1:333165/1‑150 (MQ=255)
                                                                   aataatTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGa                                                                                  >  1:4976322/1‑150 (MQ=255)
                                                                   aataatTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGa                                                                                  >  1:333133/1‑150 (MQ=255)
                                                                     taatTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAg                                                                                >  2:2877305/1‑150 (MQ=255)
                                                                       atTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTc                                                                              >  2:1719360/1‑150 (MQ=255)
                                                                       atTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTc                                                                              >  1:2872476/1‑150 (MQ=255)
                                                                        tttATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCg                                                                             <  2:2420049/150‑1 (MQ=255)
                                                                           attaAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATc                                                                          <  2:3852140/150‑1 (MQ=255)
                                                                           attaAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATc                                                                          >  1:147039/1‑150 (MQ=255)
                                                                           attaAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATc                                                                          >  2:3706497/1‑150 (MQ=255)
                                                                                gagTCTAAACAAGGGGAGCTTTGCAAGCAAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGCTATGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGgag                                                                     >  2:4062209/1‑150 (MQ=255)
                                                                                   tttAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAgg                                                                  >  2:2835818/3‑150 (MQ=255)
                                                                                      aaaCAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCtt                                                               >  1:4241964/1‑150 (MQ=255)
                                                                                      aaaCAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCtt                                                               <  2:3057829/150‑1 (MQ=255)
                                                                                      aaaCAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCtt                                                               <  1:3733059/150‑1 (MQ=255)
                                                                                         cAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCttgtt                                                            <  2:1785830/150‑1 (MQ=255)
                                                                                             gggAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATg                                                        >  2:4825193/1‑150 (MQ=255)
                                                                                                    ttGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGgtgt                                                 >  2:4260677/1‑150 (MQ=255)
                                                                                                     tGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGAATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTAACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGtt                                                >  2:774749/1‑150 (MQ=255)
                                                                                                        aaGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTc                                             <  2:2277306/150‑1 (MQ=255)
                                                                                                             aaCTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGgat                                        <  1:2424526/150‑1 (MQ=255)
                                                                                                             aaCTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGgat                                        >  1:754734/1‑150 (MQ=255)
                                                                                                             aaCTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGgat                                        >  1:3264227/1‑150 (MQ=255)
                                                                                                                  aGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGt                                   <  2:4906980/150‑1 (MQ=255)
                                                                                                                   gTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGtt                                  <  1:2991950/150‑1 (MQ=255)
                                                                                                                    tGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGttt                                 >  1:3927398/1‑150 (MQ=255)
                                                                                                                    tGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGttt                                 >  1:2768488/1‑150 (MQ=255)
                                                                                                                     gAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTc                                <  2:4564978/150‑1 (MQ=255)
                                                                                                                         ttGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACtt                            <  2:1493689/150‑1 (MQ=255)
                                                                                                                            gTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgt                         >  1:3259125/1‑150 (MQ=255)
                                                                                                                             tGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgtt                        <  1:4386218/150‑1 (MQ=255)
                                                                                                                               aaaaTCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgtttg                      <  1:3515877/150‑1 (MQ=255)
                                                                                                                               aaaaTCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgtttg                      <  2:3286491/150‑1 (MQ=255)
                                                                                                                               aaaaTCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgtttg                      >  1:848873/1‑150 (MQ=255)
                                                                                                                                 aaTCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgtttgtt                    >  1:1802384/1‑150 (MQ=255)
                                                                                                                                 aaTCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgtttgtt                    <  1:4920195/150‑1 (MQ=255)
                                                                                                                                  aTCAGTGTTTACCCGCCATAAGGCTGAGTATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAgtttgttt                   >  2:5134076/1‑150 (MQ=255)
                                                                                                                                    cAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGc                 >  1:975387/1‑150 (MQ=255)
                                                                                                                                     aGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGcc                >  1:1291865/1‑150 (MQ=255)
                                                                                                                                     aGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGcc                <  2:1403771/150‑1 (MQ=255)
                                                                                                                                      gtgtTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCg               <  2:818863/150‑1 (MQ=255)
                                                                                                                                      gtgtTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCg               >  1:161350/1‑150 (MQ=255)
                                                                                                                                         tttACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTAt            <  1:111259/150‑1 (MQ=255)
                                                                                                                                          ttACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATc           >  1:1570772/1‑150 (MQ=255)
                                                                                                                                              ccGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGccc       <  2:4780178/150‑1 (MQ=255)
                                                                                                                                               cGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCg      >  2:4174646/1‑150 (MQ=255)
                                                                                                                                                 ccATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGc    >  2:2383754/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGA  >  CP009273/4529853‑4530149

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: