Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,909,758 C→A 81.4% E12* (GAA→TAA)  proQ ← RNA chaperone, putative ProP translation regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,909,7580CA81.4% 177.6 / 41.2 70E12* (GAA→TAA) proQRNA chaperone, putative ProP translation regulator
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base A (32/25);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

AGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCG  >  CP009273/1909612‑1909904
                                                                                                                                                   |                                                                                                                                                  
aGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTAttt                                                                                                                                                 <  2:247610/150‑1 (MQ=255)
   ggATCGCAATTGCGTTTTGCTCAGGTTCATTTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTAc                                                                                                                                               >  1:343348/1‑150 (MQ=255)
   ggATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACt                                                                                                                                              >  2:696131/1‑150 (MQ=255)
       cGCAATTGCGTTTTGCTCAGGTTCATTTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCt                                                                                                                                           <  2:322164/150‑1 (MQ=255)
         cAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTAtt                                                                                                                                        <  1:291291/150‑1 (MQ=255)
                gTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACt                                                                                                                                   >  2:208518/1‑148 (MQ=255)
                gTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACt                                                                                                                                   >  1:268156/1‑148 (MQ=255)
                 ttttGTTCAGGTTCATTCCCCAAGCAACACGATCGACCAAATCCTGAAAAAAACCAATTTTCAGCGGACGCGCTCCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTATTACCTGTAATCCAACTAAg                                                                                                                                >  2:837119/1‑150 (MQ=255)
                      cTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTg                                                                                                                            <  2:322467/149‑1 (MQ=255)
                       tCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGAt                                                                                                                          >  1:367013/1‑150 (MQ=255)
                            ttCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTcc                                                                                                                     <  1:65545/150‑1 (MQ=255)
                                     cccAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAtt                                                                                                             >  1:283308/1‑150 (MQ=255)
                                     cccAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAtt                                                                                                             <  2:582479/150‑1 (MQ=255)
                                      ccAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAttt                                                                                                            >  2:464629/1‑150 (MQ=255)
                                      ccAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAttt                                                                                                            >  1:679393/1‑150 (MQ=255)
                                          cAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTg                                                                                                        <  2:620498/150‑1 (MQ=255)
                                                gATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTAc                                                                                                   >  2:806013/1‑149 (MQ=255)
                                                    gACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACg                                                                                               >  1:295089/1‑149 (MQ=255)
                                                    gACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACg                                                                                               <  2:640743/149‑1 (MQ=255)
                                                      ccAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGAc                                                                                            >  1:186517/1‑150 (MQ=255)
                                                      ccAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGAc                                                                                            >  1:668553/1‑150 (MQ=255)
                                                            ccTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGc                                                                                       >  2:19251/1‑149 (MQ=255)
                                                                aaaaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaaca                                                                                  <  1:19251/150‑1 (MQ=255)
                                                                aaaaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaac                                                                                   >  1:877731/1‑149 (MQ=255)
                                                                 aaaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaac                                                                                   >  1:885206/1‑148 (MQ=255)
                                                                  aaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAg                                                                                >  2:86102/1‑150 (MQ=255)
                                                                      aCCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAgcg                                                                              >  2:281914/1‑148 (MQ=255)
                                                                       ccGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAg                                                                           <  2:706723/150‑1 (MQ=255)
                                                                       ccGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAg                                                                           <  2:658949/150‑1 (MQ=255)
                                                                               cAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGaa                                                                    >  2:777177/1‑149 (MQ=255)
                                                                                aGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACa                                                                  <  2:343348/150‑1 (MQ=255)
                                                                                  cGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAgg                                                                >  1:500914/1‑150 (MQ=255)
                                                                                    gACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCg                                                              <  1:134201/150‑1 (MQ=255)
                                                                                        cgcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCaa                                                          <  1:389534/150‑2 (MQ=255)
                                                                                         gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCAt                                                          >  2:876776/1‑149 (MQ=255)
                                                                                         gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCAt                                                          >  1:4378/1‑149 (MQ=255)
                                                                                         gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCAt                                                          >  1:897290/1‑149 (MQ=255)
                                                                                           ttCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATta                                                        >  2:425584/1‑149 (MQ=255)
                                                                                           ttCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATta                                                        >  1:608737/1‑149 (MQ=255)
                                                                                             cACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATtat                                                       <  2:616590/148‑1 (MQ=255)
                                                                                               ccTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAAt                                                    >  2:243988/1‑149 (MQ=255)
                                                                                                    cGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCt                                              <  1:1046400/150‑1 (MQ=255)
                                                                                                       aCTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATc                                           <  1:850202/150‑1 (MQ=255)
                                                                                                       aCTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTAt                                            >  1:1017126/1‑149 (MQ=255)
                                                                                                       aCTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATc                                           <  2:994351/150‑1 (MQ=255)
                                                                                                              cAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAAt                                    <  2:268156/150‑1 (MQ=255)
                                                                                                              cAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAAt                                    >  1:130863/1‑150 (MQ=255)
                                                                                                                     gAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGt                             <  2:808204/150‑1 (MQ=255)
                                                                                                                      aaaaCGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGtt                            <  1:330054/150‑1 (MQ=255)
                                                                                                                        aaCGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTAt                          <  1:464629/150‑1 (MQ=255)
                                                                                                                        aaCGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTAt                          >  2:561031/1‑150 (MQ=255)
                                                                                                                          cGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATcc                        <  2:414787/150‑1 (MQ=255)
                                                                                                                             tCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGtt                     >  1:331666/1‑150 (MQ=255)
                                                                                                                              cGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGtt                     >  1:943638/1‑149 (MQ=255)
                                                                                                                               ggCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGa                   >  2:561416/1‑150 (MQ=255)
                                                                                                                               ggCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGa                   <  1:1020065/150‑1 (MQ=255)
                                                                                                                                 ccAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGAtt                 <  1:549152/150‑1 (MQ=255)
                                                                                                                                  cAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTa                <  2:367013/150‑1 (MQ=255)
                                                                                                                                    gAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATc              >  2:657196/1‑150 (MQ=255)
                                                                                                                                    gAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATc              <  2:668553/150‑1 (MQ=255)
                                                                                                                                    gAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATc              >  1:811223/1‑150 (MQ=255)
                                                                                                                                      aaCGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTa                           >  1:507324/1‑135 (MQ=255)
                                                                                                                                      aaCGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCt            <  2:500914/150‑1 (MQ=255)
                                                                                                                                       aCGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTAt                          <  2:507324/135‑1 (MQ=255)
                                                                                                                                       aCGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTg           >  2:240009/1‑150 (MQ=255)
                                                                                                                                        cgcgATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGc          <  1:574549/150‑1 (MQ=255)
                                                                                                                                           gATTACTTATTTACTGCTTTTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGAc       >  1:148022/1‑150 (MQ=255)
                                                                                                                                           gATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGAc       >  2:173625/1‑150 (MQ=255)
                                                                                                                                            aTTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACg      >  2:577495/1‑150 (MQ=255)
                                                                                                                                               aCTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCt    <  1:281914/149‑1 (MQ=255)
                                                                                                                                                cttattTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTc   <  1:837119/149‑1 (MQ=255)
                                                                                                                                                cttattTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCg  <  2:1058980/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCG  >  CP009273/1909612‑1909904

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: