Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,359,686 C→A 30.2% E182* (GAA→TAA)  yjeH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,359,6860CA30.2% 578.2 / 223.8 307E182* (GAA→TAA) yjeHputative transporter
Reads supporting (aligned to +/- strand):  ref base C (110/104);  new base A (46/47);  total (156/151)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 308 total aligned reads displayed.

GACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAC  >  CP009273/4359537‑4359833
                                                                                                                                                     |                                                                                                                                                   
  cTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGa                                                                                                                                                   >  1:4530639/1‑150 (MQ=255)
     aTCCTGCCAGCAGCAGTCCAATCATCAAAGCTCGAGGATATTCACGCTCTTGATTTTTATATTCCGAGGCGAGTTGGGCAATTGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCgtttt                                                                                                                                                <  2:3107314/150‑5 (MQ=255)
       ccTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTa                                                                                                                                              <  2:165789/150‑1 (MQ=255)
        cTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAc                                                                                                                                             <  1:3810432/150‑1 (MQ=255)
        cTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAc                                                                                                                                             <  2:4186108/150‑1 (MQ=255)
        cTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAc                                                                                                                                             <  2:2988587/150‑1 (MQ=255)
         tGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAcc                                                                                                                                            <  1:2688309/150‑1 (MQ=255)
         tGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAcc                                                                                                                                            <  2:3977443/150‑1 (MQ=255)
         tGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAACTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAcc                                                                                                                                            <  1:5327878/150‑1 (MQ=255)
         tGCCAGCAGCAGACCAATAATCAAAGCACGAGGAAAATCAAGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTAcc                                                                                                                                            <  1:677744/150‑1 (MQ=255)
           ccAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAg                                                                                                                                          >  2:239789/1‑150 (MQ=255)
              gcagcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTg                                                                                                                                       <  2:4407893/150‑1 (MQ=255)
              gcagcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTg                                                                                                                                       <  2:770307/150‑1 (MQ=255)
               tagcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGc                                                                                                                                      >  2:1369752/2‑150 (MQ=255)
                 gcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCg                                                                                                                                    <  1:2834905/150‑1 (MQ=255)
                 gcagACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGATCGATATTACCAGGTGCCg                                                                                                                                    >  2:519104/1‑150 (MQ=255)
                   agACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGAGATTACCAGGTGCCGGa                                                                                                                                  >  1:3557818/1‑150 (MQ=255)
                    gACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGaa                                                                                                                                 >  1:1083923/1‑150 (MQ=255)
                       cAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAgg                                                                                                                              <  2:1020131/150‑1 (MQ=255)
                         atcatcAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAgggg                                                                                                                            >  2:1095859/1‑150 (MQ=255)
                          tcatcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGa                                                                                                                           <  2:1270071/150‑1 (MQ=255)
                           catcaAAGCACGAGGAAAATCACTCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGat                                                                                                                          <  2:1930567/150‑1 (MQ=255)
                            atcaAAGCACGATGAAAATAACTCTCTGGATTTTTAAATTCCGAGGAGATATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGata                                                                                                                         <  1:12401/150‑1 (MQ=255)
                            atcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGata                                                                                                                         <  2:2732021/150‑1 (MQ=255)
                             tcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGatat                                                                                                                        <  2:4476275/150‑1 (MQ=255)
                             tcaAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGatat                                                                                                                        <  1:4487162/150‑1 (MQ=255)
                               aaaGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTc                                                                                                                      >  1:2589109/1‑150 (MQ=255)
                                   cACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAg                                                                                                                  <  1:3017771/150‑1 (MQ=255)
                                      gAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGtt                                                                                                               >  2:5236162/1‑150 (MQ=255)
                                       aGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGttt                                                                                                              >  1:4917409/1‑150 (MQ=255)
                                          aaaaTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGat                                                                                                           >  2:5131767/1‑150 (MQ=255)
                                           aaaTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGata                                                                                                          >  1:1066872/1‑150 (MQ=255)
                                           aaaTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGata                                                                                                          >  2:3032637/1‑150 (MQ=255)
                                            aaTCTCGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGatat                                                                                                         <  2:3879496/150‑1 (MQ=255)
                                            aaTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGatat                                                                                                         <  1:208523/150‑1 (MQ=255)
                                             aTCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATc                                                                                                        >  1:4065510/1‑150 (MQ=255)
                                                    tctGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATcgcccgcc                                                                                                 <  1:1343100/150‑1 (MQ=255)
                                                     ctGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCgcccgccc                                                                                                <  2:4568982/150‑1 (MQ=255)
                                                     ctGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCgcccgccc                                                                                                >  2:4672450/1‑150 (MQ=255)
                                                        gATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCccacc                                                                                             >  2:4328107/1‑150 (MQ=255)
                                                            ttaaaaTTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGAt                                                                                         >  2:319558/4‑150 (MQ=255)
                                                                aaTTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCg                                                                                     <  1:3896131/150‑1 (MQ=255)
                                                                 aTTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGa                                                                                    >  1:3634325/1‑150 (MQ=255)
                                                                   tCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACa                                                                                  <  1:4751464/150‑1 (MQ=255)
                                                                   tCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACa                                                                                  >  1:1574080/1‑150 (MQ=255)
                                                                    ccGAGGCGAGATGGGCAAATGCCTCCAGACCGACAACACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACaa                                                                                 <  2:4686219/150‑1 (MQ=255)
                                                                      gAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATCACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATc                                                                               >  2:3382573/1‑150 (MQ=255)
                                                                           gagaTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAgcgc                                                                          >  1:4861272/1‑150 (MQ=255)
                                                                           gagaTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAgcgc                                                                          >  2:231615/1‑150 (MQ=255)
                                                                            agaTGGGCAAATGCCTACAGACCGACAAAACACCAGAACATCACAGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGACGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGcgcg                                                                         >  2:4642142/1‑150 (MQ=255)
                                                                             gaTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcga                                                                        >  2:3522994/1‑150 (MQ=255)
                                                                             gaTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcga                                                                        <  2:1173109/150‑1 (MQ=255)
                                                                              aTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcgac                                                                       <  1:944498/150‑1 (MQ=255)
                                                                              aTGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcgac                                                                       <  1:4623728/150‑1 (MQ=255)
                                                                               tGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcgacg                                                                      >  2:625057/1‑150 (MQ=255)
                                                                               tGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcgacg                                                                      <  1:2263353/150‑1 (MQ=255)
                                                                                gggCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcgacga                                                                     >  1:1479002/1‑150 (MQ=255)
                                                                                gggCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGcgacga                                                                     >  2:1972324/1‑150 (MQ=255)
                                                                                    aaaTGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAg                                                                 >  1:1355167/1‑150 (MQ=255)
                                                                                    aaaTGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAg                                                                 >  2:5055419/1‑150 (MQ=255)
                                                                                     aaggCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGt                                                                <  2:1311363/147‑1 (MQ=255)
                                                                                                aCCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAAc                                                     >  1:4373410/1‑150 (MQ=255)
                                                                                                 ccGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACg                                                    <  1:3386828/150‑1 (MQ=255)
                                                                                                  cGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACgg                                                   >  1:4943968/1‑150 (MQ=255)
                                                                                                   gACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGt                                                  >  1:2931148/1‑150 (MQ=255)
                                                                                                     cAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGttt                                                >  2:2075188/1‑150 (MQ=255)
                                                                                                     cAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGttt                                                >  2:2908726/1‑150 (MQ=255)
                                                                                                      aaaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTg                                               <  2:3685695/150‑1 (MQ=255)
                                                                                                      aaaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTg                                               >  2:1678309/1‑150 (MQ=255)
                                                                                                       aaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGt                                              <  2:4821895/150‑1 (MQ=255)
                                                                                                       aaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGt                                              >  1:1002048/1‑150 (MQ=255)
                                                                                                            ccAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGAtt                                         <  2:183209/150‑1 (MQ=255)
                                                                                                             cAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTa                                        <  2:3847735/150‑1 (MQ=255)
                                                                                                              aGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAg                                       >  2:1755665/1‑150 (MQ=255)
                                                                                                              aGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAg                                       <  1:1972324/150‑1 (MQ=255)
                                                                                                               gAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGc                                      <  1:2217340/150‑1 (MQ=255)
                                                                                                                    tCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTgg                                 <  1:1183220/150‑1 (MQ=255)
                                                                                                                      aCTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGaa                               >  1:1747722/1‑150 (MQ=255)
                                                                                                                      aCTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGaa                               >  1:1910971/1‑150 (MQ=255)
                                                                                                                      aCTGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGaa                               >  1:5141439/1‑150 (MQ=255)
                                                                                                                        tGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACt                             >  1:3111493/1‑150 (MQ=255)
                                                                                                                        tGATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACt                             <  2:2582731/150‑1 (MQ=255)
                                                                                                                         gATAACGCAGCAAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTg                            >  1:443384/1‑150 (MQ=255)
                                                                                                                              cGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCAcc                       <  1:1995813/150‑1 (MQ=255)
                                                                                                                               gcagcaAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAATGGGATATTCGCAGGTTTGATATCGCCCGCACACCAGATAGCGACAATCAGCGCGACTATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCAAcg                      >  1:1813745/1‑150 (MQ=255)
                                                                                                                                 agcaAATAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACcgcg                    >  1:3492130/1‑150 (MQ=255)
                                                                                                                                   caAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAg                  >  2:3063029/1‑150 (MQ=255)
                                                                                                                                    aaaTAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGt                 >  2:4978463/1‑150 (MQ=255)
                                                                                                                                    aaaTAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGt                 >  1:2526883/1‑150 (MQ=255)
                                                                                                                                     aaTAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTa                <  1:4726816/150‑1 (MQ=255)
                                                                                                                                     aaTAACCCGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTa                >  2:4144776/1‑150 (MQ=255)
                                                                                                                                       tAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTAcc              >  1:1616716/1‑150 (MQ=255)
                                                                                                                                        aaCCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCg             >  1:176733/1‑150 (MQ=255)
                                                                                                                                          cccGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGat           <  2:3912410/150‑1 (MQ=255)
                                                                                                                                          cccGGTAAGTTAGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACAAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGat           >  1:2893418/1‑150 (MQ=255)
                                                                                                                                            cGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGatat         <  1:758823/150‑1 (MQ=255)
                                                                                                                                              gTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCTGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATAc       >  2:3408325/1‑150 (MQ=255)
                                                                                                                                              gTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATAc       >  2:3718780/1‑150 (MQ=255)
                                                                                                                                                   ttCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCacac  >  2:4278213/1‑150 (MQ=255)
                                                                                                                                                   ttCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCacac  >  1:4225596/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAC  >  CP009273/4359537‑4359833

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: