Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,621,636 A→G 5.0% L221P (CTG→CCG)  yfgF ← cyclic‑di‑GMP phosphodiesterase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,621,6360AG 5.0% 226.9 / 2.6 80L221P (CTG→CCG) yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (23/52);  new base G (2/2);  total (26/54)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.28e-01

CCAGAGCAGCGAAATCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTA  >  CP009273/2621488‑2621783
                                                                                                                                                    |                                                                                                                                                   
ccAGAGCAGCTAAATCATCTTATAACCATAGCGCATCGCTCACCACATCAACAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCa                                                                                                                                                     <  1:165591/149‑1 (MQ=255)
      cagcGAAATCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGcc                                                                                                                                               >  1:121472/1‑149 (MQ=255)
         cGAAATCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTa                                                                                                                                                  >  2:897462/1‑143 (MQ=255)
           aaaTCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAgc                                                                                                                                          >  1:1113278/1‑149 (MQ=255)
           aaaTCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAgc                                                                                                                                          >  1:472177/1‑149 (MQ=255)
             aTCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCa                                                                                                                                              <  2:1110728/143‑1 (MQ=255)
             aTCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGCGGCCGGAAGAAGGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAgcgc                                                                                                                                        <  1:1269789/149‑1 (MQ=255)
             aTCAGCTTATAACCATAGCGCAGCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAgcgc                                                                                                                                        <  1:492032/149‑1 (MQ=255)
               cAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGa                                                                                                                                      <  1:901490/149‑1 (MQ=255)
                   ttATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAActc                                                                                                                                  >  1:404693/1‑149 (MQ=255)
                   ttATAACCATAGCGCATCGCTCCCAACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGAGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATACAAGTAGCAAGAGAGCGAAcgc                                                                                                                                  >  1:966636/1‑147 (MQ=255)
                       aaCCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCtttt                                                                                                                              <  1:458908/149‑1 (MQ=255)
                         ccATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTGTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCtttttt                                                                                                                            <  1:74714/149‑1 (MQ=255)
                          cctAGCGCATCGCTCCCCACATCATCAGGGGCAGCGATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTg                                                                                                                           <  1:724323/147‑1 (MQ=255)
                          cATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTg                                                                                                                           <  1:1165852/149‑1 (MQ=255)
                                gcATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCCGTAGCCAGAGCGCGAACTCTTTTTTGGTGACt                                                                                                                     >  1:600922/1‑149 (MQ=255)
                                      cTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCa                                                                                                               <  1:264992/149‑1 (MQ=255)
                                             catcatCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCGGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACct                                                                                                        >  1:56130/1‑149 (MQ=255)
                                             catcatCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCGGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACct                                                                                                        >  1:56043/1‑149 (MQ=255)
                                              atcatcAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACctg                                                                                                       >  1:1031799/1‑149 (MQ=255)
                                               tcatcaGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAAATCTTTTTTGGTGACTTTGGCATCAACctgc                                                                                                      >  1:212739/1‑149 (MQ=255)
                                                    aGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACctg                                                                                                       <  2:632096/143‑1 (MQ=255)
                                                    aGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACctg                                                                                                       <  2:767960/143‑1 (MQ=255)
                                                         cagcaATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGtt                                                                                                  >  2:657768/1‑143 (MQ=255)
                                                          agcaATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTg                                                                                           <  1:355194/149‑1 (MQ=255)
                                                          agcaATAAGGACAAGGTATAATTGGGGTTAAAAATTTTGCTTTTTTCACTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTg                                                                                           <  1:914018/149‑1 (MQ=255)
                                                           caaataatGACAAGGTATAATTGGTGCAAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGc                                                                                          <  1:22630/147‑1 (MQ=255)
                                                               taatGAGAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAAtt                                                                                            <  2:875180/143‑1 (MQ=255)
                                                                 atGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGaataa                                                                                    <  1:798125/149‑1 (MQ=255)
                                                                      aatgTATAATTGGTGCTAAAAATTGTGCTTTTGCCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCt                                                                               <  1:358285/146‑1 (MQ=255)
                                                                         gtataatagaaaaaaaaaaTTGTGCTTTTTACATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAataata                                                                                  <  2:191507/129‑1 (MQ=255)
                                                                         gTATAATTGGTGCTAAAAACTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAataata                                                                                  <  2:925944/143‑1 (MQ=255)
                                                                          tataATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGt                                                                           <  1:921300/149‑1 (MQ=255)
                                                                                     cTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATgg                                                                >  1:356059/1‑149 (MQ=255)
                                                                                      tAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACGTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGa                                                               <  1:254121/149‑1 (MQ=255)
                                                                                        aaaaTTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCGGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATaaaa                                                                   <  2:47318/143‑1 (MQ=255)
                                                                                        aaaaTTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCGGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATaaaa                                                                   <  2:47234/143‑1 (MQ=255)
                                                                                        aaaaTTGTGCTTTTTTCATTTAACCGCCTGCATCACAGCAACAATAAAGCACCCAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATaaaa                                                                   <  2:1086054/143‑1 (MQ=255)
                                                                                         aaaTTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGAttt                                                            <  1:443385/149‑1 (MQ=255)
                                                                                           aTTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCg                                                          <  1:134383/149‑1 (MQ=255)
                                                                                               tacctttttCACTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACt                                                      <  1:486556/145‑1 (MQ=255)
                                                                                                 cTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCg                                                          <  2:523543/143‑1 (MQ=255)
                                                                                                 cTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGAGCGAACTATTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGGTGCGTAAATAAAATGGAGTTCGCACTag                                                    >  1:268972/1‑148 (MQ=255)
                                                                                                  ttttttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGGACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTGCGTAAATAAAATGGATTTCGc                                                         >  2:253815/1‑143 (MQ=255)
                                                                                                  ttttttCATTTAACGGCATGCATAACAGAAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCCTCAACCTGCTGTTTTAAGTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTAcc                                                   >  1:1161934/1‑149 (MQ=255)
                                                                                                     tttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATGATCCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACt                                                      <  2:743262/143‑1 (MQ=255)
                                                                                                     tttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTGAAAAAAATGGATTGCGCACt                                                      >  2:587291/1‑143 (MQ=255)
                                                                                                      ttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGTCTTTGGCATCAACCTGCTGTTTTATTTGCGATTAATAGCTACGTATATAAAATGGATTTCGCACTa                                                     >  2:447033/1‑143 (MQ=255)
                                                                                                      ttCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGaa                                               >  1:989540/1‑149 (MQ=255)
                                                                                                        cATTTAACGGCATGCATGACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATCAAATGGATTTCGCACTACCCGAatg                                             <  1:654255/149‑1 (MQ=255)
                                                                                                        cATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAatg                                             <  1:545399/149‑1 (MQ=255)
                                                                                                          tttAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAatgat                                           <  1:513282/149‑1 (MQ=255)
                                                                                                            tAACGGCAGGCACAACAGTAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGaa                                               <  2:794891/143‑1 (MQ=255)
                                                                                                             aaCGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGaa                                        <  1:621153/149‑1 (MQ=255)
                                                                                                               cGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCAAAACGCCTTTATTGTATTTCGCACTACCCGAATGATGAAGt                                      <  1:958998/149‑1 (MQ=255)
                                                                                                                  gatggatAACAGCAATAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATCGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAAtgatga                                         <  2:54987/142‑1 (MQ=255)
                                                                                                                    tgcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGAGTTCGCACTACCCGAATGATGAAGTAGcgc                                 >  1:574712/1‑147 (MQ=255)
                                                                                                                     gcatAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGt                                      <  2:765089/143‑1 (MQ=255)
                                                                                                                       atAACAGCAACAATAAAGCACCTAAAGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGc                              >  1:655053/1‑149 (MQ=255)
                                                                                                                        tAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCg                             <  1:244867/149‑1 (MQ=255)
                                                                                                                             gCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGc                              <  2:136893/143‑1 (MQ=255)
                                                                                                                              caacaaTAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTCTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGacac                       <  1:1297406/149‑1 (MQ=255)
                                                                                                                                 caagtaaGCAACTAATGCCAGTGGCCAGAGCGCGACCTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTTCGAATAATAGCTACGTAAATCAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGa                    <  1:259150/144‑1 (MQ=255)
                                                                                                                                 caaTAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGa                          <  2:435229/143‑1 (MQ=255)
                                                                                                                                 caaTAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGa                          <  2:1001004/143‑1 (MQ=255)
                                                                                                                                  aaTAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGCTGACTTTGGCATCAACCTGCTGTTTTAATTGCGCATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAACGATGAAGTAGCACAGCGGGACACCGAt                   <  1:695570/149‑1 (MQ=255)
                                                                                                                                   aTAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGAAGCACAGCGGGaca                        >  2:261838/1‑143 (MQ=255)
                                                                                                                                   aTAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGACATCAACCTGCAGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCAGGacc                        >  2:89165/1‑142 (MQ=255)
                                                                                                                                    tAAAGCACTTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCGACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGacac                       <  2:137934/143‑1 (MQ=255)
                                                                                                                                     aaaGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAg                <  1:1256265/149‑1 (MQ=255)
                                                                                                                                      tggCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTATTTGCGAATAATAGGTACGTGAATAAAAGGGATTTGGCAGTAGCGGGATGTTGGAGTGGCACGGGGGGAGACCg                     >  2:522824/3‑143 (MQ=255)
                                                                                                                                         cACCTAATGCCAGTAGCCAGAGCGTGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATTATGAAGTAGCACAGCGGGACACCGATCAGATTa            <  1:1120035/149‑1 (MQ=255)
                                                                                                                                             tAATGCCAGTAGCCAGAGGGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACcca        <  1:1109230/149‑1 (MQ=255)
                                                                                                                                             tAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGAt              <  2:30500/143‑1 (MQ=255)
                                                                                                                                             gaaTGCCAGTAGCCAGGGCGCGAGCTCTTTTTTGGTGACTGTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGAt              <  2:456781/142‑1 (MQ=255)
                                                                                                                                              aaTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACccac       <  1:762079/149‑1 (MQ=255)
                                                                                                                                               aTGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACccacc      <  1:494322/149‑1 (MQ=255)
                                                                                                                                               aTGCCAGTAACCCGATCGCGAACTCGTTTTTGTGGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTa            >  2:1140023/1‑143 (MQ=255)
                                                                                                                                                tGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACccacca     <  1:715742/149‑1 (MQ=255)
                                                                                                                                                  ccAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTAccc         <  2:379875/143‑1 (MQ=255)
                                                                                                                                                   cAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGAGGAAGTAGCACAGCGGGACACCGATCCGATTACCCACCAGTa  >  1:1198390/1‑149 (MQ=255)
                                                                                                                                                    agtagtCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACccac       <  2:1050721/143‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGAGCAGCGAAATCAGCTTATAACCATAGCGCATCGCTCCCCACATCATCAGGGGCAGCAATAATGACAAGGTATAATTGGTGCTAAAAATTGTGCTTTTTTCATTTAACGGCATGCATAACAGCAACAATAAAGCACCTAATGCCAGTAGCCAGAGCGCGAACTCTTTTTTGGTGACTTTGGCATCAACCTGCTGTTTTAATTGCGAATAATAGCTACGTAAATAAAATGGATTTCGCACTACCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTA  >  CP009273/2621488‑2621783

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