Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,335,195 G→T 18.9% 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,1950GT18.9% 881.5 / 149.9 307L61I (CTC→ATC) yjdFDUF2238 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (114/135);  new base T (32/26);  total (146/161)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 321 total aligned reads displayed.

CAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTG  >  CP009273/4335046‑4335343
                                                                                                                                                     |                                                                                                                                                    
cAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATagag                                                                                                                                                      <  2:3932753/150‑1 (MQ=255)
   ccccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGagcag                                                                                                                                                   >  1:1253464/1‑150 (MQ=255)
       tGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAgcggc                                                                                                                                               >  1:4099453/1‑150 (MQ=255)
         gAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGt                                                                                                                                             >  2:922018/1‑150 (MQ=255)
         gAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGt                                                                                                                                             <  2:429632/150‑1 (MQ=255)
          aaaaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGtt                                                                                                                                            <  1:1967206/150‑1 (MQ=255)
             aaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAAc                                                                                                                                         <  2:40483/150‑1 (MQ=255)
             aaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAAc                                                                                                                                         >  2:1465667/1‑150 (MQ=255)
                 gCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGAt                                                                                                                                     <  2:1154033/150‑1 (MQ=255)
                 gCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGAt                                                                                                                                     >  2:1825351/1‑150 (MQ=255)
                  ccccAGCTTGTCATACGGATTACGACTCAACCCCCGCCATTCCTGCACCGCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGGGACGGAAAATGAGCGGAGAGAGCAGCGGCGTTAACGGATa                                                                                                                                    >  2:3392367/1‑150 (MQ=255)
                  ccccAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATa                                                                                                                                    >  2:3443822/1‑150 (MQ=255)
                   cccAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAATGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATaa                                                                                                                                   >  1:2450491/1‑150 (MQ=255)
                   cccAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATaa                                                                                                                                   <  1:723036/150‑1 (MQ=255)
                   cccAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATaa                                                                                                                                   <  2:1235200/150‑1 (MQ=255)
                   cccAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATaa                                                                                                                                   >  1:358039/1‑150 (MQ=255)
                            tCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACGTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCg                                                                                                                          >  2:4745556/1‑150 (MQ=255)
                               tACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTg                                                                                                                       >  2:3015554/1‑150 (MQ=255)
                               tACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTg                                                                                                                       >  2:2281875/1‑150 (MQ=255)
                                 cGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGTGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGc                                                                                                                     <  1:2596922/150‑1 (MQ=255)
                                   gATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCaa                                                                                                                   <  1:3509800/150‑1 (MQ=255)
                                   gATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCaa                                                                                                                   <  1:372241/150‑1 (MQ=255)
                                   gATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCaa                                                                                                                   <  2:3762890/150‑1 (MQ=255)
                                   gATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCaa                                                                                                                   <  1:605804/150‑1 (MQ=255)
                                      tACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCa                                                                                                                <  2:2223647/150‑1 (MQ=255)
                                      tACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCa                                                                                                                >  1:541452/1‑150 (MQ=255)
                                          aCTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAg                                                                                                            >  2:4345398/1‑150 (MQ=255)
                                          aCTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAg                                                                                                            <  1:2265509/150‑1 (MQ=255)
                                          aCTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAg                                                                                                            >  2:572470/1‑150 (MQ=255)
                                           cTCAACCCCAGCCAGTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGc                                                                                                           <  1:3695946/150‑1 (MQ=255)
                                              aaCCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGc                                                                                                        <  1:630176/150‑1 (MQ=255)
                                               tccccAGCCTGTCCTGCACCTCATTACCATCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGct                                                                                                       <  2:3141742/149‑2 (MQ=255)
                                                   cAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGAc                                                                                                   <  1:1066109/150‑1 (MQ=255)
                                                    aGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACa                                                                                                  <  2:3906794/150‑1 (MQ=255)
                                                     gCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACaa                                                                                                 >  2:4230413/1‑150 (MQ=255)
                                                            cTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCAc                                                                                          <  2:1396804/150‑1 (MQ=255)
                                                              gCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCg                                                                                        <  1:4103752/150‑1 (MQ=255)
                                                               cACCTCAAAACCCACGGGGACTTTCGCGGAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCgg                                                                                       <  2:3647290/150‑1 (MQ=255)
                                                                aCCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgt                                                                                      <  1:4329520/150‑1 (MQ=255)
                                                                aCCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgt                                                                                      <  1:86869/150‑1 (MQ=255)
                                                                 ccTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgtg                                                                                     >  1:4792091/1‑150 (MQ=255)
                                                                   tCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTc                                                                                   <  2:876491/150‑1 (MQ=255)
                                                                       aaCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACtt                                                                               <  2:454481/150‑1 (MQ=255)
                                                                        aCCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTc                                                                              <  1:2882148/150‑1 (MQ=255)
                                                                         ccAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTcc                                                                             >  1:1221269/1‑150 (MQ=255)
                                                                          cAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCa                                                                            >  2:4441002/1‑150 (MQ=255)
                                                                          cAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCa                                                                            <  1:4307889/150‑1 (MQ=255)
                                                                           aaCGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCAt                                                                           <  2:3773306/150‑1 (MQ=255)
                                                                             cGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCa                                                                         >  1:3042854/1‑150 (MQ=255)
                                                                                 gACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCa                                                                     >  2:1796742/1‑150 (MQ=255)
                                                                                     ttCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTg                                                                 <  2:4635797/150‑1 (MQ=255)
                                                                                     ttCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTg                                                                 >  2:1318534/1‑150 (MQ=255)
                                                                                     ttCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTg                                                                 <  2:2191181/150‑1 (MQ=255)
                                                                                      tCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGa                                                                >  2:2935723/1‑150 (MQ=255)
                                                                                         cgTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGtt                                                             >  2:1159315/1‑150 (MQ=255)
                                                                                           tAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTg                                                           >  1:341661/1‑150 (MQ=255)
                                                                                            aGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGt                                                          >  1:3775293/1‑150 (MQ=255)
                                                                                            aGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGt                                                          <  2:2831378/150‑1 (MQ=255)
                                                                                              gtgtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCa                                                        <  2:815522/150‑1 (MQ=255)
                                                                                                 tATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATgg                                                     >  1:103611/1‑150 (MQ=255)
                                                                                                  aTTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGTGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATggg                                                    >  2:2335265/1‑150 (MQ=255)
                                                                                                  aTTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATggg                                                    >  1:4559419/1‑150 (MQ=255)
                                                                                                  aTTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATggg                                                    <  2:4761092/150‑1 (MQ=255)
                                                                                                  aTTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATggg                                                    <  2:4188857/150‑1 (MQ=255)
                                                                                                  aTTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATggg                                                    >  2:1033990/1‑150 (MQ=255)
                                                                                                      tCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGa                                                >  1:1494541/1‑150 (MQ=255)
                                                                                                       ccgccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGaa                                               <  2:710055/150‑1 (MQ=255)
                                                                                                          ccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATg                                            >  1:3889543/1‑150 (MQ=255)
                                                                                                          ccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATg                                            <  1:3195716/150‑1 (MQ=255)
                                                                                                             aCCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCg                                         <  1:3939023/150‑1 (MQ=255)
                                                                                                             aCCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCg                                         >  1:2254177/1‑150 (MQ=255)
                                                                                                              ccAGCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCgg                                        >  1:2420428/1‑150 (MQ=255)
                                                                                                                aTCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGta                                      <  2:1594522/150‑1 (MQ=255)
                                                                                                                aTCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGta                                      <  1:3055370/150‑1 (MQ=255)
                                                                                                                   aGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATaa                                   >  2:413121/1‑150 (MQ=255)
                                                                                                                      atgatgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATc                                <  1:815152/150‑1 (MQ=255)
                                                                                                                       tgatgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCa                               >  2:337941/1‑150 (MQ=255)
                                                                                                                         atgGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTAACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAgg                             >  1:4753025/1‑150 (MQ=255)
                                                                                                                             cGTGAAGGAAAATGAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTa                         <  1:1825351/150‑1 (MQ=255)
                                                                                                                              gTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTaa                        <  1:1689541/150‑1 (MQ=255)
                                                                                                                              gTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTaa                        <  1:1661003/150‑1 (MQ=255)
                                                                                                                               tGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAAc                       <  1:2117616/150‑1 (MQ=255)
                                                                                                                                  aGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGtt                    <  1:4128135/150‑1 (MQ=255)
                                                                                                                                   ggAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTa                   <  1:2666010/150‑1 (MQ=255)
                                                                                                                                    gAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAg                  >  1:4538034/1‑150 (MQ=255)
                                                                                                                                    gAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAg                  <  2:2486253/150‑1 (MQ=255)
                                                                                                                                    gAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAg                  <  2:612727/150‑1 (MQ=255)
                                                                                                                                       aaTGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGt               >  1:4738719/1‑150 (MQ=255)
                                                                                                                                         tGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCa             >  1:2748439/1‑150 (MQ=255)
                                                                                                                                          gAGCGTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAg            >  2:2063374/1‑150 (MQ=255)
                                                                                                                                             cGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgc         <  1:1931951/150‑1 (MQ=255)
                                                                                                                                              gTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcg        >  1:1286013/1‑150 (MQ=255)
                                                                                                                                              gTATAGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcg        >  2:3217255/1‑150 (MQ=255)
                                                                                                                                              gTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcg        <  2:2494320/150‑1 (MQ=255)
                                                                                                                                               tataGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgcgc       <  1:3301921/150‑1 (MQ=255)
                                                                                                                                                ataGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCt      >  1:3420027/1‑150 (MQ=255)
                                                                                                                                                ataGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCt      >  2:230046/1‑150 (MQ=255)
                                                                                                                                                 taGATCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTg     <  2:1730014/150‑1 (MQ=255)
                                                                                                                                                   gagCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGgt   <  1:2325499/150‑1 (MQ=255)
                                                                                                                                                    agcagcGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACACTAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGgtg  <  2:2623012/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTG  >  CP009273/4335046‑4335343

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: