Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,454,297 T→G 100% E39D (GAA→GAC yfcZ ← UPF0381 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,454,2970TG99.2% 910.2 / ‑4.8 224E39D (GAA→GACyfcZUPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (106/116);  total (107/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 229 total aligned reads displayed.

CAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCAA  >  CP009273/2454153‑2454441
                                                                                                                                                |                                                                                                                                                
cAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCt                                                                                                                                               >  2:2689786/1‑148 (MQ=255)
cAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCt                                                                                                                                               >  1:3062782/1‑148 (MQ=255)
cAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCt                                                                                                                                               >  1:2234592/1‑148 (MQ=255)
   cATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCtctgcgtctgcg                                                                                                                                            >  1:536366/1‑148 (MQ=255)
   cATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCtctgcgtctgcg                                                                                                                                            >  1:2784177/1‑148 (MQ=255)
      ttCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCgg                                                                                                                                         >  2:857043/1‑148 (MQ=255)
        cAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGACTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGtt                                                                                                                                       >  1:1063714/1‑148 (MQ=255)
         aGCTTTGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGttt                                                                                                                                      <  2:1369157/148‑1 (MQ=255)
                cAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAaca                                                                                                                               >  2:1134750/1‑148 (MQ=255)
                  ggCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACAcg                                                                                                                             >  2:2071743/1‑148 (MQ=255)
                  ggCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACAcg                                                                                                                             <  1:1772703/148‑1 (MQ=255)
                   gCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACAcgc                                                                                                                            >  1:1576534/1‑148 (MQ=255)
                    cAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACAcgcg                                                                                                                           >  2:1189397/1‑148 (MQ=255)
                     aaaaGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCgg                                                                                                                          <  1:3309071/148‑1 (MQ=255)
                          tGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTa                                                                                                                     >  2:871156/1‑148 (MQ=255)
                          tGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTa                                                                                                                     >  2:2693918/1‑148 (MQ=255)
                           gAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTaa                                                                                                                    <  1:2914244/148‑1 (MQ=255)
                            aaaTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAg                                                                                                                   >  1:2872440/1‑148 (MQ=255)
                             aaTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGa                                                                                                                  <  1:761150/148‑1 (MQ=255)
                             aaTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGa                                                                                                                  <  1:1282807/148‑1 (MQ=255)
                             aaTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGGCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGa                                                                                                                  >  1:289158/1‑148 (MQ=255)
                              aTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGAc                                                                                                                 >  1:1898914/1‑148 (MQ=255)
                              aTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGAc                                                                                                                 >  2:2929954/1‑148 (MQ=255)
                                  atatCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCgg                                                                                                             <  1:699712/148‑1 (MQ=255)
                                   tatCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGt                                                                                                            >  1:1705470/1‑148 (MQ=255)
                                    atCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTg                                                                                                           <  1:179063/148‑1 (MQ=255)
                                    atCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTg                                                                                                           >  1:1442904/1‑148 (MQ=255)
                                     tCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGc                                                                                                          >  1:2690962/1‑148 (MQ=255)
                                     tCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGc                                                                                                          <  1:2420757/148‑1 (MQ=255)
                                     tCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTGTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCAGTAAGACGCGGTGc                                                                                                          <  2:1988002/148‑1 (MQ=255)
                                       cAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAg                                                                                                        <  1:3206943/148‑1 (MQ=255)
                                        aGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGt                                                                                                       <  2:2806795/148‑1 (MQ=255)
                                         ggCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTc                                                                                                      >  1:3085794/1‑148 (MQ=255)
                                         ggCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTc                                                                                                      >  2:2524620/1‑148 (MQ=255)
                                          gcgcACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCg                                                                                                     <  1:2773467/148‑1 (MQ=255)
                                           cgcACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCgg                                                                                                    <  1:1399313/148‑1 (MQ=255)
                                           cgcACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGGTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCgg                                                                                                    >  1:3157695/1‑148 (MQ=255)
                                              acacCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGt                                                                                                 <  2:431953/148‑1 (MQ=255)
                                              acacCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGt                                                                                                 <  2:432991/148‑1 (MQ=255)
                                               cacCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGtt                                                                                                >  2:1050285/1‑148 (MQ=255)
                                                acCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTg                                                                                               <  1:2621293/148‑1 (MQ=255)
                                                 ccGTCAGACTCTTCAGTGAAAGTCGTGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGt                                                                                              <  1:1003746/148‑1 (MQ=255)
                                                 ccGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGt                                                                                              <  2:2872440/148‑1 (MQ=255)
                                                  cGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTc                                                                                             <  1:2368172/148‑1 (MQ=255)
                                                     cAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCAt                                                                                          <  2:240852/148‑1 (MQ=255)
                                                       gACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATaa                                                                                        >  2:1953155/1‑148 (MQ=255)
                                                        aCTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAAt                                                                                       >  2:3072848/1‑148 (MQ=255)
                                                          tctTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATgg                                                                                     >  2:1434736/1‑148 (MQ=255)
                                                                   aaaGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACAt                                                                            <  1:1110631/148‑1 (MQ=255)
                                                                     aGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATcc                                                                          >  2:640330/1‑148 (MQ=255)
                                                                       tCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCAt                                                                        >  1:2250137/1‑148 (MQ=255)
                                                                       tCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCAt                                                                        <  1:1300402/148‑1 (MQ=255)
                                                                       tCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCAt                                                                        <  2:1689728/148‑1 (MQ=255)
                                                                         ggggTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATAc                                                                      <  1:1374713/148‑1 (MQ=255)
                                                                                 tttGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCaa                                                              >  1:1608912/1‑148 (MQ=255)
                                                                                   tGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAAc                                                            <  1:2308421/148‑1 (MQ=255)
                                                                                     cAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCg                                                          >  1:578779/1‑148 (MQ=255)
                                                                                     cAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCg                                                          <  2:3085794/148‑1 (MQ=255)
                                                                                     cAAGGTTCAGATACCACGCTACGGGCTTTTTCAGTCAGCGCAGCCAGCGTTTGCTATGAGTATTCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATAAATACAGCAGCAAACCg                                                          >  2:1317656/1‑148 (MQ=255)
                                                                                        ggTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCgggg                                                       <  2:233010/148‑1 (MQ=255)
                                                                                        ggTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCgggg                                                       <  1:2686715/148‑1 (MQ=255)
                                                                                          ttCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGtt                                                     <  1:35906/148‑1 (MQ=255)
                                                                                           tCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGttt                                                    <  2:1484971/148‑1 (MQ=255)
                                                                                           tCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGttt                                                    <  1:1430480/148‑1 (MQ=255)
                                                                                           tCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGttt                                                    >  2:1122008/1‑148 (MQ=255)
                                                                                           tCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGttt                                                    <  2:2688193/148‑1 (MQ=255)
                                                                                           tCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGttt                                                    >  2:2209288/1‑148 (MQ=255)
                                                                                            cAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtc                                                   <  1:1965692/148‑1 (MQ=255)
                                                                                            cAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtc                                                   <  1:1767754/148‑1 (MQ=255)
                                                                                                ttCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcatca                                               >  1:1556132/1‑148 (MQ=255)
                                                                                                ttCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcatca                                               <  2:2935260/148‑1 (MQ=255)
                                                                                                ttCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcatca                                               >  1:3048812/1‑148 (MQ=255)
                                                                                                ttCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcatca                                               >  1:3048749/1‑148 (MQ=255)
                                                                                                ttCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcatca                                               <  2:1993584/148‑1 (MQ=255)
                                                                                                 tCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAg                                              >  2:557052/1‑148 (MQ=255)
                                                                                                  ccACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCTGGGTTTCATCAGc                                             <  2:309887/148‑1 (MQ=255)
                                                                                                        tACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACa                                       >  2:3003265/1‑148 (MQ=255)
                                                                                                        tACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACa                                       <  1:1337554/148‑1 (MQ=255)
                                                                                                        tACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACa                                       <  2:2881954/148‑1 (MQ=255)
                                                                                                         aCGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACAt                                      <  2:463491/148‑1 (MQ=255)
                                                                                                         aCGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACAt                                      >  2:2570973/1‑148 (MQ=255)
                                                                                                            ggCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTa                                   <  2:2884431/148‑1 (MQ=255)
                                                                                                               tttttCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTc                                <  1:3132420/148‑1 (MQ=255)
                                                                                                               tttttCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTc                                <  1:1189403/148‑1 (MQ=255)
                                                                                                                 tttCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCAt                              <  2:1412168/148‑1 (MQ=255)
                                                                                                                  ttCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATc                             <  1:114001/148‑1 (MQ=255)
                                                                                                                   tcagtcagCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCg                            <  2:776802/148‑1 (MQ=255)
                                                                                                                    cagtcagCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGt                           <  2:2494327/148‑1 (MQ=255)
                                                                                                                       tcagCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGa                        >  2:1724549/1‑148 (MQ=255)
                                                                                                                        cagCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGAt                       >  1:3185591/1‑148 (MQ=255)
                                                                                                                           cgccgcCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTc                    >  2:751477/1‑148 (MQ=255)
                                                                                                                           cgccgcCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTc                    <  2:1202978/148‑1 (MQ=255)
                                                                                                                                 cAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGt              >  1:2354564/1‑148 (MQ=255)
                                                                                                                                       ttGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGc        >  2:3282668/1‑148 (MQ=255)
                                                                                                                                       ttGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGc        >  1:143629/1‑148 (MQ=255)
                                                                                                                                       ttGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGc        >  2:2663005/1‑148 (MQ=255)
                                                                                                                                          ctctGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCAcc     >  1:3239154/1‑148 (MQ=255)
                                                                                                                                             tGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCaa  <  1:1187237/148‑1 (MQ=255)
                                                                                                                                             tGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCaa  <  1:248809/148‑1 (MQ=255)
                                                                                                                                             tGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCaa  <  2:1775996/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCAA  >  CP009273/2454153‑2454441

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: