Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,263 A→T 8.2% L19* (TTA→TAA)  nanC ← N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,529,2630AT 8.2% 780.6 / 47.9 232L19* (TTA→TAA) nanCN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base A (88/125);  new base T (7/12);  total (95/137)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 237 total aligned reads displayed.

TATTATCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATA  >  CP009273/4529114‑4529412
                                                                                                                                                     |                                                                                                                                                     
tattatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTa                                                                                                                                                       <  2:3630624/150‑1 (MQ=255)
tattatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTa                                                                                                                                                       <  2:3024729/150‑1 (MQ=255)
   tatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATg                                                                                                                                                    <  2:2016554/150‑1 (MQ=255)
   tatCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATg                                                                                                                                                    >  1:609711/1‑150 (MQ=255)
    atCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATgg                                                                                                                                                   >  2:3285232/1‑150 (MQ=255)
     tCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATggg                                                                                                                                                  <  2:3343967/150‑1 (MQ=255)
     tCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATggg                                                                                                                                                  <  2:2283469/150‑1 (MQ=255)
      cATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATgggg                                                                                                                                                 >  2:2065711/1‑150 (MQ=255)
      cATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATgggg                                                                                                                                                 >  2:686294/1‑150 (MQ=255)
      cATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATgggg                                                                                                                                                 <  1:2910765/150‑1 (MQ=255)
       aTGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGa                                                                                                                                                <  1:790531/150‑1 (MQ=255)
        tGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGAc                                                                                                                                               <  1:1019252/150‑1 (MQ=255)
        tGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGAc                                                                                                                                               >  2:2230769/1‑150 (MQ=255)
         gAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACg                                                                                                                                              >  1:1770797/1‑150 (MQ=255)
           aTGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGAGACGaa                                                                                                                                            <  2:3514009/150‑1 (MQ=255)
             ggTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaaaa                                                                                                                                          >  1:4129917/1‑150 (MQ=255)
             ggTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGaaaa                                                                                                                                          <  1:4461656/150‑1 (MQ=255)
              gTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATGGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAg                                                                                                                                         <  1:3902263/150‑1 (MQ=255)
              gTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAg                                                                                                                                         <  1:3749841/150‑1 (MQ=255)
               tATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGc                                                                                                                                        >  1:68857/1‑150 (MQ=255)
               tATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGc                                                                                                                                        <  1:4404130/150‑1 (MQ=255)
                 ttCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGcac                                                                                                                                      <  1:1922617/150‑1 (MQ=255)
                   ccAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAg                                                                                                                                    >  2:2156968/1‑150 (MQ=255)
                   ccAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAg                                                                                                                                    <  1:3146680/150‑1 (MQ=255)
                    cAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGt                                                                                                                                   <  2:1092776/150‑1 (MQ=255)
                     aGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGta                                                                                                                                  <  1:1409672/150‑1 (MQ=255)
                      ggTATTACTTTCCATGCTTGCCCACCATTCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaa                                                                                                                                 >  1:2568688/1‑150 (MQ=255)
                      ggTATTACTTTCCATGCTTGCCCACCATTCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaa                                                                                                                                 >  1:300124/1‑150 (MQ=255)
                       gTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaat                                                                                                                                <  1:2030859/150‑1 (MQ=255)
                        tATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGtaata                                                                                                                               >  1:4369271/1‑150 (MQ=255)
                         aTTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGtaataa                                                                                                                              >  1:370063/1‑150 (MQ=255)
                          ttACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTaataat                                                                                                                             >  1:1574584/1‑150 (MQ=255)
                           tACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAataata                                                                                                                            <  1:48437/150‑1 (MQ=255)
                            aCTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATAc                                                                                                                           >  1:4319413/1‑150 (MQ=255)
                                tCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCa                                                                                                                       >  2:1133341/1‑150 (MQ=255)
                                  cATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGa                                                                                                                     <  2:2152262/150‑1 (MQ=255)
                                  cATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGa                                                                                                                     <  1:1501847/150‑1 (MQ=255)
                                   aTGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGaa                                                                                                                    >  1:1781153/1‑150 (MQ=255)
                                      cTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGt                                                                                                                 <  2:2736457/150‑1 (MQ=255)
                                      cTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGt                                                                                                                 >  1:333329/1‑150 (MQ=255)
                                          cccACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAttt                                                                                                             <  2:1213744/150‑1 (MQ=255)
                                           ccaccaTCCATTTTGCCGTCCCTCACTGACATTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAtttt                                                                                                            <  2:1368716/150‑1 (MQ=255)
                                           ccaccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTAtttt                                                                                                            <  2:3652371/150‑1 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                           <  1:4198379/150‑1 (MQ=255)
                                            caccaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTa                                                                                                           <  1:1323508/150‑1 (MQ=255)
                                             accaTCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAg                                                                                                          <  1:2789335/150‑1 (MQ=255)
                                               catccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGcc                                                                                                        <  1:322528/150‑1 (MQ=255)
                                               catccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGcc                                                                                                        <  2:3233469/150‑1 (MQ=255)
                                                 tccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtt                                                                                                      <  1:70507/150‑1 (MQ=255)
                                                  ccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCttt                                                                                                     >  2:3884256/1‑150 (MQ=255)
                                                  ccatTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCttt                                                                                                     >  1:1946558/1‑150 (MQ=255)
                                                   catTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtttt                                                                                                    <  1:3052972/150‑1 (MQ=255)
                                                   catTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCtttt                                                                                                    <  1:1518045/150‑1 (MQ=255)
                                                        tGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATa                                                                                               >  1:4069825/1‑150 (MQ=255)
                                                         gCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATaa                                                                                              <  1:3713406/150‑1 (MQ=255)
                                                             tCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAAttt                                                                                          >  1:2215604/1‑150 (MQ=255)
                                                                 tCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACt                                                                                      >  1:1007/1‑150 (MQ=255)
                                                                   aCTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACGACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCa                                                                                    <  2:2459584/150‑1 (MQ=255)
                                                                    cTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAt                                                                                   <  2:4062365/150‑1 (MQ=255)
                                                                    cTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAt                                                                                   >  2:2673659/1‑150 (MQ=255)
                                                                      gACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAttt                                                                                 <  1:1433773/150‑1 (MQ=255)
                                                                      gACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAttt                                                                                 <  1:2729855/150‑1 (MQ=255)
                                                                      gACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGAAATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCAttt                                                                                 >  1:2672501/1‑150 (MQ=255)
                                                                        cTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGt                                                                               >  1:636579/1‑150 (MQ=255)
                                                                         tttGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTa                                                                              >  1:572547/1‑150 (MQ=255)
                                                                          ttGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAg                                                                             <  2:2609977/150‑1 (MQ=255)
                                                                             agtcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTTGGAt                                                                          >  2:3676337/1‑150 (MQ=255)
                                                                             agtcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGAt                                                                          <  1:1115068/150‑1 (MQ=255)
                                                                             agtcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGAt                                                                          <  1:3986289/150‑1 (MQ=255)
                                                                             agtcgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGAt                                                                          >  2:1712297/1‑150 (MQ=255)
                                                                                cgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACa                                                                       >  1:206818/1‑150 (MQ=255)
                                                                                cgagtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACa                                                                       >  1:3919652/1‑150 (MQ=255)
                                                                                  agtcTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGa                                                                     <  1:4471618/150‑1 (MQ=255)
                                                                                     ctcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAg                                                                  >  2:1599749/1‑150 (MQ=255)
                                                                                     ctcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAg                                                                  <  2:2750808/150‑1 (MQ=255)
                                                                                      tcATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGc                                                                 <  2:3025256/150‑1 (MQ=255)
                                                                                                   ttCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCg                                                    >  1:4525286/1‑150 (MQ=255)
                                                                                                      cACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCAt                                                 <  1:3456542/150‑1 (MQ=255)
                                                                                                          aCGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAAc                                             <  2:2972989/150‑1 (MQ=255)
                                                                                                           cGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACa                                            <  2:513471/150‑1 (MQ=255)
                                                                                                            gATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACaa                                           >  2:2395428/1‑150 (MQ=255)
                                                                                                                    ccacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCAtt                                   <  2:1176015/150‑1 (MQ=255)
                                                                                                                     cacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAACAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATGTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACACTTAGCAttt                                  <  2:2505923/150‑1 (MQ=255)
                                                                                                                           cgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAGTACAAAGCCCGCATAAACAATTAGCATTTAtgttg                            >  1:2172370/1‑150 (MQ=255)
                                                                                                                             tCCAGTGTCGCAGCCTGAGAAATTTATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTtgtg                          <  1:184921/150‑1 (MQ=255)
                                                                                                                              ccAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTgtgt                         <  1:842717/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAat                      <  2:3860036/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAat                      <  2:4325573/150‑1 (MQ=255)
                                                                                                                                   gtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAatat                    <  1:2106597/150‑1 (MQ=255)
                                                                                                                                   gtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAatat                    <  1:1683547/150‑1 (MQ=255)
                                                                                                                                       cAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAttttt                <  2:412146/150‑1 (MQ=255)
                                                                                                                                        aGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAtttttt               >  2:4301444/1‑150 (MQ=255)
                                                                                                                                        aGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAtttttt               >  1:854098/1‑150 (MQ=255)
                                                                                                                                        aGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAtttttt               >  2:3256330/1‑150 (MQ=255)
                                                                                                                                         gCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAttttttt              <  2:4138534/150‑1 (MQ=255)
                                                                                                                                          ccTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTg             <  2:1423118/150‑1 (MQ=255)
                                                                                                                                            tGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGcc           >  2:1046982/1‑150 (MQ=255)
                                                                                                                                                aaaTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGc       >  2:738205/1‑150 (MQ=255)
                                                                                                                                                 aaTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCt      >  1:1199134/1‑150 (MQ=255)
                                                                                                                                                 aaTTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCt      <  1:2286242/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTATCATGAATGGTATTCCAGGTATTACTTTCCATGCTTGCCCACCATCCATTTTGCCATCCCTCACTGACTTTGAGTCGAGTCTCATAGGCGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATA  >  CP009273/4529114‑4529412

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: