Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,789,381 C→A 100% L15I (CTC→ATC)  waaC → ADP‑heptose:LPS heptosyl transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,789,3810CA98.1% 1053.4 / 6.0 332L15I (CTC→ATC) waaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base A (144/180);  total (149/183)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 336 total aligned reads displayed.

AAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGT  >  CP009273/3789234‑3789527
                                                                                                                                                   |                                                                                                                                                  
aaaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtatc                                                                                                                                                  <  1:3635891/150‑1 (MQ=255)
 aaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATcc                                                                                                                                                 >  2:1927391/1‑150 (MQ=255)
 aaGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATcc                                                                                                                                                 >  2:1330679/1‑150 (MQ=255)
    gTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATa                                                                                                                                              <  2:1859989/150‑1 (MQ=255)
     tGACGCTGCGGAGGTTTATCACCATTGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATAc                                                                                                                                             <  1:2288532/150‑1 (MQ=255)
     tGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATAc                                                                                                                                             <  2:4058432/150‑1 (MQ=255)
      gccGCTGCGGAGGGTTATCCCCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACg                                                                                                                                            <  2:1997238/148‑1 (MQ=255)
          cTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGc                                                                                                                                        >  1:200728/1‑150 (MQ=255)
           tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGcc                                                                                                                                       >  2:1097557/1‑150 (MQ=255)
              ggAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGc                                                                                                                                    >  2:2136413/1‑150 (MQ=255)
                aGGGTTATCACCTGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGcac                                                                                                                                  <  2:4231526/150‑1 (MQ=255)
                 tgtttATCACCAGAGATTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGTGTTTTGATCGTTAAAACATAGTCGATGGGCGATGTTATCCATACGTTGCCCGcact                                                                                                                                 <  1:1001732/147‑1 (MQ=255)
                 gggTTATCACCAGAGCTTACTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGcact                                                                                                                                 <  1:854382/150‑1 (MQ=255)
                  ggTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGcactc                                                                                                                                <  2:2766675/150‑1 (MQ=255)
                     tATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGcactcact                                                                                                                             <  2:2481825/150‑1 (MQ=255)
                      aTCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTg                                                                                                                            >  1:1588420/1‑150 (MQ=255)
                       tCACCAGATCTTAATCGAAATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATAGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGa                                                                                                                           >  2:711844/1‑150 (MQ=255)
                        ccccAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGAt                                                                                                                          <  2:123612/148‑1 (MQ=255)
                        cACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGAt                                                                                                                          <  2:189507/150‑1 (MQ=255)
                        cACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGAt                                                                                                                          <  2:1414368/150‑1 (MQ=255)
                            agagCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGccc                                                                                                                      <  1:1363587/150‑1 (MQ=255)
                             gagCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCca                                                                                                                     <  2:1718435/150‑1 (MQ=255)
                             gagCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCca                                                                                                                     <  2:2934299/150‑1 (MQ=255)
                                cTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCcagca                                                                                                                  >  2:2523815/1‑150 (MQ=255)
                                 ttAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCcagcag                                                                                                                 >  2:3829439/1‑150 (MQ=255)
                                   aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGc                                                                                                               >  1:5322331/1‑150 (MQ=255)
                                      cGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAAt                                                                                                            >  1:289825/1‑150 (MQ=255)
                                       gACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATc                                                                                                           >  1:3800377/1‑150 (MQ=255)
                                          aTTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGTAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCa                                                                                                        <  1:3326427/150‑1 (MQ=255)
                                          aTTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCa                                                                                                        <  1:1430343/150‑1 (MQ=255)
                                           ttACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAg                                                                                                       >  2:3969289/1‑150 (MQ=255)
                                            tACTCCCCATCGCGTACTGTAAGAACTCAACTCTCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCTATGTTATACATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCatg                                                                                                      <  1:1104116/150‑3 (MQ=255)
                                              cTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGa                                                                                                    >  1:4398212/1‑150 (MQ=255)
                                               tCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAt                                                                                                   <  2:5183810/150‑1 (MQ=255)
                                                ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                                  <  2:475379/150‑1 (MQ=255)
                                                ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                                  <  1:4173148/150‑1 (MQ=255)
                                                ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                                  <  1:194407/150‑1 (MQ=255)
                                                ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                                  <  1:5213703/150‑1 (MQ=255)
                                                ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                                  <  1:3134/150‑1 (MQ=255)
                                                 cccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTa                                                                                                 <  2:4442516/150‑1 (MQ=255)
                                                 cccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTa                                                                                                 <  1:2954651/150‑1 (MQ=255)
                                                   cAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAg                                                                                               <  2:3740357/150‑1 (MQ=255)
                                                    aGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGt                                                                                              >  2:4164642/1‑150 (MQ=255)
                                                          tACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACt                                                                                        <  1:371670/150‑1 (MQ=255)
                                                                aagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtgg                                                                                  <  1:2251029/150‑1 (MQ=255)
                                                                aagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtgg                                                                                  >  1:3033898/1‑150 (MQ=255)
                                                                aagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtgg                                                                                  <  2:2726548/150‑1 (MQ=255)
                                                                aagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtgg                                                                                  >  2:5069186/1‑150 (MQ=255)
                                                                 agaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTTATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtggt                                                                                 >  2:2210007/1‑150 (MQ=255)
                                                                   aaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGtggtgg                                                                               <  2:2416111/150‑1 (MQ=255)
                                                                    aCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGga                                                                              <  2:4700733/150‑1 (MQ=255)
                                                                    aCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGga                                                                              >  1:3411891/1‑150 (MQ=255)
                                                                     cTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaa                                                                             <  1:3980156/150‑1 (MQ=255)
                                                                      tCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaag                                                                            <  2:3537396/150‑1 (MQ=255)
                                                                       cAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaaga                                                                           >  2:564423/1‑150 (MQ=255)
                                                                       cAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaaga                                                                           >  2:4389534/1‑150 (MQ=255)
                                                                        aaCGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaagaa                                                                          <  2:2103160/150‑1 (MQ=255)
                                                                          cgcgCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAgg                                                                        <  1:4483090/150‑1 (MQ=255)
                                                                          cgcgCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAgg                                                                        <  2:3493571/150‑1 (MQ=255)
                                                                           gcgcTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAggg                                                                       <  1:3566737/150‑1 (MQ=255)
                                                                            cgcTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGt                                                                      <  1:4855264/150‑1 (MQ=255)
                                                                            cgcTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGt                                                                      <  1:2265558/150‑1 (MQ=255)
                                                                              cTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTc                                                                    <  1:4813086/150‑1 (MQ=255)
                                                                               tATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCg                                                                   <  1:3239322/150‑1 (MQ=255)
                                                                                  tGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGcac                                                                <  1:263957/150‑1 (MQ=255)
                                                                                     tACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGa                                                             >  2:3086416/1‑150 (MQ=255)
                                                                                      aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                            <  1:1371505/150‑1 (MQ=255)
                                                                                      aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                            <  2:1686009/150‑1 (MQ=255)
                                                                                      aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                            >  2:651850/1‑150 (MQ=255)
                                                                                      aCAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAt                                                            <  2:4782268/150‑1 (MQ=255)
                                                                                       cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAtt                                                           >  2:1479057/1‑150 (MQ=255)
                                                                                       cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAtt                                                           <  1:499971/150‑1 (MQ=255)
                                                                                       cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAtt                                                           >  2:5254383/1‑150 (MQ=255)
                                                                                       cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAtt                                                           <  1:2693275/150‑1 (MQ=255)
                                                                                          gagGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcct                                                        <  2:817/150‑1 (MQ=255)
                                                                                             gAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttcc                                                     <  1:1330679/150‑1 (MQ=255)
                                                                                             gAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttcc                                                     <  1:1927391/150‑1 (MQ=255)
                                                                                             gAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttcc                                                     <  1:1727281/150‑1 (MQ=255)
                                                                                              aaGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATtccttcct                                                    <  1:4704221/150‑1 (MQ=255)
                                                                                                   tGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGACGAAGGGTTCGCACAGATTCCTTCCTGGCAc                                               >  2:1708260/1‑150 (MQ=255)
                                                                                                       ggATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTg                                           >  1:4029385/1‑150 (MQ=255)
                                                                                                           gCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGt                                       >  1:5148138/1‑150 (MQ=255)
                                                                                                            cGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGtt                                      <  2:2263962/150‑1 (MQ=255)
                                                                                                             gggTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTg                                     >  2:3292649/1‑150 (MQ=255)
                                                                                                                   tGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAg                               <  1:3929756/150‑1 (MQ=255)
                                                                                                                   tGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAg                               >  1:3663919/1‑150 (MQ=255)
                                                                                                                         ttAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTc                         <  1:4735585/150‑1 (MQ=255)
                                                                                                                           aaaaCATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCt                       >  1:990193/1‑150 (MQ=255)
                                                                                                                                 tcgtcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCa                 <  1:3373359/150‑1 (MQ=255)
                                                                                                                                  cgtcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCaa                >  1:1840146/1‑150 (MQ=255)
                                                                                                                                  cgtcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCaa                >  1:4568196/1‑150 (MQ=255)
                                                                                                                                   gtcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCTGTTGAGCGAGTTATTCCTGTGGCAAt               >  1:4245018/1‑150 (MQ=255)
                                                                                                                                   gtcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt               >  1:2400833/1‑150 (MQ=255)
                                                                                                                                   gtcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt               <  2:3128295/150‑1 (MQ=255)
                                                                                                                                     cgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATAc             >  1:1096100/1‑150 (MQ=255)
                                                                                                                                      gATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACg            >  2:1302305/1‑150 (MQ=255)
                                                                                                                                        tGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTc          <  2:5017638/150‑1 (MQ=255)
                                                                                                                                           gCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCt       <  1:92082/150‑1 (MQ=255)
                                                                                                                                               tGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCg   >  1:2554507/1‑150 (MQ=255)
                                                                                                                                               tGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCg   >  1:4976359/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGGCGT  >  CP009273/3789234‑3789527

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: