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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,184,6540AT14.5% 689.4 / 94.0 239N222Y (AAT→TAT) yqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base A (105/99);  new base T (15/20);  total (120/119)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

ATAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGG  >  CP009273/3184505‑3184803
                                                                                                                                                     |                                                                                                                                                     
aTAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTa                                                                                                                                                       >  2:49261/1‑150 (MQ=255)
aTAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTa                                                                                                                                                       <  2:3061104/150‑1 (MQ=255)
 tAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTaa                                                                                                                                                      >  1:1473444/1‑150 (MQ=255)
 tAAAGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTaa                                                                                                                                                      >  1:4609646/1‑150 (MQ=255)
   aaGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATg                                                                                                                                                    <  2:3869133/150‑1 (MQ=255)
    aGTAAATAAAAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGt                                                                                                                                                   >  1:2901077/1‑150 (MQ=255)
      taaataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTg                                                                                                                                                 <  2:1923142/150‑1 (MQ=255)
       aaataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGa                                                                                                                                                >  1:3769552/1‑150 (MQ=255)
        aataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGAt                                                                                                                                               <  1:1116104/150‑1 (MQ=255)
        aataaaAAACTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGAt                                                                                                                                               >  1:4289087/1‑150 (MQ=255)
            aaaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTgg                                                                                                                                           <  1:2765531/150‑1 (MQ=255)
             aaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGa                                                                                                                                          >  2:4642122/1‑150 (MQ=255)
             aaaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGa                                                                                                                                          <  1:3404306/150‑1 (MQ=255)
              aaaCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGaa                                                                                                                                         >  1:1764702/1‑150 (MQ=255)
                aCTTGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAAcc                                                                                                                                       >  2:3203493/1‑150 (MQ=255)
                  ttGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAg                                                                                                                                     <  2:4549448/150‑1 (MQ=255)
                   tGTATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGc                                                                                                                                    <  1:3066591/150‑1 (MQ=255)
                     tATTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTg                                                                                                                                  >  2:1035345/1‑150 (MQ=255)
                      aTTACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGa                                                                                                                                 >  1:2471491/1‑150 (MQ=255)
                       ttACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGAt                                                                                                                                >  1:2968711/1‑150 (MQ=255)
                       ttACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGAt                                                                                                                                >  1:2968981/1‑150 (MQ=255)
                       ttACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGAt                                                                                                                                >  1:2138727/1‑150 (MQ=255)
                        tACCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATg                                                                                                                               >  2:2336366/1‑150 (MQ=255)
                         aCCTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATgg                                                                                                                              <  1:1189405/150‑1 (MQ=255)
                          ccTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGc                                                                                                                             <  2:914101/150‑1 (MQ=255)
                           cTGACGAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCg                                                                                                                            <  1:769379/150‑1 (MQ=255)
                                gAATATATAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTc                                                                                                                       >  1:4463463/1‑150 (MQ=255)
                                  atatatAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGa                                                                                                                     <  2:3253824/150‑1 (MQ=255)
                                  atatatAGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGa                                                                                                                     <  2:3251746/150‑1 (MQ=255)
                                   tatataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGAc                                                                                                                    <  2:1186676/150‑1 (MQ=255)
                                      ataGCACTTGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACttt                                                                                                                 <  1:3779219/150‑1 (MQ=255)
                                              tGAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaat                                                                                                         <  1:2578513/150‑1 (MQ=255)
                                               gAATTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaata                                                                                                        <  1:1804126/150‑1 (MQ=255)
                                                aaTTTAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataa                                                                                                       >  2:847369/1‑150 (MQ=255)
                                                  tttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataaat                                                                                                     <  2:865575/150‑1 (MQ=255)
                                                  tttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataaat                                                                                                     <  2:619630/150‑1 (MQ=255)
                                                  tttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAaaataaat                                                                                                     >  1:1481826/1‑150 (MQ=255)
                                                   ttAAAGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAAtt                                                                                                    <  1:3531273/150‑1 (MQ=255)
                                                     aaaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTcc                                                                                                  <  1:1554218/150‑1 (MQ=255)
                                                      aaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTccc                                                                                                 >  2:3646930/1‑150 (MQ=255)
                                                      aaGGTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTccc                                                                                                 >  1:454947/1‑150 (MQ=255)
                                                        ggTATGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGa                                                                                               <  2:361326/150‑1 (MQ=255)
                                                           aTGAGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTACCGAACa                                                                                            <  1:3261220/150‑1 (MQ=255)
                                                              aGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAAtt                                                                                         >  2:2679270/1‑150 (MQ=255)
                                                              aGCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAAtt                                                                                         <  1:555145/150‑1 (MQ=255)
                                                               gCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTa                                                                                        <  2:1211956/150‑1 (MQ=255)
                                                               gCGGAGCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTa                                                                                        <  1:4024443/150‑1 (MQ=255)
                                                                    gCAGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAAt                                                                                   >  1:2528364/1‑150 (MQ=255)
                                                                      aGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATAc                                                                                 <  2:1827490/150‑1 (MQ=255)
                                                                      aGGTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATAc                                                                                 >  2:348722/1‑150 (MQ=255)
                                                                        gTGCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCa                                                                               <  1:939059/150‑1 (MQ=255)
                                                                          gCTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGc                                                                             <  2:21560/150‑1 (MQ=255)
                                                                           cTATAGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCg                                                                            <  2:4442187/150‑1 (MQ=255)
                                                                             ataGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtg                                                                          >  2:2650057/1‑150 (MQ=255)
                                                                             ataGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtg                                                                          <  2:1402216/150‑1 (MQ=255)
                                                                               aGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgag                                                                        <  1:1706321/150‑1 (MQ=255)
                                                                               aGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgag                                                                        >  1:52147/1‑150 (MQ=255)
                                                                               aGAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGAATGCTTTGTTAATGTTGATCTGGAACCAGATGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCttgag                                                                        >  1:802146/1‑150 (MQ=255)
                                                                                gAAGTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATATCTTTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCgtgagt                                                                       <  1:434641/150‑1 (MQ=255)
                                                                                   gTAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGtgagtgag                                                                    >  2:2550145/1‑150 (MQ=255)
                                                                                    tAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTgaga                                                                   >  2:3791476/1‑150 (MQ=255)
                                                                                    tAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTgaga                                                                   <  1:2702325/150‑1 (MQ=255)
                                                                                    tAGATAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTgaga                                                                   <  2:3154001/150‑1 (MQ=255)
                                                                                        tAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGtt                                                               <  1:3252951/150‑1 (MQ=255)
                                                                                        tAAAAATCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGtt                                                               <  1:736956/150‑1 (MQ=255)
                                                                                          aaaaTCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTa                                                             >  2:4115473/1‑150 (MQ=255)
                                                                                           aaaTCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAg                                                            >  2:402929/1‑150 (MQ=255)
                                                                                             aTCTGACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCg                                                          >  1:3465227/1‑150 (MQ=255)
                                                                                                 gACCTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGt                                                      <  1:4449855/150‑1 (MQ=255)
                                                                                                   ccTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCa                                                    <  1:223679/150‑1 (MQ=255)
                                                                                                    cTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCAt                                                   <  2:1588744/150‑1 (MQ=255)
                                                                                                    cTTCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCAt                                                   >  2:3820664/1‑150 (MQ=255)
                                                                                                     ttCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATg                                                  >  2:2522375/1‑150 (MQ=255)
                                                                                                     ttCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATg                                                  <  2:2307930/150‑1 (MQ=255)
                                                                                                      tCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGa                                                 <  2:4185987/150‑1 (MQ=255)
                                                                                                      tCAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGa                                                 <  2:4170978/150‑1 (MQ=255)
                                                                                                       cAGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGAc                                                >  2:2222884/1‑150 (MQ=255)
                                                                                                        aGGATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGAcc                                               >  2:4247345/1‑150 (MQ=255)
                                                                                                         ggATTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCa                                              <  1:2080812/150‑1 (MQ=255)
                                                                                                           aTTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCaaa                                            <  2:4599363/150‑1 (MQ=255)
                                                                                                           aTTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCaaa                                            >  1:1093264/1‑150 (MQ=255)
                                                                                                              cgcgGATTAAATAACATTCATGTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGAt                                         >  2:3034387/1‑150 (MQ=255)
                                                                                                               gcgGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATc                                        >  1:651845/1‑150 (MQ=255)
                                                                                                                    ttAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGt                                   <  2:4419205/150‑1 (MQ=255)
                                                                                                                       aaTAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGTGTGCg                                <  2:2442906/150‑1 (MQ=255)
                                                                                                                          aaCATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGcc                             <  1:3568296/150‑1 (MQ=255)
                                                                                                                           aCATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCt                            >  1:2023266/1‑150 (MQ=255)
                                                                                                                            cATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTg                           <  2:3307230/150‑1 (MQ=255)
                                                                                                                             aTTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGt                          <  1:3075963/150‑1 (MQ=255)
                                                                                                                                 aTGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTg                      >  1:3971171/1‑150 (MQ=255)
                                                                                                                                  tGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCATAGTGAGTGAGACGTTTAGCGTAGTCATGACAAAAGATCAGGGTGCGTCCTGTACTGa                     >  2:4568288/1‑150 (MQ=255)
                                                                                                                                   gTCCTTGACTGCTTTGTTTATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAg                    >  2:3506611/1‑150 (MQ=255)
                                                                                                                                   gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAg                    <  2:2925954/150‑1 (MQ=255)
                                                                                                                                     ccTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCa                  <  1:1344229/150‑1 (MQ=255)
                                                                                                                                      cTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAg                 >  2:2820008/1‑150 (MQ=255)
                                                                                                                                      cTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAg                 >  1:3980517/1‑150 (MQ=255)
                                                                                                                                        tGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGtt               >  2:1894968/1‑150 (MQ=255)
                                                                                                                                              cTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAatat         <  2:4149144/150‑1 (MQ=255)
                                                                                                                                               tttGTTAATGTTGATCTGGAACCAGTTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATAtt        <  1:3194487/150‑1 (MQ=255)
                                                                                                                                                     aaTGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAggg  >  1:2105322/1‑150 (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: ‑