Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,184,654 A→T 10.2% N222Y (AAT→TAT)  yqiI → fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,184,6540AT10.2% 785.4 / 63.2 244N222Y (AAT→TAT) yqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base A (97/122);  new base T (16/9);  total (113/131)
Fisher's exact test for biased strand distribution p-value = 8.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 253 total aligned reads displayed.

ATAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGG  >  CP009273/3184505‑3184803
                                                                                                                                                     |                                                                                                                                                     
aTAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTa                                                                                                                                                       <  1:652552/150‑1 (MQ=255)
 tAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTaa                                                                                                                                                      >  1:4133517/1‑150 (MQ=255)
 tAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTaa                                                                                                                                                      <  2:1267753/150‑1 (MQ=255)
  aaaGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAAt                                                                                                                                                     <  1:2970411/150‑1 (MQ=255)
   aaGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATg                                                                                                                                                    >  2:36424/1‑150 (MQ=255)
     gTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGtt                                                                                                                                                  <  1:1750727/150‑1 (MQ=255)
      taaataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTg                                                                                                                                                 >  1:3303072/1‑150 (MQ=255)
      taaataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTg                                                                                                                                                 >  1:2744975/1‑150 (MQ=255)
      taaataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTg                                                                                                                                                 >  2:3707514/1‑150 (MQ=255)
        aataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGAt                                                                                                                                               >  2:2659220/1‑150 (MQ=255)
           aaaaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTg                                                                                                                                            <  1:1984119/150‑1 (MQ=255)
           aaaaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTGTATGTTGATCTg                                                                                                                                            >  1:809568/1‑150 (MQ=255)
              aaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGaa                                                                                                                                         >  2:2935180/1‑150 (MQ=255)
               aaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAAc                                                                                                                                        <  1:430640/150‑1 (MQ=255)
                    gTATTACCTGACGACTATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCt                                                                                                                                   <  1:1763974/150‑1 (MQ=255)
                    gTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTCTGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTATGTTAATGTTGATCTGGAACCAGCt                                                                                                                                   >  2:539040/1‑150 (MQ=255)
                    gTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCt                                                                                                                                   <  1:2292001/150‑1 (MQ=255)
                     tATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTg                                                                                                                                  >  2:141732/1‑150 (MQ=255)
                     tATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACAGTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTg                                                                                                                                  <  1:2152398/150‑1 (MQ=255)
                      aTTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGa                                                                                                                                 >  1:3894196/1‑150 (MQ=255)
                      aTTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGa                                                                                                                                 >  1:251035/1‑150 (MQ=255)
                      aTTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGa                                                                                                                                 >  2:2760397/1‑150 (MQ=255)
                       ttACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGAt                                                                                                                                <  2:2674265/150‑1 (MQ=255)
                       ttACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGAt                                                                                                                                >  2:4120985/1‑150 (MQ=255)
                        tACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATg                                                                                                                               >  1:567356/1‑150 (MQ=255)
                        tACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATg                                                                                                                               >  1:88688/1‑150 (MQ=255)
                        tACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATg                                                                                                                               >  1:1451514/1‑150 (MQ=255)
                         aCCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATgg                                                                                                                              <  1:1827436/150‑1 (MQ=255)
                         aCCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATgg                                                                                                                              <  2:2217549/150‑1 (MQ=255)
                         aCCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATgg                                                                                                                              >  1:3349519/1‑150 (MQ=255)
                          ccTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGc                                                                                                                             >  1:2459478/1‑150 (MQ=255)
                          ccTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGc                                                                                                                             <  1:3116697/150‑1 (MQ=255)
                            tGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGt                                                                                                                           >  1:698776/1‑150 (MQ=255)
                             gACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGtt                                                                                                                          <  2:200602/150‑1 (MQ=255)
                               cGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGt                                                                                                                        >  1:1142145/1‑150 (MQ=255)
                                gAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTc                                                                                                                       <  1:2787336/150‑1 (MQ=255)
                                gAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTc                                                                                                                       <  2:2990623/150‑1 (MQ=255)
                                   tatataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGAc                                                                                                                    >  1:2099295/1‑150 (MQ=255)
                                   tatataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGAc                                                                                                                    <  1:3541040/150‑1 (MQ=255)
                                      ataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTCTAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACttt                                                                                                                 >  2:2502422/1‑150 (MQ=255)
                                       taGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTg                                                                                                                >  2:1906949/1‑150 (MQ=255)
                                        aGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTgg                                                                                                               >  2:1700763/1‑150 (MQ=255)
                                           aCTTGAATTTTAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATATAAATCTGACCTTCAGGTTTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTTGAACCAGCTGTTGGCGTTGTCGACTTTGGTaa                                                                                                            <  1:2180609/150‑1 (MQ=255)
                                           aCTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTaa                                                                                                            <  1:1881396/150‑1 (MQ=255)
                                           aCTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTaa                                                                                                            >  1:267260/1‑150 (MQ=255)
                                               gAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaata                                                                                                        <  1:317863/150‑1 (MQ=255)
                                                 aTTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataaa                                                                                                      <  1:1262308/150‑1 (MQ=255)
                                                   ttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAAtt                                                                                                    <  1:2563833/150‑1 (MQ=255)
                                                     aaaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAACATCTGACCGTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTcc                                                                                                  <  2:1881126/150‑1 (MQ=255)
                                                     aaaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTcc                                                                                                  >  1:1254103/1‑150 (MQ=255)
                                                       aGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCg                                                                                                >  1:1029375/1‑150 (MQ=255)
                                                         gTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGaa                                                                                              >  1:3638710/1‑150 (MQ=255)
                                                           aTGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACa                                                                                            <  2:3996802/150‑1 (MQ=255)
                                                           aTGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACa                                                                                            <  2:3189016/150‑1 (MQ=255)
                                                              aGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAAtt                                                                                         <  1:2363564/150‑1 (MQ=255)
                                                               gCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTa                                                                                        <  2:3033563/150‑1 (MQ=255)
                                                                 ggAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTaaa                                                                                      <  2:96502/150‑1 (MQ=255)
                                                                   aGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTaaaaa                                                                                    <  1:1017861/150‑1 (MQ=255)
                                                                     cAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATa                                                                                  <  1:1606017/150‑1 (MQ=255)
                                                                     cAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATa                                                                                  <  2:3903244/150‑1 (MQ=255)
                                                                     cAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATa                                                                                  >  2:150864/1‑150 (MQ=255)
                                                                       ggTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATAcc                                                                                >  1:243649/1‑150 (MQ=255)
                                                                        gTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCa                                                                               <  2:3181360/150‑1 (MQ=255)
                                                                         tGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAg                                                                              <  2:3818285/150‑1 (MQ=255)
                                                                             ataGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtg                                                                          >  2:2283906/1‑150 (MQ=255)
                                                                                gAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgagt                                                                       >  2:2147068/1‑150 (MQ=255)
                                                                                  aGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgagtga                                                                     <  2:1073534/150‑1 (MQ=255)
                                                                                  aGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgagtga                                                                     >  1:2686513/1‑150 (MQ=255)
                                                                                         aaaaaTCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGGTTTGTTAATGTTGATCTGGAACCAGCTGTTGGCGTTGTCGACTTTGGTAAAATAAAGTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGttt                                                              <  2:1292346/150‑1 (MQ=255)
                                                                                         aaaaaTCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGttt                                                              <  1:3840767/150‑1 (MQ=255)
                                                                                                  aCCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTc                                                     >  1:406280/1‑150 (MQ=255)
                                                                                                   ccTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCa                                                    >  2:2856474/1‑150 (MQ=255)
                                                                                                    cTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCac                                                   <  2:831863/150‑2 (MQ=255)
                                                                                                      tCAGGATTCGCGGATTAAATAACATTCATGTCCGTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGa                                                 >  2:1250814/1‑150 (MQ=255)
                                                                                                        aGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTCGTCATGAcc                                               >  2:59585/1‑150 (MQ=255)
                                                                                                        aGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGAcc                                               >  2:3864707/1‑150 (MQ=255)
                                                                                                             tCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGa                                          <  1:4007140/150‑1 (MQ=255)
                                                                                                              cgcgGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGAt                                         >  2:2927257/1‑150 (MQ=255)
                                                                                                              cgcgGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGAt                                         >  1:2295212/1‑150 (MQ=255)
                                                                                                                  gATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCgg                                     <  1:2328320/150‑1 (MQ=255)
                                                                                                                     tAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTg                                  <  1:3235998/150‑1 (MQ=255)
                                                                                                                      aaaTAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGc                                 <  2:2873827/150‑1 (MQ=255)
                                                                                                                           aCATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCt                            <  2:3848473/150‑1 (MQ=255)
                                                                                                                            cATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTg                           <  1:2749576/150‑1 (MQ=255)
                                                                                                                            cATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTg                           >  2:1035091/1‑150 (MQ=255)
                                                                                                                             aTTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGt                          <  2:3757260/150‑1 (MQ=255)
                                                                                                                             aTTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGt                          <  1:3011858/150‑1 (MQ=255)
                                                                                                                                cATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACt                       <  1:163271/150‑1 (MQ=255)
                                                                                                                                cATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAAGTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACt                       <  1:1094945/150‑1 (MQ=255)
                                                                                                                                   gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAg                    >  2:639367/1‑150 (MQ=255)
                                                                                                                                       ttGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGt                <  1:978042/150‑1 (MQ=255)
                                                                                                                                        tGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGtt               >  2:115769/1‑150 (MQ=255)
                                                                                                                                         gACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGttt              <  1:2036549/150‑1 (MQ=255)
                                                                                                                                          aCTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTa             >  2:3373227/1‑150 (MQ=255)
                                                                                                                                          aCTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTa             <  2:2682661/150‑1 (MQ=255)
                                                                                                                                             gCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAata          >  2:2364245/1‑150 (MQ=255)
                                                                                                                                               tttGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAtt        <  1:1415376/150‑1 (MQ=255)
                                                                                                                                                ttGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAttt       <  2:797143/150‑1 (MQ=255)
                                                                                                                                                  gTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTa     <  2:1536721/150‑1 (MQ=255)
                                                                                                                                                     aaTGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAggg  <  1:2760439/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGG  >  CP009273/3184505‑3184803

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