Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,940,829 A→T 6.7% intergenic (+7/‑27) metZ → / → metW tRNA‑Met/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,940,8290AT 6.7% 605.5 / 2.9 219intergenic (+7/‑27)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (58/135);  new base T (7/11);  total (69/150)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

ATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAG  >  CP009273/2940680‑2940977
                                                                                                                                                     |                                                                                                                                                     
aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaaa                                                                                                                                                       <  2:1261712/150‑1 (MQ=255)
  ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                                                                     <  2:213943/150‑1 (MQ=255)
    gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                                                                                                   <  2:4822604/150‑1 (MQ=255)
      gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                                 <  1:2708623/150‑1 (MQ=255)
       aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATATtgatg                                                                                                                                                >  2:4217224/1‑150 (MQ=255)
        tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                               <  2:695296/150‑1 (MQ=255)
        tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                                               >  1:5089410/1‑150 (MQ=255)
         cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                                              <  2:1401474/150‑1 (MQ=255)
            ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑a                                                                                                                                          <  1:5168990/150‑1 (MQ=255)
            ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑a                                                                                                                                          <  1:5167355/150‑1 (MQ=255)
            ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑a                                                                                                                                          <  1:4748268/150‑1 (MQ=255)
             gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa                                                                                                                                         >  2:1617629/1‑150 (MQ=255)
              tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCAACTGAAAATATGATGAAGT‑aaa                                                                                                                                        >  1:2786995/1‑150 (MQ=21)
              tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                                                                        <  1:242641/150‑1 (MQ=255)
              tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                                                                        >  1:2118253/1‑150 (MQ=255)
              tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGATGAATT‑aaa                                                                                                                                        >  1:2118293/1‑150 (MQ=255)
               aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                                       <  2:2217186/150‑1 (MQ=255)
               aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                                       <  2:4492794/150‑1 (MQ=255)
                aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaagc                                                                                                                                     <  2:1914404/150‑4 (MQ=255)
                aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGc                                                                                                                                      <  1:1573913/150‑1 (MQ=255)
                 aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aatcac                                                                                                                                     >  2:4310540/1‑146 (MQ=255)
                  aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAg                                                                                                                                    <  1:866173/150‑1 (MQ=255)
                   aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGt                                                                                                                                   <  2:381694/150‑1 (MQ=255)
                   aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGt                                                                                                                                   <  2:2072892/150‑1 (MQ=255)
                      aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGTTGAAGT‑AAAGCAGTACg                                                                                                                                >  1:4570571/1‑150 (MQ=255)
                           gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAc                                                                                                                           <  1:2691834/150‑1 (MQ=255)
                            cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGaag                                                                                                                          >  1:3698710/1‑148 (MQ=255)
                            cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTATAAGAATA‑AAACCAGTAAGATGAcc                                                                                                                          >  2:3073244/1‑149 (MQ=18)
                              tGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcgcg                                                                                                                        <  1:2052606/150‑1 (MQ=255)
                              tGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                                                        <  2:3123018/150‑1 (MQ=255)
                                aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCggg                                                                                                                      <  2:329256/150‑1 (MQ=255)
                                     gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGa                                                                                                                 <  1:1657482/150‑1 (MQ=255)
                                       gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagc                                                                                                               >  1:4924025/1‑150 (MQ=255)
                                        aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagca                                                                                                              >  2:816595/1‑150 (MQ=255)
                                            gAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                                                          <  1:2437858/150‑1 (MQ=255)
                                              gTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAATT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg                                                                                                        >  2:269800/1‑150 (MQ=255)
                                              gTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg                                                                                                        <  1:4671688/150‑1 (MQ=255)
                                               tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                                                       >  1:3530738/1‑150 (MQ=255)
                                                 tcttGCGCATCCCCGGACGCTGGGTGGATCAGCCTTGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAg                                                                                                     <  2:3760119/147‑1 (MQ=17)
                                                  aGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTATGACCTAGA‑TGAGCTGTAATGTGACGCGGGGGGGGGCTGCCTGGTAGc                                                                                                    >  2:4092431/1‑150 (MQ=14)
                                                  aGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTAGCACCACGT‑AGTGCAGTCCGGTGACGCGGGATGATGCACCCTGGAAga                                                                                                    >  2:4751699/1‑149 (MQ=14)
                                                   gTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCTTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCTGCCCACGCAACCAATTAAAATTTGATGAATT‑ATAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                                                   <  1:2730954/150‑1 (MQ=17)
                                                    tGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtc                                                                                                  <  1:5349275/150‑1 (MQ=255)
                                                         aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGTCGCGGGGTGGAGCAGCATTGTAGCTCGTCgg                                                                                             >  1:3340000/1‑150 (MQ=255)
                                                          tCCACGGACGCGGGGTGGATCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCGGGTCGCTCGGCggg                                                                                            >  2:1178958/1‑150 (MQ=255)
                                                             aCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTc                                                                                         <  1:821777/150‑1 (MQ=255)
                                                             aCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTc                                                                                         <  1:3775277/150‑1 (MQ=255)
                                                               ggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCAt                                                                                       <  2:4563973/150‑1 (MQ=255)
                                                                 aCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATaa                                                                                     <  1:2759330/150‑1 (MQ=40)
                                                                    cgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAccc                                                                                  >  1:2572439/1‑150 (MQ=40)
                                                                      gggTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGa                                                                                <  2:2392467/150‑1 (MQ=255)
                                                                        gTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg                                                                              <  1:1781554/150‑1 (MQ=40)
                                                                        gTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg                                                                              <  1:4507228/150‑1 (MQ=40)
                                                                           gatcagaCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtc                                                                           <  2:518025/147‑1 (MQ=255)
                                                                               agcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGtcgtcg                                                                       <  2:1421540/150‑1 (MQ=40)
                                                                                   tGGTAGCTCGTCGGTCACTTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTc                                                                   <  2:5044530/150‑1 (MQ=255)
                                                                                     gTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑CAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCaa                                                                 <  1:92471/150‑1 (MQ=255)
                                                                                     gTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCaa                                                                 <  1:90328/150‑1 (MQ=40)
                                                                                       aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAAt                                                               <  2:2507651/150‑1 (MQ=255)
                                                                                        gCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATc                                                              <  2:2799572/150‑1 (MQ=40)
                                                                                         cTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATcc                                                             <  2:3709414/150‑1 (MQ=40)
                                                                                          tcgtcgGTCTCATAACCCGAAGGTCGTCGGTTCATCTCCGTCCCCCTCAACCAATAAAAATTTGCTGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCg                                                            <  2:3389008/150‑1 (MQ=11)
                                                                                          tcgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCg                                                            <  1:1617629/150‑1 (MQ=40)
                                                                                                ggCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAACTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc                                                      <  1:632665/150‑1 (MQ=38)
                                                                                                   tCATAACCCGAAGGTCGTCGGTTCACATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCa                                                   <  2:3472918/150‑1 (MQ=255)
                                                                                                     aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGA‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGACGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAc                                                 >  2:3138190/1‑150 (MQ=255)
                                                                                                       aaCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCa                                               >  1:3149763/1‑150 (MQ=40)
                                                                                                          ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACGAAATACATTTTGATGAATT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATc                                            <  2:3324679/150‑1 (MQ=14)
                                                                                                           cGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTTTATGAAGT‑AATGGAGTACGGTGTGGGGGGGTGGAGCAGCCTGGTTGCTCGACGGGCTCATAACACGAAGGGGGTAGGTTCAAATCCGGCCCCCGCAAACAATcg                                           >  2:5150904/1‑149 (MQ=14)
                                                                                                           cGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAGGAAGT‑CAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCa                                           >  2:963436/1‑150 (MQ=255)
                                                                                                              aGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AATGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAt                                        >  2:3452864/1‑150 (MQ=255)
                                                                                                              aGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAt                                        >  2:925952/1‑150 (MQ=255)
                                                                                                              aGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAt                                        <  1:137948/150‑1 (MQ=255)
                                                                                                                 tcgtcgGTTCAAATCCGGCCCCCGCAACCAATTACACTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtg                                     <  1:1477178/150‑1 (MQ=255)
                                                                                                                 tcgtcgGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtg                                     <  1:4203990/150‑1 (MQ=255)
                                                                                                                   gtcgGTTCAAATCCGGCCCCCGCAACCAATTAAATTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATACAAttagta                                   >  2:1271942/1‑146 (MQ=11)
                                                                                                                   gtcgGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATaaaagtacaa                                   >  1:2044577/1‑144 (MQ=255)
                                                                                                                    tcgGTTCAAATCCGGACCCCGCAACCAATTAAAATTTGATGAAGT‑CAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatg                                  <  2:5224670/150‑1 (MQ=255)
                                                                                                                     cgGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatga                                 >  2:1105553/1‑150 (MQ=255)
                                                                                                                      ggTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                <  1:4557417/150‑1 (MQ=255)
                                                                                                                        ttCAAATCCGGCCCCCTCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGGAACCAATCAAATTTGATGAAGt                              <  1:4310540/150‑1 (MQ=255)
                                                                                                                          cataTCCGGCCCGCTCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaa                            <  2:3090457/147‑1 (MQ=255)
                                                                                                                          cAAATCCGGCCCCCGCAAGCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGGGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaa                            <  1:4780971/150‑1 (MQ=255)
                                                                                                                           aaaTCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAATATGATGAAGTatg                           >  2:1514874/1‑148 (MQ=255)
                                                                                                                            taTCCGGCCCCCGCAACCAATTAATTTTTGATGTAGT‑ATAGGAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                          <  2:4924025/149‑1 (MQ=17)
                                                                                                                             atccttcccccGCAACCAATTATAATTTGTTGAATA‑TACGGATTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAg                         <  1:3994745/144‑1 (MQ=17)
                                                                                                                                  gaccccGCAACCAATTAAAATTGGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                    <  2:3090584/148‑1 (MQ=255)
                                                                                                                                  gCCCCCGCAACCAATTAAAAATTGATGAAGT‑AAAGCATTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                    <  2:333224/150‑1 (MQ=255)
                                                                                                                                       ctcAACCAATAATAATTTGATGAATT‑AAAGCAGTACTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTg               <  2:385972/148‑1 (MQ=255)
                                                                                                                                        tcAACCAATTAAATTTTGATGAATT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa              <  2:2118260/149‑1 (MQ=255)
                                                                                                                                         cAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCg                                                     >  1:3633008/1‑110 (MQ=255)
                                                                                                                                          aaCCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg            <  2:2572439/150‑1 (MQ=255)
                                                                                                                                          aaCCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGAGGAAGTAAAAGCAGTACGGTGAcg            >  1:136726/1‑150 (MQ=255)
                                                                                                                                            ccAATTATCATTCGTAGAATT‑CAAGCACTACTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg          <  2:3122858/150‑1 (MQ=14)
                                                                                                                                            ccAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCACTCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg          >  1:4544557/1‑150 (MQ=255)
                                                                                                                                              aaTTCAAATTTGATGAAGT‑AATGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCggg        <  2:5089410/150‑1 (MQ=255)
                                                                                                                                                 tAAAATTTGATGAAGT‑AAAGCAGTAAGGTGACGCGGGGTGGAGCAGCCTTGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTAAAATCCGGCCCCCGCAACCAAAAAAATTTCAAGAATAAAGAGAATTAAGTTGACGCTGGGTg     >  2:3561402/1‑150 (MQ=14)
                                                                                                                                                  ctaaTTTGATGACGT‑CAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAACTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTgg    <  2:2966915/148‑1 (MQ=255)
                                                                                                                                                   aaaTTTGATGAAGT‑AAAGCAGTACTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGa   <  1:1178958/150‑1 (MQ=255)
                                                                                                                                                    aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGag  >  1:2022805/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAG  >  CP009273/2940680‑2940977

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