Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,335,195 G→T 21.0% L61I (CTC→ATC)  yjdF ← DUF2238 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,335,1950GT21.0% 828.9 / 166.0 305L61I (CTC→ATC) yjdFDUF2238 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (133/108);  new base T (36/28);  total (169/136)
Fisher's exact test for biased strand distribution p-value = 8.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 322 total aligned reads displayed.

AGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGT  >  CP009273/4335047‑4335344
                                                                                                                                                    |                                                                                                                                                     
 gCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCAGTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagca                                                                                                                                                     <  2:3905210/150‑1 (MQ=255)
   cccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagcagc                                                                                                                                                   <  1:4115580/150‑1 (MQ=255)
     cTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAgcgg                                                                                                                                                 <  2:1895381/150‑1 (MQ=255)
     cTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAgcgg                                                                                                                                                 >  2:4783232/1‑150 (MQ=255)
       ggAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAgcggcg                                                                                                                                               >  1:4895391/1‑150 (MQ=255)
          aaaaGAGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTACGTGTATTGTCCGCCGACCATCAGGAAGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTa                                                                                                                                            >  2:3334943/1‑150 (MQ=255)
           aaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGACAATGAGCGTATAGAGCAGCGGCGTTaa                                                                                                                                           >  1:4769476/1‑150 (MQ=255)
              gTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACgg                                                                                                                                        <  1:4240299/150‑1 (MQ=255)
              gTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACgg                                                                                                                                        <  1:4579145/150‑1 (MQ=255)
                     aGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGt                                                                                                                                 <  1:2587018/150‑1 (MQ=255)
                     aGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGt                                                                                                                                 <  1:1149101/150‑1 (MQ=255)
                          gTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCACCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGc                                                                                                                            <  2:1595547/150‑1 (MQ=255)
                           tCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCg                                                                                                                           <  1:284188/150‑1 (MQ=255)
                           tCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCg                                                                                                                           <  2:2902864/150‑1 (MQ=255)
                            cATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCgg                                                                                                                          <  2:478592/150‑1 (MQ=255)
                            cATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCgg                                                                                                                          >  2:3504393/1‑150 (MQ=255)
                             aTACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGt                                                                                                                         >  2:2180637/1‑150 (MQ=255)
                              tACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTg                                                                                                                        <  2:416961/150‑1 (MQ=255)
                              tACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTg                                                                                                                        >  1:3703571/1‑150 (MQ=255)
                                cGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGc                                                                                                                      <  1:1848686/150‑1 (MQ=255)
                                   aTTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAg                                                                                                                   <  1:1738560/150‑1 (MQ=255)
                                    ttACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGc                                                                                                                  <  2:1840685/150‑1 (MQ=255)
                                     tACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCa                                                                                                                 <  1:2257083/150‑1 (MQ=255)
                                          cTCAACCCCAGCCAGTCCTGCACCTCAAAACCATCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGc                                                                                                            <  1:1590134/150‑1 (MQ=255)
                                           tCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCt                                                                                                           <  1:3204/150‑1 (MQ=255)
                                            cAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTg                                                                                                          <  1:1693444/150‑1 (MQ=255)
                                             aaCCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGc                                                                                                         >  1:4797280/1‑150 (MQ=255)
                                             aaCCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGc                                                                                                         >  2:1343732/1‑150 (MQ=255)
                                              aCCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCa                                                                                                        >  2:1525012/1‑150 (MQ=255)
                                               ccccAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCAc                                                                                                       <  1:79248/150‑1 (MQ=255)
                                                cccAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACg                                                                                                      >  2:622266/1‑150 (MQ=255)
                                                  cAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGAc                                                                                                    <  1:1497506/150‑1 (MQ=255)
                                                   aGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACa                                                                                                   >  1:2835993/1‑150 (MQ=255)
                                                     ccATTCCTGCACCTCAAAACCACCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACaat                                                                                                 <  1:1000207/150‑1 (MQ=255)
                                                         tCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATc                                                                                             <  1:3882021/150‑1 (MQ=255)
                                                           cTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCAc                                                                                           >  2:3092582/1‑150 (MQ=255)
                                                            tGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGACGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCAcc                                                                                          >  2:1105491/1‑150 (MQ=255)
                                                              cACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCgg                                                                                        >  2:95991/1‑150 (MQ=255)
                                                                ccTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGACAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgtg                                                                                      >  1:3016880/1‑150 (MQ=255)
                                                                ccTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgtg                                                                                      >  2:101147/1‑150 (MQ=255)
                                                                 cTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgtgt                                                                                     >  1:2172231/1‑150 (MQ=255)
                                                                   cAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCa                                                                                   <  1:918683/150‑1 (MQ=255)
                                                                     aaaCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACt                                                                                 >  1:1054318/1‑150 (MQ=255)
                                                                      aaCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACtt                                                                                >  1:4977351/1‑150 (MQ=255)
                                                                       aCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTc                                                                               <  1:3113526/150‑1 (MQ=255)
                                                                       aCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTc                                                                               <  1:3602626/150‑1 (MQ=255)
                                                                        ccAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTcc                                                                              >  1:1997148/1‑150 (MQ=255)
                                                                          aaCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCAt                                                                            <  1:4279077/150‑1 (MQ=255)
                                                                            cGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCa                                                                          >  2:1875686/1‑150 (MQ=255)
                                                                            cGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCa                                                                          <  2:4423312/150‑1 (MQ=255)
                                                                            cGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCa                                                                          >  1:4711121/1‑150 (MQ=255)
                                                                            cGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCa                                                                          >  1:2959356/1‑150 (MQ=255)
                                                                               ggACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGcc                                                                       >  1:165114/1‑150 (MQ=255)
                                                                                gACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCa                                                                      <  2:142903/150‑1 (MQ=255)
                                                                                   tttCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGt                                                                   >  2:714459/1‑150 (MQ=255)
                                                                                   tttCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGt                                                                   >  1:3753611/1‑150 (MQ=255)
                                                                                      cgcgTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAg                                                                >  1:3640660/1‑150 (MQ=255)
                                                                                         gTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGttt                                                             <  1:4667505/150‑1 (MQ=255)
                                                                                           aGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGt                                                           >  1:3367435/1‑150 (MQ=255)
                                                                                             gtgtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTAACTTCCATCAGCCAGGTGAGTTTGTCa                                                         <  2:3967290/150‑1 (MQ=255)
                                                                                              tgtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCAt                                                        <  1:2936844/150‑1 (MQ=255)
                                                                                              tgtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCAt                                                        <  1:564638/150‑1 (MQ=255)
                                                                                                tATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATgg                                                      <  1:2247665/150‑1 (MQ=255)
                                                                                                  ttGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGACTGGCGGTGGCAAGCACTAGCTGCACGACAATAATCACCGGTGTCACTTCCAGCAGCCAGGTGAGTGTGACATGGGc                                                    >  1:3164090/1‑150 (MQ=255)
                                                                                                  ttGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGc                                                    <  1:1101029/150‑1 (MQ=255)
                                                                                                   tGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGcc                                                   <  2:5016215/150‑1 (MQ=255)
                                                                                                   tGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGGGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGcc                                                   <  2:4765239/150‑1 (MQ=255)
                                                                                                    gTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCg                                                  <  2:4661552/150‑1 (MQ=255)
                                                                                                       cgccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGaaa                                               >  2:3532605/1‑150 (MQ=255)
                                                                                                        gccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAAt                                              >  2:4572232/1‑150 (MQ=255)
                                                                                                        gccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAAt                                              <  1:643475/150‑1 (MQ=255)
                                                                                                        gccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAAt                                              >  1:4448110/1‑150 (MQ=255)
                                                                                                           gACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTAACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGcc                                           >  1:407294/1‑150 (MQ=255)
                                                                                                            aCCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCg                                          >  2:3484095/1‑150 (MQ=255)
                                                                                                             ccATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCgg                                         <  2:3612233/150‑1 (MQ=255)
                                                                                                                 cAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGtata                                     <  2:2505286/150‑1 (MQ=255)
                                                                                                                  aGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATaa                                    >  1:2083962/1‑150 (MQ=255)
                                                                                                                    gatgatGGCGTTAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCTGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTAATGGGCCTAAATGCCGGTATAAAt                                  <  2:4830837/150‑1 (MQ=255)
                                                                                                                     atgatgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATc                                 <  2:754553/150‑1 (MQ=255)
                                                                                                                      tgatgGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCa                                <  2:2955306/150‑1 (MQ=255)
                                                                                                                      tgatgGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCa                                <  2:2925715/150‑1 (MQ=255)
                                                                                                                       gatgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTCTAAATCAg                               >  2:2447333/1‑150 (MQ=255)
                                                                                                                           gCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGAtt                           <  2:3986023/150‑1 (MQ=255)
                                                                                                                           gCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGAtt                           >  1:2542772/1‑150 (MQ=255)
                                                                                                                           gCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGAtt                           <  1:1315186/150‑1 (MQ=255)
                                                                                                                              tGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAAc                        <  2:4081811/150‑1 (MQ=255)
                                                                                                                              tGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAAc                        <  2:922038/150‑1 (MQ=255)
                                                                                                                              tGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAAc                        <  2:4932745/150‑1 (MQ=255)
                                                                                                                                aaGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGt                      >  2:821200/1‑150 (MQ=255)
                                                                                                                                 aGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGtt                     <  1:4338756/150‑1 (MQ=255)
                                                                                                                                    aaaaTGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGc                  <  1:527098/150‑1 (MQ=255)
                                                                                                                                     aaaTGAGCGTATAGAGCAGCGGCGTTAAAGGATAACGTCTGGAGGTGGCAAGCAGTAGCTGCACGACAATAATCACAGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCg                 >  2:1189666/1‑150 (MQ=255)
                                                                                                                                        tGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCa              <  2:1692374/150‑1 (MQ=255)
                                                                                                                                        tGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCa              <  2:3845062/150‑1 (MQ=255)
                                                                                                                                        tGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCa              <  1:2331720/150‑1 (MQ=255)
                                                                                                                                           gCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTg           <  2:1636840/150‑1 (MQ=255)
                                                                                                                                              tataGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcgc        >  2:2790488/1‑150 (MQ=255)
                                                                                                                                              tataGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcgc        >  2:968543/1‑150 (MQ=255)
                                                                                                                                                taGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTg      >  2:2085358/1‑150 (MQ=255)
                                                                                                                                                    gcagcGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGgtgt  >  2:231207/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGT  >  CP009273/4335047‑4335344

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: