Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,279,021 G→A 100% A247V (GCG→GTG)  bcr ← bicyclomycin/cysteine/sulfonamide efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,279,0210GA98.6% 307.9 / ‑4.3 93A247V (GCG→GTG) bcrbicyclomycin/cysteine/sulfonamide efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (44/47);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCGCCAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCGCGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGA  >  NC_000913/2278925‑2279117
                                                                                                |                                                                                                
aaCGCGCCAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGa                                                                                               >  2:17519/1‑100 (MQ=255)
  cgcgCCAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCacgacg                                                                                             <  1:490579/100‑1 (MQ=255)
  cgcgCCAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCacgacg                                                                                             >  1:1596279/1‑100 (MQ=255)
      ccAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAGGTAATAACCAAAGTTTTCCGGCACGACGTGGt                                                                                         <  2:1759573/100‑1 (MQ=255)
      ccAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGt                                                                                         <  2:1180526/100‑1 (MQ=255)
       cAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGtt                                                                                        <  2:260763/100‑1 (MQ=255)
          tgctgcGGACGATGCGGCTGTGAAAGATCGTCATCACTATCAGGAGAACAGTGGTGAACTCGGTGTAATTGCCATAGTTTTCCGTCGCGACGTGGTTAAt                                                                                     <  2:598117/100‑1 (MQ=255)
           gcggcgGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAAtt                                                                                    >  1:949126/1‑100 (MQ=255)
           gcggcgGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAAtt                                                                                    <  2:666505/100‑1 (MQ=255)
            cggcggACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAAttt                                                                                   >  2:1656794/1‑100 (MQ=255)
              gcggACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCa                                                                                 >  1:1232917/1‑100 (MQ=255)
               cggACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCaa                                                                                <  1:1273680/100‑1 (MQ=255)
                ggACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAat                                                                               <  2:552207/100‑1 (MQ=255)
                  aCGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAatat                                                                             <  1:1487138/100‑1 (MQ=255)
                   cGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAAtata                                                                            <  1:2476665/100‑1 (MQ=255)
                     aaGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATaaa                                                                          >  2:420550/1‑100 (MQ=255)
                     aaGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATaaa                                                                          <  2:1621730/100‑1 (MQ=255)
                     aaGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATaaa                                                                          <  1:2063368/100‑1 (MQ=255)
                      aGCGGCTGTTAAGGATGGTCAACCCGAACAGAAAAACAATGTTTAGGGCAGAGTAATAACCGAAGTTTTAAGGCACGATGAGATTAAACTAACTATaaa                                                                          >  2:2229147/1‑99 (MQ=255)
                      aGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATaaac                                                                         <  2:528591/100‑1 (MQ=255)
                          gCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATaaacaaac                                                                     <  2:445601/100‑1 (MQ=255)
                            aaatAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACgg                                                                   <  2:1155459/97‑1 (MQ=255)
                             gTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGt                                                                  <  2:1304488/100‑1 (MQ=255)
                              ttAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTc                                                                 <  2:749644/100‑1 (MQ=255)
                                  aGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGgc                                                             <  1:2513098/100‑1 (MQ=255)
                                     tGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCt                                                          <  2:2073645/100‑1 (MQ=255)
                                     tGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCt                                                          >  1:1651402/1‑100 (MQ=255)
                                      ggTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCtt                                                         <  1:369643/100‑1 (MQ=255)
                                        tcatcaCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTaa                                                       >  2:382394/1‑100 (MQ=255)
                                         catcaCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAg                                                      <  2:12066/100‑1 (MQ=255)
                                         catcaCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAg                                                      <  2:2091380/100‑1 (MQ=255)
                                         catcaCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGCTAATTTCAATATAAACACACGGTCCGGCGCTTaa                                                       >  1:499008/1‑99 (MQ=255)
                                           tcaCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGaa                                                    <  2:958094/100‑1 (MQ=255)
                                           tcaCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGaa                                                    >  1:2123179/1‑100 (MQ=255)
                                           tcaCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGaa                                                    >  1:240654/1‑100 (MQ=255)
                                            caCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAAt                                                   <  1:1702920/100‑1 (MQ=255)
                                            caCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAAt                                                   >  1:778817/1‑100 (MQ=255)
                                             aCGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTATCCGGCACGACGGGGTTAATATAAATATAAACAAACGGTCCGGCGCGTAAGAATg                                                  >  2:1270494/1‑100 (MQ=255)
                                              cGAACAGAAAAACAATGTTTAGCGCAAATTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATga                                                 <  2:1734134/100‑1 (MQ=255)
                                               gAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATgag                                                <  1:389344/100‑1 (MQ=255)
                                                aaCAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATgaga                                               >  1:1337422/1‑100 (MQ=255)
                                                aaCAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATgaga                                               <  1:745630/100‑1 (MQ=255)
                                                aaCAGAAAAACAAGGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCGCGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATgaga                                               <  1:501266/100‑1 (MQ=255)
                                                  cAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAAc                                             <  2:950847/100‑1 (MQ=255)
                                                   aGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACa                                            <  1:305373/100‑1 (MQ=255)
                                                    gAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACAt                                           >  2:456343/1‑100 (MQ=255)
                                                      aaaaCAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATcc                                         >  1:1541158/1‑100 (MQ=255)
                                                      aaaaCAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATcc                                         >  1:1541240/1‑100 (MQ=255)
                                                      aaaaCAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATcc                                         <  2:1337413/100‑1 (MQ=255)
                                                       aaaCAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATccc                                        >  1:134088/1‑100 (MQ=255)
                                                       aaaCAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATccc                                        >  2:2142976/1‑100 (MQ=255)
                                                         aCAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCgg                                      >  2:1966707/1‑100 (MQ=255)
                                                           aaTGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCa                                    >  2:2207691/1‑100 (MQ=255)
                                                              gTTTAGCGCAAAGTAATAACCTAAGTTTTCCGGCACGACGTCGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCaa                                   >  2:600418/1‑98 (MQ=255)
                                                              gTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAg                                 >  1:1774271/1‑100 (MQ=255)
                                                              gTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAg                                 <  2:1149146/100‑1 (MQ=255)
                                                              gTTTAGCGCAAAGGAATGACCAAAAATTCCCGGCACGACGTCGTTAATTTCAATATAAACAACCGGTCCGGCGCTTAAGAATGAGAACATCCCGGCcaa                                  >  2:600400/1‑99 (MQ=255)
                                                                ttAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCt                               >  2:1273688/1‑100 (MQ=255)
                                                                ttAGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCt                               <  1:1263995/100‑1 (MQ=255)
                                                                  aGCGCAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGa                             <  1:553948/100‑1 (MQ=255)
                                                                   gcgcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGaa                            <  1:2107590/100‑1 (MQ=255)
                                                                   gcgcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGaa                            <  2:610923/100‑1 (MQ=255)
                                                                   gcgcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGaa                            >  1:1196094/1‑100 (MQ=255)
                                                                    cgcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGaaa                           >  1:1085111/1‑100 (MQ=255)
                                                                    cgcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTGCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGaaa                           >  2:775048/1‑100 (MQ=255)
                                                                     gcAAAGTAATAACCAGAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAAc                          <  1:1425486/100‑1 (MQ=255)
                                                                     gcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAAc                          >  2:1479800/1‑100 (MQ=255)
                                                                     gcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAAc                          >  1:1704894/1‑100 (MQ=255)
                                                                     gcAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAAc                          <  1:39920/100‑1 (MQ=255)
                                                                     gcAAAGTAATAAACAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAAc                          <  2:1342630/100‑1 (MQ=255)
                                                                      cAAAGTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAAcc                         >  1:2009278/1‑100 (MQ=255)
                                                                        aaatcataaCCAAACTTTTCCCACACGACCTGGTTACTTTCAATCTAAACCACCGGTCAGGCGCTTACGAATAAGCACATCCCGGCAAAGCTGAAACCAc                       <  2:212138/95‑1 (MQ=255)
                                                                          gTAATAACCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTg                     >  1:765216/1‑100 (MQ=255)
                                                                            aataaCCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGc                   >  2:2162614/1‑100 (MQ=255)
                                                                            aataaCCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGc                   <  2:2178043/100‑1 (MQ=255)
                                                                               aaCCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAg                >  1:1048532/1‑100 (MQ=255)
                                                                               aaCCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAg                <  2:134086/100‑1 (MQ=255)
                                                                               aaCCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAg                <  1:2440057/100‑1 (MQ=255)
                                                                                aCCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGc               >  2:711235/1‑100 (MQ=255)
                                                                                aCCAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGc               <  1:253364/100‑1 (MQ=255)
                                                                                 ccAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCa              <  1:1485467/100‑1 (MQ=255)
                                                                                 ccAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCa              >  2:329939/1‑100 (MQ=255)
                                                                                 ccAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCa              >  1:1311175/1‑100 (MQ=255)
                                                                                 ccAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCa              >  1:2458704/1‑100 (MQ=255)
                                                                                  cAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCa              >  1:274438/1‑99 (MQ=255)
                                                                                  cAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCAt             >  2:2252531/1‑100 (MQ=255)
                                                                                  cAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCAt             >  2:1854271/1‑100 (MQ=255)
                                                                                  cAAAGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCAt             >  2:848145/1‑100 (MQ=255)
                                                                                   aaaGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATg            <  1:898373/100‑1 (MQ=255)
                                                                                   aaaGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATg            <  2:1651390/100‑1 (MQ=255)
                                                                                    aaGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGt           <  1:420554/100‑1 (MQ=255)
                                                                                     aGTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa          <  2:154053/100‑1 (MQ=255)
                                                                                      gTTTTCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAg         >  2:1519563/1‑100 (MQ=255)
                                                                                         ttCCGGCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTc      <  2:2271666/100‑1 (MQ=255)
                                                                                             ggCACGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGa  >  2:449464/1‑100 (MQ=255)
                                                                                                |                                                                                                
AACGCGCCAATGCGGCGGACGAAGCGGCTGTTAAAGATGGTCATCACGAACAGAAAAACAATGTTTAGCGCAAAGTAATAACCAAAGTTTTCCGGCGCGACGTGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGA  >  NC_000913/2278925‑2279117

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