Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,471 T→A 53.0% L674Q (CTG→CAG)  zntA → zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,4710TA53.0% 38.6 / 312.2 217L674Q (CTG→CAG) zntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base T (41/61);  new base A (55/60);  total (96/121)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 218 total aligned reads displayed.

CGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGT  >  NC_000913/3608373‑3608569
                                                                                                  |                                                                                                  
cGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTg                                                                                                   >  1:1014973/1‑100 (MQ=255)
cGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTg                                                                                                   >  1:1792866/1‑100 (MQ=255)
 gggATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTgg                                                                                                  <  1:2670340/100‑1 (MQ=255)
 gggATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAgg                                                                                                  <  2:889995/100‑1 (MQ=255)
 gggATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAgg                                                                                                  <  2:2820007/100‑1 (MQ=255)
 gggATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAgg                                                                                                  <  2:4089318/100‑1 (MQ=255)
  ggATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGc                                                                                                 <  1:1210895/100‑1 (MQ=255)
  ggATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGc                                                                                                 >  1:3047553/1‑100 (MQ=255)
      tGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCAcgc                                                                                             >  1:2494834/1‑100 (MQ=255)
       gCAATGGGTAGCGGCACAGACGTGGCGCTGGAGACCGCCTACTCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCAcgcg                                                                                            <  2:2897453/100‑1 (MQ=255)
       gCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCAcgcg                                                                                            <  1:1277562/100‑1 (MQ=255)
        cAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACgcgc                                                                                           >  1:1355609/1‑100 (MQ=255)
         aaTGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCCAATGATTGAACTGGCACGCGcc                                                                                          <  2:1072854/100‑1 (MQ=255)
         aaTGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGcc                                                                                          <  1:1659835/100‑1 (MQ=255)
         aaTGCGTAGCGGCCCAGACGTGGCGCCGCAAACCGCCGACGCAGCATTAACCCATACCCACCTGCGCGGCCTCGTGCACATGATTAAACTCCCACGCGcc                                                                                          <  2:3507423/100‑1 (MQ=255)
             ggTAGCGGCACAGACGTGGCGCTGGAAACAGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTc                                                                                      <  2:2474830/100‑1 (MQ=255)
              gTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCa                                                                                     >  1:352312/1‑100 (MQ=255)
               tAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATCGAACTGGCACGCGCCACTCAc                                                                                    >  2:2174440/1‑100 (MQ=255)
                aGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCGGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACg                                                                                   <  2:2232795/100‑1 (MQ=255)
                 gCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGc                                                                                  >  1:4024017/1‑100 (MQ=255)
                  cGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCCCGcc                                                                                 >  2:2453989/1‑100 (MQ=255)
                    gCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCaa                                                                               >  2:3868610/1‑100 (MQ=255)
                     cacaGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATCATTGAACTGGCACGCGCCACTCACGCCAat                                                                              <  2:2727264/100‑1 (MQ=255)
                      acaGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAata                                                                             <  2:1200420/100‑1 (MQ=255)
                      acaGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAata                                                                             <  2:1014967/100‑1 (MQ=255)
                       caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAatat                                                                            >  1:387755/1‑100 (MQ=255)
                       caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAatat                                                                            >  2:2845375/1‑100 (MQ=255)
                       caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAatat                                                                            >  2:1480871/1‑100 (MQ=255)
                       caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAatat                                                                            >  2:1798537/1‑100 (MQ=255)
                       caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAatat                                                                            >  1:1380489/1‑100 (MQ=255)
                         gACGTGGCGCTGGACACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTCGCACAAATGATTGAACAAGCACGCGCCACTCACGCCAATATcc                                                                          >  2:815306/1‑100 (MQ=255)
                         gACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATcc                                                                          >  1:3185677/1‑100 (MQ=255)
                          aCGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCg                                                                         <  1:23805/100‑1 (MQ=255)
                           cGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGc                                                                        <  1:1865578/100‑1 (MQ=255)
                           cGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGc                                                                        >  2:3393503/1‑100 (MQ=255)
                           cGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGc                                                                        <  1:1114654/100‑1 (MQ=255)
                            gTGGCGCTGGAACCCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGcc                                                                       <  2:2581147/100‑1 (MQ=255)
                                 gatGGAAACCGCCGACGCAGCATTAACCCATAACAACCTGCGCCGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAAc                                                                  <  1:829167/98‑1 (MQ=255)
                                  cTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACa                                                                 <  1:2959631/100‑1 (MQ=255)
                                   tGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATAGCCGCCAGAACAt                                                                >  2:38809/1‑100 (MQ=255)
                                    ggAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACAt                                                                >  2:3471827/1‑99 (MQ=255)
                                    ggAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATc                                                               <  1:747536/100‑1 (MQ=255)
                                     gAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCa                                                              <  1:2633057/100‑1 (MQ=255)
                                      aaaCCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAACTGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCAc                                                             <  2:1950149/100‑1 (MQ=255)
                                      aaaCCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCAc                                                             <  2:3951798/100‑1 (MQ=255)
                                          cgccgACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTAtt                                                         >  1:519451/1‑100 (MQ=255)
                                            ccgACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTgc                                                       <  2:344895/100‑1 (MQ=255)
                                              gACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGTGCCACTCACGCCAATATCCGCCAGAACATCACTATTgcgc                                                     >  1:1325472/1‑100 (MQ=255)
                                              gACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTgcgc                                                     <  2:3047534/100‑1 (MQ=255)
                                               aCGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCt                                                    <  1:2819127/100‑1 (MQ=255)
                                                cGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTg                                                   <  2:4023991/100‑1 (MQ=255)
                                                 gcagcaTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTgg                                                  >  1:3280768/1‑100 (MQ=255)
                                                  cagcaTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTggg                                                 <  1:1653333/100‑1 (MQ=255)
                                                    gcaTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTggg                                                 >  2:1605982/1‑98 (MQ=255)
                                                    gcaTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGc                                               <  1:3352843/100‑1 (MQ=255)
                                                      aTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTg                                             >  1:1985637/1‑100 (MQ=255)
                                                         aaCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGaaa                                          >  2:2540164/1‑100 (MQ=255)
                                                         aaCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGaaa                                          >  2:2775676/1‑100 (MQ=255)
                                                          aCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAg                                         <  2:2445412/100‑1 (MQ=255)
                                                          aCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAg                                         <  1:1480883/100‑1 (MQ=255)
                                                           cccATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAgg                                        <  1:2907450/100‑1 (MQ=255)
                                                           cccATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAgg                                        <  1:1506834/100‑1 (MQ=255)
                                                           cccAGAGCCCCCTGCGCGGCCTAGGGCAAATGATTTAACAGGCGCGGGCGACTCACGCCAAGATCCGCCAGAACATCACTATTGCGCTGGGGc                                               <  2:3282846/93‑1 (MQ=255)
                                                              aTAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGAt                                     >  2:543558/1‑100 (MQ=255)
                                                              aTAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGAt                                     <  1:104564/100‑1 (MQ=255)
                                                              aTAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGAt                                     <  1:3113711/100‑1 (MQ=255)
                                                                aaCCGCCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATcn                                   <  2:3638112/100‑2 (MQ=255)
                                                                aaCCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCt                                   >  2:2050019/1‑100 (MQ=255)
                                                                aaCCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCt                                   <  2:1271626/100‑1 (MQ=255)
                                                                 accaccTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCtt                                  >  2:1696954/1‑100 (MQ=255)
                                                                 accaccTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCtt                                  >  2:1869302/1‑100 (MQ=255)
                                                                 accaccTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCtt                                  >  2:3541763/1‑100 (MQ=255)
                                                                   caccTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTTAAAGGGATCTTcc                                <  2:2558374/100‑1 (MQ=255)
                                                                    accTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCt                               >  1:68947/1‑100 (MQ=255)
                                                                    accTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACCTCACTATTGCGCTGGGGCTGAAAGGGATCTTcc                                >  2:1314702/1‑99 (MQ=255)
                                                                    accTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGagatcgga                               >  1:3023806/1‑93 (MQ=255)
                                                                    accTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCt                               >  1:2602967/1‑100 (MQ=255)
                                                                      cTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCg                             <  1:3224260/100‑1 (MQ=255)
                                                                       ttcgcgGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCt                               <  1:3714019/95‑1 (MQ=255)
                                                                         cgcgGCCGGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTca                          <  2:3372702/100‑1 (MQ=255)
                                                                          gcgGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCATTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTcac                         <  1:2922246/100‑1 (MQ=255)
                                                                             gCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTcaccac                      <  1:613149/100‑1 (MQ=255)
                                                                              ccTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACg                     >  1:3564619/1‑100 (MQ=255)
                                                                              ccTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACg                     >  1:2113263/1‑100 (MQ=255)
                                                                               cTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCGTCCTCGTCACCACGc                    >  1:4063851/1‑100 (MQ=255)
                                                                               cTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGc                    >  2:3035484/1‑100 (MQ=255)
                                                                                tGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCt                   <  1:3948080/100‑1 (MQ=255)
                                                                                 ggTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTg                  >  1:3825513/1‑100 (MQ=255)
                                                                                 ggTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTg                  <  2:494328/100‑1 (MQ=255)
                                                                                 ggTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTg                  >  1:3326897/1‑100 (MQ=255)
                                                                                   tGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGt                 >  2:1344816/1‑99 (MQ=255)
                                                                                      aaaTGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAgg             <  2:685579/100‑1 (MQ=255)
                                                                                      aaaTGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAgg             <  2:3161974/100‑1 (MQ=255)
                                                                                      aaaTGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAgg             <  1:2769289/100‑1 (MQ=255)
                                                                                      aaaTGATTGAAATGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAgg             <  1:232878/100‑1 (MQ=255)
                                                                                           aTTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGa        >  2:1668177/1‑100 (MQ=255)
                                                                                            ttGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGAc       >  1:4015017/1‑100 (MQ=255)
                                                                                               aaCAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGCTGACCg     >  2:846911/1‑99 (MQ=255)
                                                                                                 cTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGt  >  1:2144918/1‑100 (MQ=255)
                                                                                                  tGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCgg    >  2:2085632/1‑97 (MQ=255)
                                                                                                  |                                                                                                  
CGGGATTGCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGT  >  NC_000913/3608373‑3608569

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