Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,422,549 G→A 5.2% V222M (GTG→ATG)  PP_3038 → TraC domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,422,5490GA 5.2% 233.6 / 2.3 76V222M (GTG→ATG) PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (40/32);  new base A (2/2);  total (42/34)
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.31e-01

GGGGCGGGACAAGTCCTACTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTGGAATCCTGGTGA  >  NC_002947/3422401‑3422691
                                                                                                                                                    |                                                                                                                                              
ggggCGGGACAAGTCCTACTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGt                                                                                                                                               <  1:736391/150‑1 (MQ=255)
 gggCGGGACAAGTCCTACTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTg                                                                                                                                              <  2:260039/150‑1 (MQ=255)
       gACAAGTCCTACTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGGGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATc                                                                                                                                        <  2:606278/150‑1 (MQ=255)
            gTCCTACTGGCCGCCTCGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCCCCCAGGCGGTGGCCATCGCCtt                                                                                                                                   <  2:40303/150‑1 (MQ=255)
                tACTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGAt                                                                                                                               <  1:369405/150‑1 (MQ=255)
                tACTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGAt                                                                                                                               <  2:371076/150‑1 (MQ=255)
                  cTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGc                                                                                                                             >  1:632817/1‑150 (MQ=255)
                   tGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCtt                                                                                                                                   >  2:275151/1‑143 (MQ=255)
                   tGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCtt                                                                                                                                   <  1:275151/143‑1 (MQ=255)
                     gCCGCATGGCATGGCAAAGGACGGTGATTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGTCTGTGCTGGGGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCggg                                                                                                                          >  2:431261/1‑150 (MQ=255)
                                   ccAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCg                                                                                                                            >  1:241735/1‑134 (MQ=255)
                                   ccAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCg                                                                                                                            <  2:241735/134‑1 (MQ=255)
                                          tGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCaa                                                                                                     <  2:632817/150‑1 (MQ=255)
                                               cTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCata                                                                                                >  1:1411043/1‑150 (MQ=255)
                                                tttCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCatat                                                                                               >  2:1192930/1‑150 (MQ=255)
                                                  tCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGc                                                                                             <  2:236359/150‑1 (MQ=255)
                                                      caTGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGTTCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGa                                                                                         >  1:784975/1‑150 (MQ=255)
                                                         gATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACg                                                                                      >  2:948897/1‑150 (MQ=255)
                                                            tggtggCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCtt                                                                                   >  1:393009/1‑150 (MQ=255)
                                                                 gCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTg                                                                              >  2:786516/1‑150 (MQ=255)
                                                                   cATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTggt                                                                            >  1:747756/1‑150 (MQ=255)
                                                                       cTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTggtcggt                                                                        >  2:1511561/1‑150 (MQ=255)
                                                                        tGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTggtcggtc                                                                       <  1:1511561/150‑1 (MQ=255)
                                                                        tGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTggtcggtc                                                                       >  1:1355752/1‑150 (MQ=255)
                                                                           gCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCCCCGGGGCTAGCCCGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATACGTGAACGCTTCCATGGTCGGTcccc                                                                    >  2:140814/1‑148 (MQ=255)
                                                                             caggTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatca                                                                  >  2:730330/3‑150 (MQ=255)
                                                                                gTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatcatca                                                               <  2:747756/150‑1 (MQ=255)
                                                                                    gCCTGTGCTGGTGTGTGAGCGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATg                                                           <  2:784975/150‑1 (MQ=255)
                                                                                     ccTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCgg                                                                                                                           >  2:118863/1‑85 (MQ=255)
                                                                                     ccTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCgg                                                                                                                           <  1:118863/85‑1 (MQ=255)
                                                                                     ccTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAAc                                                                                       <  1:670502/121‑1 (MQ=255)
                                                                                     ccTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAAc                                                                                       >  2:670502/1‑121 (MQ=255)
                                                                                     ccTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGc                                                          >  2:731539/1‑150 (MQ=255)
                                                                                       tgtgCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCCCATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCca                                                                     <  1:1252543/137‑1 (MQ=255)
                                                                                       tgtgCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCtt                                                                                   >  2:1464951/1‑123 (MQ=255)
                                                                                       tgtgCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCtt                                                                                   <  1:1464951/123‑1 (MQ=255)
                                                                                       tgtgCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCca                                                                     >  2:1252543/1‑137 (MQ=255)
                                                                                       tgtgCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCcg                                                        >  2:709030/1‑150 (MQ=255)
                                                                                       tgtgCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCcg                                                        >  2:521817/1‑150 (MQ=255)
                                                                                            tGGTGTGTGAGGGCTACGCCACTGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGTTCGGTCCATCATCATATGCCGCGATg                                                   >  2:1440546/1‑150 (MQ=255)
                                                                                              gtgtgtGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGc                                                                                                                             >  1:1116105/1‑74 (MQ=255)
                                                                                              gtgtgtGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGc                                                                                                                             <  2:1116105/74‑1 (MQ=255)
                                                                                               tgtgtgAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATgacg                                                >  1:213643/1‑150 (MQ=255)
                                                                                                  gtgAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTccc                                                                                <  2:316518/115‑1 (MQ=255)
                                                                                                  gtgAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTccc                                                                                >  1:316518/1‑115 (MQ=255)
                                                                                                          tACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGAcc                                     >  1:1537205/1‑150 (MQ=255)
                                                                                                              ccACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAAc                                                                                                                      >  2:586297/1‑65 (MQ=255)
                                                                                                              ccACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAAc                                                                                                                      <  1:586297/65‑1 (MQ=255)
                                                                                                               cACCGGGGCTAGCCTGCACATGGCGACCTGCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCAACATATGCCGCGATGACGACTGGAAGACCAAGCg                                >  2:416480/1‑150 (MQ=255)
                                                                                                               cACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCg                                >  2:1159349/1‑150 (MQ=255)
                                                                                                                aCCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCt                                                                                                                                    <  2:576712/49‑1 (MQ=255)
                                                                                                                aCCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCt                                                                                                                                    >  1:576712/1‑49 (MQ=255)
                                                                                                                aCCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCgg                               >  1:156547/1‑150 (MQ=255)
                                                                                                                 ccGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGc                              <  1:521817/150‑1 (MQ=255)
                                                                                                                        tAGGCTGCACATTGGGACCGCCCAGGCGGTGGCCATGGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGGTTCCCTGGTCGGTCCATCATCATATGCCGCGGTGACGAGTGGAAGACCAAGCGGCCGACggg                       >  2:774629/1‑150 (MQ=255)
                                                                                                                        tAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACggg                       >  2:675170/1‑150 (MQ=255)
                                                                                                                                acaTGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGAGGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCt               <  2:393009/150‑1 (MQ=255)
                                                                                                                                   tGGCGGCCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCCGGCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTGGa            <  1:140814/150‑1 (MQ=255)
                                                                                                                                    ggCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAATACCAAGCGGCCGACGGGCGAACCCTGGaa           <  1:1192930/150‑1 (MQ=255)
                                                                                                                                          cTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTg              <  2:1031645/141‑1 (MQ=255)
                                                                                                                                          cTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTg              >  1:1031645/1‑141 (MQ=255)
                                                                                                                                           taccAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGGCGCCAAGCATATGCGTGAACGCTTCCCTGGTTGGTCCATCATCATATGCCGCGAGGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTGGAATCCTGGt    <  1:1383805/148‑1 (MQ=255)
                                                                                                                                           tCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatcatc                                                                >  2:1133995/1‑90 (MQ=255)
                                                                                                                                           tCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatcatc                                                                >  2:1133306/1‑90 (MQ=255)
                                                                                                                                           tCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatcatc                                                                <  1:1133995/90‑1 (MQ=255)
                                                                                                                                           tCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatcatc                                                                <  1:1133306/90‑1 (MQ=255)
                                                                                                                                             ccAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatcatca                                                               >  1:956483/1‑89 (MQ=255)
                                                                                                                                             ccAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCcatcatca                                                               <  2:956483/89‑1 (MQ=255)
                                                                                                                                             ccAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCTACGGGCGAACCCTGGAATCCTGGTGa  >  1:821212/1‑150 (MQ=255)
                                                                                                                                             ccAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTGGAATCCTGGTGa  <  1:1440546/150‑1 (MQ=255)
                                                                                                                                              aaGGCGATGGCCACCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGc                                 >  2:1535076/2‑118 (MQ=255)
                                                                                                                                              aaGGCGATGGCCACCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGc                                 >  2:1534500/2‑118 (MQ=255)
                                                                                                                                              aaGGCGATGGCCACCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGc                                 <  1:1535076/117‑1 (MQ=255)
                                                                                                                                              aaGGCGATGGCCACCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGc                                 <  1:1534500/117‑1 (MQ=255)
                                                                                                                                               aGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTGGAATCCTgg     >  2:668047/1‑145 (MQ=255)
                                                                                                                                               aGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTGGAATCCTgg     <  1:668047/145‑1 (MQ=255)
                                                                                                                                                    gTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGcc                             >  2:462279/1‑116 (MQ=255)
                                                                                                                                                    gTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGcc                             <  1:462279/116‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGGGCGGGACAAGTCCTACTGGCCGCATGGCATGGCCAAGGATGGTGCTTTCCACATGATTGGTGGCCATCCTGAGCCCGGTGAGCCTGTGCTGGTGTGTGAGGGCTACGCCACCGGGGCTAGCCTGCACATGGCGACCTCCCAGGCGGTGGCCATCGCCTTCGATGCGGGCAACCTGATGGCGGTCGCCAAGCATATGCGTGAACGCTTCCCTGGTCGGTCCATCATCATATGCCGCGATGACGACTGGAAGACCAAGCGGCCGACGGGCGAACCCTGGAATCCTGGTGA  >  NC_002947/3422401‑3422691

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