Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,602 T→C 9.9% pseudogene (954/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,460,6020TC 9.9% 175.7 / 2.4 74pseudogene (954/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base C (4/4);  total (41/33)
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.86e-01

GACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATA  >  CP009273/1460467‑1460744
                                                                                                                                       |                                                                                                                                              
gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACGAGGCCATCGTCAATGACGAAGGCGAAGGCACCATCACCAACGGTGGCACCGGGACCCAGATTAATGgtggtg                                                                                                                                         >  2:816045/1‑143 (MQ=17)
gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCa                                                                                                                                   >  1:514239/1‑149 (MQ=255)
 aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgatga                                                                                                                                        <  2:419157/143‑1 (MQ=255)
    gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGc                                                                                                                               >  1:708497/1‑149 (MQ=255)
      ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                                                                                                             >  1:1038476/1‑149 (MQ=255)
       ggacaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGCGCACCATCACCAACGGGGGCACCGGCACGCAGATTAACGGTGATGACGCCAc                                                                                                                                  >  2:251584/5‑143 (MQ=17)
        aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGTAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                           <  1:787431/149‑1 (MQ=255)
          aaCAAGGGAACCATGACCGCCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaa                                                                                                                         <  1:569025/149‑1 (MQ=255)
           aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATGGACGGTGACAAGGCCATCGTCAATAACGAAGGCGCGAGCACCATCACCAACGGTGGCACCGGACACCAGATTACACGTGCTGACGCCACCGCa                                                                                                                              >  2:31676/1‑143 (MQ=17)
           aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGGGATGACGCCACGGCa                                                                                                                              >  2:716455/1‑143 (MQ=255)
            cAAGGGCACCATGACCGTCACCGATCCGGAGTCCCTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACCCAGATTAATGGTGATGACGCCCCGGCAaacaaca                                                                                                                       >  1:65077/1‑149 (MQ=255)
            cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATAGGTATCCAGATCGACCTTGACAAGGCCATCGACAATAACGAAGGCGAAAGCACCATCACCAACGAAGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                                                                                                       >  1:484001/1‑149 (MQ=18)
               gggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACg                                                                                                                    <  1:311095/149‑1 (MQ=255)
                   aCCATGACCGTCACCGATCCGTAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTCATAGGAATGACTCCACGGCAaacaataa                                                                                                                      <  2:521015/143‑1 (MQ=18)
                   aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                      >  2:1085970/1‑143 (MQ=255)
                    ccATGACCGTCACCGACCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAACACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                     <  2:186585/143‑1 (MQ=255)
                     cATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                              <  1:1185933/149‑1 (MQ=255)
                        gACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACACCGGCAAAacca                                                                                                           >  1:162973/1‑149 (MQ=255)
                             tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacc                                                                                                            >  2:1055/1‑143 (MQ=255)
                               aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCCACAACGACGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCGCGCGTCGTGCTCATGACGCCACGGCAAGCAACAACGGCAAAacgac                                                                                                          >  2:118694/1‑143 (MQ=18)
                               aCCGATCCGGAGTCCATGGGTATCAAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCCGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAaccac                                                                                                          <  2:557919/143‑1 (MQ=255)
                                ggatctcGGAGCCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                         <  2:694159/137‑1 (MQ=255)
                                ccatctcGGAGTCCATCGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCATACAACGACGCCAAAaccacc                                                                                                         <  2:267187/137‑1 (MQ=11)
                                ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGACATCGACAAAAACGAAGGCGAGAGCACCATCACCAACCGTGGCACCGGCACCCAGATTAATGGTGATGACGCAACGGACCACACCAAAGACAACACCACGGATGat                                                                                                   >  1:1094532/1‑148 (MQ=255)
                                 cGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGGTGACg                                                                                                  >  1:93909/1‑149 (MQ=255)
                                  gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                 <  1:417601/149‑1 (MQ=255)
                                  gATCCGGAGTCCACGGGTATCCAGATCGACGGTGACAAGGCCATCGTCCATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                 <  1:599687/149‑1 (MQ=255)
                                   aTCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                                                <  1:891411/149‑1 (MQ=255)
                                    tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                     <  2:384049/143‑1 (MQ=255)
                                    tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccaagctgg                                                                                                     >  2:176749/1‑137 (MQ=255)
                                      cGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACCACGGCAAAACCACGGTTGAAGGCAAg                                                                                             >  1:1098304/1‑149 (MQ=255)
                                       ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTGGACGGCAggg                                                                                            >  1:306229/1‑149 (MQ=255)
                                       ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGCGAACGGGAAgg                                                                                            >  1:802725/1‑149 (MQ=18)
                                          gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                         <  1:730276/149‑1 (MQ=255)
                                           ctcaggggTATCCAGATCGACGGTGACAAGGCCATCGTCAACAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                              <  2:706700/138‑1 (MQ=255)
                                              aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAACGTTGACGGGCAGGGTTCGAc                                                                                     >  1:816121/1‑149 (MQ=21)
                                                  gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGAGGCCACGGCAAACAACCAAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                                                 >  1:498122/1‑149 (MQ=255)
                                                  gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTcc                                                                                       >  2:1086920/1‑143 (MQ=255)
                                                    aTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                     <  2:1011842/143‑1 (MQ=255)
                                                      ccAGATCGACGGTGACAAGGCTATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                             >  1:921740/1‑149 (MQ=255)
                                                      ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                             >  1:521874/1‑149 (MQ=255)
                                                        aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                           >  1:1103751/1‑149 (MQ=255)
                                                         ggtCGACGGCGACCAGGCCATCGTCAACAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGGGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                <  2:428195/141‑1 (MQ=11)
                                                            cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAACAGAACAGTTGCCGGCAGGGaaaccgcggtggcgg                                                                             >  2:877757/1‑129 (MQ=18)
                                                             gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGCAATCaa                                                                      >  1:1242952/1‑149 (MQ=255)
                                                              aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                     >  1:499559/1‑149 (MQ=255)
                                                               cGGTGACAAGGCCATCGTCAATAGCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGATGACGGCAAGGATTCCACCGGTCCGGGAATCAACg                                                                    >  1:725170/1‑149 (MQ=255)
                                                                      aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCCAGGATTCCACCGGTACGGAAATCAACGGtaataa                                                             >  1:249650/1‑149 (MQ=255)
                                                                      aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaataa                                                             >  1:1146910/1‑149 (MQ=255)
                                                                      aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaataa                                                             >  1:273200/1‑149 (MQ=255)
                                                                           cATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACAAAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCAACGTTGACGGCAAGGATTCCAACGGTACGGAAATCCACGCtaata                                                              >  2:406648/1‑143 (MQ=14)
                                                                             tCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAACCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                      <  1:1026302/149‑1 (MQ=255)
                                                                              cGTCAATAACGAAGGCAAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAGTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACg                                                           >  2:800341/1‑143 (MQ=255)
                                                                                   aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                                <  1:1584/149‑1 (MQ=255)
                                                                                             cGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCg                                      <  1:1229300/149‑1 (MQ=255)
                                                                                                gagCACCATCACCAACGGTGGCACCGGCACCCAGATGAATGGGGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAATCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATc                                   <  1:917993/149‑1 (MQ=255)
                                                                                                 agCACTAACGCCAACGGTGGCACCGGCACCCAGATCAAGGGTGAGGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt                                  <  1:879871/149‑1 (MQ=18)
                                                                                                    accaccAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                               <  1:416702/149‑1 (MQ=255)
                                                                                                      cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                             >  1:784179/1‑149 (MQ=255)
                                                                                                          aCCAACGGTTGCACCGGCACCCATATTAATGGTGATGTCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                               <  2:157024/143‑1 (MQ=255)
                                                                                                          aCCAACGGTGGCACCGGCACCCAGATTACTGGTGATGACGCCACGGCACGCCACAACGCCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                         >  1:1224686/1‑149 (MQ=18)
                                                                                                          aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                         >  1:580491/1‑149 (MQ=255)
                                                                                                             aaCGGTGGCACCGGCACCCAGATTAATGGTGATGATGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGt                            >  2:397887/1‑143 (MQ=255)
                                                                                                             aaCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcgg                      <  1:826505/149‑1 (MQ=255)
                                                                                                                     caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGTAGTAACGGAAGAGGTATTCAGGACGGCGGTCGGGATGTcggggggg                    >  2:854295/1‑136 (MQ=17)
                                                                                                                     caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcg                    <  2:448516/143‑1 (MQ=255)
                                                                                                                       ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcggc                  <  2:1010783/143‑1 (MQ=255)
                                                                                                                       ccgccaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcggc                  <  2:218045/143‑1 (MQ=255)
                                                                                                                       cccgcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcggc                  <  2:238307/140‑1 (MQ=255)
                                                                                                                        cggcaccCAGATTAATGGTGATGAGGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGAAGTCAGCGGCGGCGGTCACg           >  1:965056/1‑149 (MQ=255)
                                                                                                                          gcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAGGTGATTCAGGACGGCGATCCGGATGTCAGcggcggcggg               >  2:665889/1‑142 (MQ=255)
                                                                                                                          gaaccCAGTGTAATGGGGACGAAGCCACGGCAAACCACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt         <  1:396094/147‑1 (MQ=17)
                                                                                                                             cccAGATTAATGGTTATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGACTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt      <  1:119088/149‑1 (MQ=255)
                                                                                                                               cAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa    <  1:137265/149‑1 (MQ=255)
                                                                                                                                       tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata  >  2:1085979/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                              
GACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATA  >  CP009273/1460467‑1460744

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