Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,048,928 C→A 32.1% E903* (GAA→TAA)  torS ← hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,048,9280CA32.1% 440.9 / 248.8 261E903* (GAA→TAA) torShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base C (74/103);  new base A (41/43);  total (115/146)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 264 total aligned reads displayed.

CTGGCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTC  >  CP009273/1048779‑1049076
                                                                                                                                                     |                                                                                                                                                    
  ggCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCa                                                                                                                                                    >  1:1399000/1‑150 (MQ=255)
  ggCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCa                                                                                                                                                    <  1:1513153/150‑1 (MQ=255)
   gccgccGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCag                                                                                                                                                   >  2:1574060/1‑150 (MQ=255)
   gccgccGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACag                                                                                                                                                   >  1:3245282/1‑150 (MQ=255)
    ccgccgCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACaga                                                                                                                                                  >  2:267143/1‑150 (MQ=255)
         gcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGc                                                                                                                                             <  2:3196484/150‑1 (MQ=255)
         gcTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGc                                                                                                                                             >  1:2891381/1‑150 (MQ=255)
          cggCCAGCGGGTTCTTTTATGCTCTGCCATGATTTTGCCAGTCCGTTAAGTCGTATACCACATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTCTGCAACCATGCTTCCAGAGCGGCa                                                                                                                                            <  2:4253552/148‑1 (MQ=255)
          cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCa                                                                                                                                            <  1:480956/150‑1 (MQ=255)
          cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCa                                                                                                                                            <  1:360276/150‑1 (MQ=255)
          cTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCa                                                                                                                                            <  1:5028851/150‑1 (MQ=255)
           tGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCaa                                                                                                                                           >  1:75414/1‑150 (MQ=255)
           tGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCaa                                                                                                                                           >  2:3121112/1‑150 (MQ=255)
           tGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCaa                                                                                                                                           <  2:3366651/150‑1 (MQ=255)
              cAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAacac                                                                                                                                        >  2:4866430/1‑150 (MQ=255)
                 cGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGc                                                                                                                                     >  1:2110260/1‑150 (MQ=255)
                   ggTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgt                                                                                                                                   >  2:222840/1‑150 (MQ=255)
                    gTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCgtg                                                                                                                                  <  2:491476/150‑1 (MQ=255)
                     ttCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCgtgt                                                                                                                                 <  1:4239576/150‑1 (MQ=255)
                     ttCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCgtgt                                                                                                                                 <  2:949574/150‑1 (MQ=255)
                          ttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAAttt                                                                                                                            >  1:1258587/1‑150 (MQ=255)
                          ttATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAAttt                                                                                                                            >  2:5038109/1‑150 (MQ=255)
                           tATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATttc                                                                                                                           >  2:3400106/1‑150 (MQ=255)
                           tATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATttc                                                                                                                           <  1:3045664/150‑1 (MQ=255)
                           tATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATttc                                                                                                                           >  2:1567734/1‑150 (MQ=255)
                           tATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGATCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATttc                                                                                                                           >  2:2535703/1‑150 (MQ=255)
                              tctcTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATttcttc                                                                                                                        <  2:3434520/150‑1 (MQ=255)
                                      aTGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                                <  1:2667267/150‑1 (MQ=255)
                                      aTGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                                <  1:203421/150‑1 (MQ=255)
                                      aTGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGaa                                                                                                                >  2:207459/1‑150 (MQ=255)
                                       tGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAg                                                                                                               <  2:3657828/150‑1 (MQ=255)
                                         aTTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAggg                                                                                                             >  1:4380913/1‑150 (MQ=255)
                                         aTTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAggg                                                                                                             <  1:2983103/150‑1 (MQ=255)
                                           tttGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAggggg                                                                                                           <  1:3502839/150‑1 (MQ=255)
                                             tGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGc                                                                                                         <  1:3587627/150‑1 (MQ=255)
                                             tGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGc                                                                                                         >  2:33267/1‑150 (MQ=255)
                                               ccAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAg                                                                                                       <  2:4970354/150‑1 (MQ=255)
                                                 aGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGAt                                                                                                     <  1:665832/150‑1 (MQ=255)
                                                    ccGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAAt                                                                                                  <  1:2307481/150‑1 (MQ=255)
                                                     cGTTAAGTTGTATACAAAATGCCACTATTCTAGTTGTTCTTAACTGGCTTATATTGATTAAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATg                                                                                                 <  1:2723979/150‑1 (MQ=255)
                                                       ttAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgc                                                                                               >  1:3115479/1‑150 (MQ=255)
                                                       ttAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgc                                                                                               <  2:578088/150‑1 (MQ=255)
                                                          aGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgc                                                                                            <  1:2809595/150‑1 (MQ=255)
                                                             tGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGgctgctgc                                                                                         >  1:3653453/1‑150 (MQ=255)
                                                                ataCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTcc                                                                                      >  1:470426/1‑150 (MQ=255)
                                                                ataCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTcc                                                                                      >  2:4161971/1‑150 (MQ=255)
                                                                 taCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCa                                                                                     <  2:266599/150‑1 (MQ=255)
                                                                  aCCACATGCCACTATTCGAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTAGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAg                                                                                    <  2:147621/150‑1 (MQ=255)
                                                                  aCCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAg                                                                                    <  2:3492662/150‑1 (MQ=255)
                                                                         gCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                                             <  2:3138229/150‑1 (MQ=255)
                                                                           cACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                                           >  1:1351489/1‑150 (MQ=255)
                                                                            aCTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                                          <  2:3354957/150‑1 (MQ=255)
                                                                             cTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATa                                                                         <  2:4454373/150‑1 (MQ=255)
                                                                               aTTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGACGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgc                                                                       >  1:1902393/1‑150 (MQ=255)
                                                                                  cTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctgg                                                                    >  2:2445191/1‑150 (MQ=255)
                                                                                   tAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                                   <  1:4035510/150‑1 (MQ=255)
                                                                                    aGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                                  <  1:386109/150‑1 (MQ=255)
                                                                                     gttgttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggctg                                                                 >  2:4038841/1‑150 (MQ=255)
                                                                                          tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                                            <  1:422716/150‑1 (MQ=255)
                                                                                             tAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                                         >  2:739844/1‑150 (MQ=255)
                                                                                              aaCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGc                                                        <  2:4812871/150‑1 (MQ=255)
                                                                                                 tGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAtt                                                     >  2:4345961/1‑150 (MQ=255)
                                                                                                   gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                                   >  1:919552/1‑150 (MQ=255)
                                                                                                   gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                                   >  1:3518504/1‑150 (MQ=255)
                                                                                                   gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                                   >  2:4696429/1‑150 (MQ=255)
                                                                                                   gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                               >  1:2629216/1‑138 (MQ=255)
                                                                                                    cTGATTTTGATTCAATTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACTCTGCGTGTTATTTCTTCGTGTGGAAGGGGGGCAGATACTGGCTGCTGCTCCTGTTGTGCGCATAGCTGGCTGGCGATATGCATTccc                                                  <  1:4323050/150‑1 (MQ=255)
                                                                                                    cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTccc                                                  >  1:2045560/1‑150 (MQ=255)
                                                                                                      gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                                <  2:4554851/150‑1 (MQ=255)
                                                                                                      gATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAg                                                <  2:4555104/150‑1 (MQ=255)
                                                                                                              ttCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGc                                        <  2:1129482/150‑1 (MQ=255)
                                                                                                                cAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCaa                                      <  1:4942033/150‑1 (MQ=255)
                                                                                                                 aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAAc                                     <  1:4500001/150‑1 (MQ=255)
                                                                                                                   aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTg                                   >  1:1858793/1‑150 (MQ=255)
                                                                                                                    tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                                  >  1:2501142/1‑150 (MQ=255)
                                                                                                                        gTTCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                              >  1:1479325/1‑150 (MQ=255)
                                                                                                                          tCAGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTa                            <  1:3275786/150‑1 (MQ=255)
                                                                                                                            aGGTCTTTCTTATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGc                          >  1:4531902/1‑150 (MQ=255)
                                                                                                                             ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt                         <  2:746070/150‑1 (MQ=255)
                                                                                                                                ctttcttATGCAACCTTGCTTCCAGAGCGGAAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                      <  1:4755494/150‑1 (MQ=255)
                                                                                                                                ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                      >  1:2024335/1‑150 (MQ=255)
                                                                                                                                ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                      <  1:1015514/150‑1 (MQ=255)
                                                                                                                                ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                      <  1:760785/150‑1 (MQ=255)
                                                                                                                                ctttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGAt                      <  1:683808/150‑1 (MQ=255)
                                                                                                                                 tttcttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTCATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg                     >  2:286848/1‑150 (MQ=255)
                                                                                                                                     ttATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGcgcg                 <  1:4459433/150‑1 (MQ=255)
                                                                                                                                      tATGCAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg                <  2:2020472/150‑1 (MQ=255)
                                                                                                                                       aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGc               <  1:5188916/150‑1 (MQ=255)
                                                                                                                                         gCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCAc             <  2:1924084/150‑1 (MQ=255)
                                                                                                                                          cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACg            <  2:3049230/150‑1 (MQ=255)
                                                                                                                                          cAACCATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACg            >  1:3333900/1‑150 (MQ=255)
                                                                                                                                           aaCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGc           >  1:22768/1‑150 (MQ=255)
                                                                                                                                             ccATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCtt         <  1:1969416/150‑1 (MQ=255)
                                                                                                                                             ccATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCtt         <  2:1113214/150‑1 (MQ=255)
                                                                                                                                              cATGCTTACAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCttt        >  2:2464510/1‑150 (MQ=255)
                                                                                                                                               aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTa       <  1:1490558/150‑1 (MQ=255)
                                                                                                                                                tGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAt      <  1:1615293/150‑1 (MQ=255)
                                                                                                                                                   ttCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtttt   >  2:2399740/1‑150 (MQ=255)
                                                                                                                                                   ttCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTAtttt   >  2:4115024/1‑150 (MQ=255)
                                                                                                                                                    tCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTtc  >  1:1161408/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGCCGCCGCTGCCAGCGGGTTCTTTTATTCTCTGCCATGATTTTGCCAGTCCGTTAAGTTGTATACCAAATGCCACTATTCTAGTTGTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGGAAGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGGCACGCTTTATTTTC  >  CP009273/1048779‑1049076

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: