Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,937 C→T 12.3% intergenic (‑619/‑837) 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,529,9370CT12.3% 874.5 / 94.6 276intergenic (‑619/‑837)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 (133/109);  new base T (19/15);  total (152/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 278 total aligned reads displayed.

AAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGG  >  CP009273/4529788‑4530085
                                                                                                                                                     |                                                                                                                                                    
aaaaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTGTAAATAATTTATTAAGAGTc                                                                                                                                                      <  2:1716660/150‑1 (MQ=255)
aaaaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTc                                                                                                                                                      >  2:4590988/1‑150 (MQ=255)
 aaaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCt                                                                                                                                                     <  2:1415855/150‑1 (MQ=255)
   aCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTaa                                                                                                                                                   <  2:4084714/150‑1 (MQ=255)
    ccAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTaaa                                                                                                                                                  >  2:3515803/1‑150 (MQ=255)
     cAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAAc                                                                                                                                                 <  2:2827367/150‑1 (MQ=255)
       gATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACaa                                                                                                                                               <  1:1297279/150‑1 (MQ=255)
       gATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACaa                                                                                                                                               <  1:3320773/150‑1 (MQ=255)
       gATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACaa                                                                                                                                               <  1:4223072/150‑1 (MQ=255)
        aTTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAg                                                                                                                                              >  1:707436/1‑150 (MQ=255)
          ttGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAggg                                                                                                                                            <  1:1204232/150‑1 (MQ=255)
           tGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAgggg                                                                                                                                           >  2:1063432/1‑150 (MQ=255)
            gCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGa                                                                                                                                          >  2:281275/1‑150 (MQ=255)
             cAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAg                                                                                                                                         >  2:3196349/1‑150 (MQ=255)
              aactCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTAATTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGc                                                                                                                                        >  1:182325/4‑150 (MQ=255)
              aaTTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGc                                                                                                                                        >  1:2679338/1‑150 (MQ=255)
              aaTTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGc                                                                                                                                        >  2:2328488/1‑150 (MQ=255)
               aTTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCt                                                                                                                                       >  2:1304563/1‑150 (MQ=255)
                ttCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCtt                                                                                                                                      >  2:2424526/1‑150 (MQ=255)
                  cGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTg                                                                                                                                    >  1:3843824/1‑150 (MQ=255)
                       cacaAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAATCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGTGGAGCTTTGCAAGc                                                                                                                               >  1:2234965/1‑150 (MQ=255)
                             atatGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTc                                                                                                                         >  2:3019891/1‑150 (MQ=255)
                               atGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAg                                                                                                                       >  2:3563265/1‑150 (MQ=255)
                                tGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGt                                                                                                                      <  1:5035831/150‑1 (MQ=255)
                                 gTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTg                                                                                                                     >  1:3004489/1‑150 (MQ=255)
                                  tCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGa                                                                                                                    <  1:3094855/150‑1 (MQ=255)
                                  tCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGa                                                                                                                    <  1:4917003/150‑1 (MQ=255)
                                   cGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAg                                                                                                                   <  2:4497935/150‑1 (MQ=255)
                                    gATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATATGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGc                                                                                                                  >  2:4429353/1‑150 (MQ=255)
                                        tttttCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTATCTATTTGTTTTATAAATAATTTATTAAGAGTCTAATCAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTgg                                                                                                              <  2:242861/150‑1 (MQ=255)
                                         ttttCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGt                                                                                                             >  1:2758962/1‑150 (MQ=255)
                                            tCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGaa                                                                                                          >  2:4038106/1‑150 (MQ=255)
                                            tCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGaa                                                                                                          >  2:1480413/1‑150 (MQ=255)
                                             cTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGaaa                                                                                                         <  1:3582772/150‑1 (MQ=255)
                                               aaGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAAt                                                                                                       <  1:3670525/150‑1 (MQ=255)
                                                agagGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATc                                                                                                      >  2:2384777/1‑150 (MQ=255)
                                                agagGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCATAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATc                                                                                                      <  2:748556/150‑1 (MQ=255)
                                                 gagGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCa                                                                                                     <  1:1117909/150‑1 (MQ=255)
                                                     aaGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAgtgt                                                                                                 >  1:2013800/1‑150 (MQ=255)
                                                        aTGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTa                                                                                              >  1:3055031/1‑150 (MQ=255)
                                                         tGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGAACTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTAc                                                                                             <  2:645692/150‑1 (MQ=255)
                                                          gCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTAcc                                                                                            >  2:3280294/1‑150 (MQ=255)
                                                                 gAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATc                                                                                     <  1:1747428/150‑1 (MQ=255)
                                                                  aaGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCAtaa                                                                                    <  2:1515555/150‑3 (MQ=255)
                                                                   aGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAg                                                                                   <  2:4345071/150‑1 (MQ=255)
                                                                       aGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTg                                                                               >  1:54374/1‑150 (MQ=255)
                                                                        gACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGa                                                                              >  1:1403236/1‑150 (MQ=255)
                                                                         aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAg                                                                             >  1:1570938/1‑150 (MQ=255)
                                                                         aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAg                                                                             <  2:2413268/150‑1 (MQ=255)
                                                                          cGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAgt                                                                            >  1:4543260/1‑149 (MQ=255)
                                                                           gAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCa                                                                           >  2:111259/1‑150 (MQ=255)
                                                                           gAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCa                                                                           >  2:3399352/1‑150 (MQ=255)
                                                                           gAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCa                                                                           >  2:4440839/1‑150 (MQ=255)
                                                                             acacTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATa                                                                         >  1:4492676/1‑150 (MQ=255)
                                                                               acTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAt                                                                       <  2:1197898/150‑1 (MQ=255)
                                                                               acTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAt                                                                       >  2:2522792/1‑150 (MQ=255)
                                                                                cTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAtt                                                                      <  1:2108128/150‑1 (MQ=255)
                                                                                cTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAtt                                                                      <  1:3260257/150‑1 (MQ=255)
                                                                                   gCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATtctc                                                                   <  1:839466/150‑1 (MQ=255)
                                                                                   gCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATtctc                                                                   <  1:3416073/150‑1 (MQ=255)
                                                                                    cGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtca                                                                  <  2:2896628/150‑1 (MQ=255)
                                                                                     gttggtATTCAAAAACTAAAGATCTTATTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTTGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcat                                                                 <  2:2212411/148‑1 (MQ=255)
                                                                                     ggtggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCtcat                                                                 >  1:1018904/1‑150 (MQ=255)
                                                                                      gtggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcatc                                                                <  2:2106203/150‑1 (MQ=255)
                                                                                       tggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcatca                                                               <  1:482189/150‑1 (MQ=255)
                                                                                        ggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTcatcat                                                              >  1:2408887/1‑150 (MQ=255)
                                                                                        ggtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTcatcat                                                              >  1:2501297/1‑150 (MQ=255)
                                                                                         gtCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATg                                                             >  1:4770984/1‑150 (MQ=255)
                                                                                          tCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGa                                                            >  2:1020638/1‑150 (MQ=255)
                                                                                             tCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAat                                                         <  2:2673265/150‑1 (MQ=255)
                                                                                              cAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAata                                                        >  1:580182/1‑150 (MQ=255)
                                                                                                   aCTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGttt                                                   <  1:2610918/150‑1 (MQ=255)
                                                                                                   aCTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTAGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATACTCATCATGAAATATGttt                                                   >  1:1513293/1‑150 (MQ=255)
                                                                                                    cTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTc                                                  <  2:1733831/150‑1 (MQ=255)
                                                                                                      aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCt                                                >  1:5129740/1‑150 (MQ=255)
                                                                                                      aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCt                                                >  2:204641/1‑150 (MQ=255)
                                                                                                      aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCt                                                <  2:1373825/150‑1 (MQ=255)
                                                                                                      aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCt                                                <  2:3937761/150‑1 (MQ=255)
                                                                                                       aaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGTATAATTCTCATCATGAAATATGTTTCCTg                                               <  1:4129616/150‑1 (MQ=255)
                                                                                                           tCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGttt                                           <  1:629012/150‑1 (MQ=255)
                                                                                                            cTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTtg                                          >  1:2953761/1‑150 (MQ=255)
                                                                                                              tAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGTATAATTCTCATCATGAAATATGTTTCCTGGTTtgtg                                        >  2:1766556/1‑150 (MQ=255)
                                                                                                               aGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTgg                                       >  2:5041496/1‑150 (MQ=255)
                                                                                                                    aaCTATTTGTTTTATAAATAATTTATTAAGAGTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGt                                  >  2:489000/1‑150 (MQ=255)
                                                                                                                    aaCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGt                                  >  1:1100766/1‑150 (MQ=255)
                                                                                                                     aCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGATCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTa                                 >  1:499548/1‑150 (MQ=255)
                                                                                                                     aCTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTa                                 >  2:4733981/1‑150 (MQ=255)
                                                                                                                        aTTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACt                              <  2:3595021/150‑1 (MQ=255)
                                                                                                                               ttATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACt                       <  2:3843824/150‑1 (MQ=255)
                                                                                                                               ttATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGGCACt                       >  2:1700064/1‑150 (MQ=255)
                                                                                                                                  taaataaTTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCt                    <  2:534087/150‑1 (MQ=255)
                                                                                                                                    aataatTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAACTATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGa                  >  1:333152/1‑150 (MQ=255)
                                                                                                                                        atTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTc              >  1:2872476/1‑150 (MQ=255)
                                                                                                                                        atTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTc              <  1:4414521/150‑1 (MQ=255)
                                                                                                                                          ttattaAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGa            >  2:4300123/1‑150 (MQ=255)
                                                                                                                                            attaAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATc          >  1:147039/1‑150 (MQ=255)
                                                                                                                                            attaAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATc          <  1:433029/150‑1 (MQ=255)
                                                                                                                                             ttaAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCt         >  1:4016651/1‑150 (MQ=255)
                                                                                                                                                  agTTTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATAAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGgaga    >  2:4199445/1‑150 (MQ=255)
                                                                                                                                                  agTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGgaga    >  1:4956757/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGG  >  CP009273/4529788‑4530085

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: