Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,253,172 G→A 9.7% intergenic (+30/+38) argG → / ← PP_1089 argininosuccinate synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,253,1720GA 9.7% 709.5 / 49.2 246intergenic (+30/+38)argG/PP_1089argininosuccinate synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (123/99);  new base A (15/9);  total (138/108)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 254 total aligned reads displayed.

AACATCGCCACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGCGGTTGG  >  NC_002947/1253023‑1253321
                                                                                                                                                     |                                                                                                                                                      
aaCATCGCCACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCgg                                                                                                                                                        >  1:234852/1‑150 (MQ=255)
     cGCCACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCgg                                                                                                                                                   <  1:2407872/150‑1 (MQ=255)
      gCCACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCggg                                                                                                                                                  <  1:2796116/150‑1 (MQ=255)
         aCCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGACCGAGCCGGGGtt                                                                                                                                              <  2:2392963/150‑1 (MQ=255)
          ccTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGACCGAGCCGGGGttt                                                                                                                                             <  2:4100840/150‑1 (MQ=255)
            tttGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGttttt                                                                                                                                           >  2:44393/1‑150 (MQ=255)
            tttGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGACCGAGCCGGGGttttc                                                                                                                                           <  2:4061279/150‑2 (MQ=255)
            tttGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGACCGAGCCGGGGttttc                                                                                                                                           <  2:2508036/150‑2 (MQ=255)
              tGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGttttttt                                                                                                                                         <  2:683157/150‑1 (MQ=255)
               gaagaaGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTa                                                                                                                                        <  1:2563441/150‑1 (MQ=255)
                    agaTGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGccc                                                                                                                                   <  2:5509297/150‑1 (MQ=255)
                    agaTGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGccc                                                                                                                                   >  1:2923790/1‑150 (MQ=255)
                    agaTGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGccc                                                                                                                                   <  1:919684/150‑1 (MQ=255)
                           ggtgCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCaaa                                                                                                                            <  2:571272/150‑1 (MQ=255)
                              gCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGc                                                                                                                         <  2:4783745/150‑1 (MQ=255)
                              gCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGc                                                                                                                         <  2:3295115/150‑1 (MQ=255)
                                cTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCgg                                                                                                                       <  1:3646291/150‑1 (MQ=255)
                                 tACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                                      >  1:5473345/1‑150 (MQ=255)
                                 tACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCggg                                                                                                                      >  1:2892508/1‑150 (MQ=255)
                                  aCAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGt                                                                                                                     >  2:2093177/1‑150 (MQ=255)
                                    aaCCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCg                                                                                                                   >  1:3314985/1‑150 (MQ=255)
                                    aaCCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCg                                                                                                                   >  2:2319352/1‑150 (MQ=255)
                                       cAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTc                                                                                                                <  1:4297058/150‑1 (MQ=255)
                                        aGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCa                                                                                                               >  2:69497/1‑150 (MQ=255)
                                          gCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAAc                                                                                                             >  1:354279/1‑150 (MQ=255)
                                            aGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGGTCAACtt                                                                                                           >  2:4927575/1‑150 (MQ=255)
                                             gATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTc                                                                                                          <  1:1080149/150‑1 (MQ=255)
                                                gctgctGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCgg                                                                                                        <  2:1687085/150‑1 (MQ=255)
                                                gctgctGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGt                                                                                                       <  2:3122448/150‑1 (MQ=255)
                                                gctgctGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGt                                                                                                       <  1:5404816/150‑1 (MQ=255)
                                                         ttCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCAcat                                                                                              >  2:4297887/1‑150 (MQ=255)
                                                         ttCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCAcat                                                                                              >  1:4522062/1‑150 (MQ=255)
                                                            atcaAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCAcatcat                                                                                           >  2:5378978/1‑150 (MQ=255)
                                                             tcaAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACatcatc                                                                                          <  1:1887852/150‑1 (MQ=255)
                                                             tcaAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAACGGCCGTTCGCTGCTCTGATTGCTGGCGTGAGCGAAAACCCCGGCTCGGTCCGGGGTTATATTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACatcatc                                                                                          <  2:3360128/150‑1 (MQ=255)
                                                              caAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCg                                                                                         <  2:2877115/150‑1 (MQ=255)
                                                               aaGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGc                                                                                        >  2:369854/1‑150 (MQ=255)
                                                               aaGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGc                                                                                        >  1:4106031/1‑150 (MQ=255)
                                                                aGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGcc                                                                                       <  2:4728178/150‑1 (MQ=255)
                                                                 gCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGccc                                                                                      <  2:2042321/150‑1 (MQ=255)
                                                                  cTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCt                                                                                     >  1:2247851/1‑150 (MQ=255)
                                                                   tCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAATCCCCGGCTCGGTCCGGGGTTTTTGTCTGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTc                                                                                    >  2:3857460/1‑150 (MQ=255)
                                                                     aaTGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGa                                                                                  <  2:1009033/150‑1 (MQ=255)
                                                                     aaTGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGa                                                                                  <  2:5200726/150‑1 (MQ=255)
                                                                     aaTGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGa                                                                                  <  2:5200708/150‑1 (MQ=255)
                                                                     aaTGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGa                                                                                  >  2:1866769/1‑150 (MQ=255)
                                                                      aTGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAg                                                                                 >  2:3823414/1‑150 (MQ=255)
                                                                         cACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgc                                                                              >  2:1176848/1‑150 (MQ=255)
                                                                          aCTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgcg                                                                             >  2:1571343/1‑150 (MQ=255)
                                                                          aCTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgcg                                                                             <  1:4869395/150‑1 (MQ=255)
                                                                           cTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCTGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgcgc                                                                            <  2:2683450/150‑1 (MQ=255)
                                                                           cTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTgcgc                                                                            <  1:1323958/150‑1 (MQ=255)
                                                                               gTATGCGCATTGCCGCCACCAAGGGCCGTTCGCTGCTCTGATTGCGGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTGTTATGCCCGCGCAAGGGCGGGTCGTTCATCGTCGGTGGCTCACATCATCGCCCTCGAGGGCGCTggt                                                                        >  2:1142781/1‑148 (MQ=255)
                                                                               gTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGAt                                                                        <  1:1866769/150‑1 (MQ=255)
                                                                                tATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCTGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATc                                                                       <  1:3452403/150‑1 (MQ=255)
                                                                                    cgcATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGaaa                                                                   <  1:4858947/150‑1 (MQ=255)
                                                                                       aTTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCt                                                                <  2:3929002/150‑1 (MQ=255)
                                                                                         tGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTgc                                                              <  2:1067553/150‑1 (MQ=255)
                                                                                            cgccAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGca                                                           >  1:3885317/1‑150 (MQ=255)
                                                                                                  caaGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGa                                                     <  2:346052/150‑1 (MQ=255)
                                                                                                    aGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACg                                                   >  2:1192593/1‑150 (MQ=255)
                                                                                                     gggCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGt                                                  <  1:3850362/150‑1 (MQ=255)
                                                                                                     gggCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGt                                                  >  2:3009011/1‑150 (MQ=255)
                                                                                                         cGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGt                                              <  2:5419849/150‑1 (MQ=255)
                                                                                                          gTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTc                                             >  2:4522963/1‑150 (MQ=255)
                                                                                                           ttCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCCCATTGACGTCTGTcc                                            >  1:2006403/1‑150 (MQ=255)
                                                                                                           ttCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGATTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTcc                                            <  2:2550447/150‑1 (MQ=255)
                                                                                                              gctgctCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGa                                         >  2:4987335/1‑150 (MQ=255)
                                                                                                                   tctGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTGCAGATTACa                                    >  1:582711/1‑150 (MQ=255)
                                                                                                                   tctGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACa                                    >  1:194331/1‑150 (MQ=255)
                                                                                                                    ctGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaa                                   <  2:2708324/150‑1 (MQ=255)
                                                                                                                      gATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaaaa                                 <  1:1796824/150‑1 (MQ=255)
                                                                                                                      gATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaaaa                                 <  2:3085122/150‑1 (MQ=255)
                                                                                                                      gATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaaaa                                 >  2:4587141/1‑150 (MQ=255)
                                                                                                                      gATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACaaaa                                 >  2:5047397/1‑150 (MQ=255)
                                                                                                                      gATTGCTGGCGTGTGCGAAAACACCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTAAAATACGGTGGATCAAATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACAATTACGTCTTTCCACAGTACaaaa                                 >  1:1574355/1‑150 (MQ=255)
                                                                                                                          gCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAgcgc                             <  2:561492/150‑1 (MQ=255)
                                                                                                                            tGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCa                           <  2:2824152/150‑1 (MQ=255)
                                                                                                                                gtgtGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGAcg                       >  1:1740437/1‑150 (MQ=255)
                                                                                                                                   tgCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACACAAGCGCCAGAcgcgc                    <  1:1192593/150‑1 (MQ=255)
                                                                                                                                   tgCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGAcgcgc                    <  2:5523754/150‑1 (MQ=255)
                                                                                                                                     aaaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGCCTGTCCAGAGTACAAAAGCGCCAGACgcgcgc                  >  1:5239688/3‑150 (MQ=255)
                                                                                                                                      gAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAGAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCt                 <  2:2247851/150‑1 (MQ=255)
                                                                                                                                      aaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGa                                         >  1:4611656/2‑126 (MQ=255)
                                                                                                                                      aaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCt                 >  1:1565379/2‑150 (MQ=255)
                                                                                                                                      aaaaaCCCCGGACCGAGCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATGACAAAAGCGCCAGACgcgcgcg                 >  2:3628989/2‑149 (MQ=255)
                                                                                                                                        aaaCCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTgg               >  2:4490464/1‑150 (MQ=255)
                                                                                                                                         aaCCCCGGACCGAGCCGGGGTTTTCTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGc              >  2:932485/1‑150 (MQ=255)
                                                                                                                                           ccccGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCATATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTc            <  2:4111912/150‑1 (MQ=255)
                                                                                                                                            cccGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCg           >  1:1405197/1‑150 (MQ=255)
                                                                                                                                            cccGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATAGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACAGTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCg           >  2:1399674/1‑150 (MQ=255)
                                                                                                                                             ccGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCgg          >  2:3553862/1‑150 (MQ=255)
                                                                                                                                             ccGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCgg          >  1:3169851/1‑150 (MQ=255)
                                                                                                                                               ggCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGc        <  2:3223626/150‑1 (MQ=255)
                                                                                                                                               ggCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCGGGCTCGGGc        >  1:2104100/1‑150 (MQ=255)
                                                                                                                                                 cTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGACTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGCgg      >  2:4240563/1‑150 (MQ=255)
                                                                                                                                                  tCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGCGGt     >  2:2203997/1‑150 (MQ=255)
                                                                                                                                                  tCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGCGGt     <  1:2053098/150‑1 (MQ=255)
                                                                                                                                                    ggTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGCGGTTg   >  1:5080421/1‑150 (MQ=255)
                                                                                                                                                     gTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGCGGTTgg  <  2:1766586/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AACATCGCCACCTTTGAAGAAGATGGTGGTGCCTACAACCAGGCAGATGCTGCTGGCTTCATCAAGCTCAATGCACTGCGTATGCGCATTGCCGCCAACAAGGGCCGTTCGCTGCTCTGATTGCTGGCGTGTGCGAAAACCCCGGCTCGGTCCGGGGTTTTTTTATGCCCGCGCAAAGGCGGGTCGTTCAACTTCGGTGGCTCACATCATCGCCCTCGAGTGCGCTGATCGAAAGCTGCGCACATTGACGTCTGTCCAGATTACAAAAGCGCCAGACGCGCGCTGGCTCGGGCGGTTGG  >  NC_002947/1253023‑1253321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: