Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,364,381 T→C Y221C (TAC→TGC)  dsbD ← thiol:disulfide interchange protein and activator of DsbC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,364,3810TC99.1% 362.0 / ‑6.6 105Y221C (TAC→TGC) dsbDthiol:disulfide interchange protein and activator of DsbC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (45/59);  total (46/59)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 106 total aligned reads displayed.

GTAAAGACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGTAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCAGAAATCA  >  NC_000913/4364285‑4364477
                                                                                                |                                                                                                
gTAAAGACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTGACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAgg                                                                                               <  1:2587582/100‑1 (MQ=255)
gTAAAGACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAgg                                                                                               >  1:1033799/1‑100 (MQ=255)
gTAAAGACGATGGCGAGGCCAAAGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAgg                                                                                               <  2:1657695/100‑1 (MQ=255)
  aaaGACGATGGCGAGGCCAATTAGCACGTATGGGTGCTGTAGCGCCGCCTGTAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTc                                                                                             <  1:620538/100‑1 (MQ=255)
  aaaGACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTc                                                                                             >  1:2563553/1‑100 (MQ=255)
  aaaGACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTc                                                                                             >  2:2291088/1‑100 (MQ=255)
  aaaGACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTc                                                                                             <  1:1717915/100‑1 (MQ=255)
     gACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAgc                                                                                          >  2:2256549/1‑100 (MQ=255)
      aCGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAgcg                                                                                         <  2:2978738/100‑1 (MQ=255)
       cGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAgcgc                                                                                        <  1:2448262/100‑1 (MQ=255)
       cGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAgcgc                                                                                        <  2:2832920/100‑1 (MQ=255)
        gATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGcc                                                                                       <  2:247471/100‑1 (MQ=255)
          tGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCCGCGCCa                                                                                      >  2:162956/1‑99 (MQ=255)
            gCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATcc                                                                                   <  1:612922/100‑1 (MQ=255)
            gCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATcc                                                                                   >  1:319553/1‑100 (MQ=255)
            gCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATcc                                                                                   >  1:1920496/1‑100 (MQ=255)
             cGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATccc                                                                                  >  2:1517111/1‑100 (MQ=255)
              gAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTTGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATcccc                                                                                 <  2:1078002/100‑1 (MQ=255)
              gAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATcccc                                                                                 <  1:161245/100‑1 (MQ=255)
              gAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATcccc                                                                                 <  2:900167/100‑1 (MQ=255)
                ggCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTg                                                                               <  1:2868956/100‑1 (MQ=255)
                ggCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTg                                                                               >  1:2745133/1‑100 (MQ=255)
                 gCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGc                                                                              <  1:1711191/100‑1 (MQ=255)
                 gCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGc                                                                              <  2:2416369/100‑1 (MQ=255)
                 gCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGc                                                                              <  1:2825719/100‑1 (MQ=255)
                 gCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGc                                                                              <  2:517588/100‑1 (MQ=255)
                    aaTGAGCACGTATGGGTGCTGTAGCACCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGcaca                                                                           <  1:2646273/100‑1 (MQ=255)
                        aGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACataaa                                                                       >  1:294927/1‑100 (MQ=255)
                           aCGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACataaataa                                                                    >  2:2557846/1‑100 (MQ=255)
                            cGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACAtaaataaa                                                                   >  1:1467764/1‑100 (MQ=255)
                              tATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAACgg                                                                 >  2:2379740/1‑100 (MQ=255)
                              tATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAgg                                                                 <  2:258978/100‑1 (MQ=255)
                              tATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAgg                                                                 <  1:1468078/100‑1 (MQ=255)
                              tATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAgg                                                                 <  1:2565457/100‑1 (MQ=255)
                              tATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAgg                                                                 <  1:220429/100‑1 (MQ=255)
                                tGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTc                                                               <  1:110891/100‑1 (MQ=255)
                                tGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTc                                                               <  2:1824536/100‑1 (MQ=255)
                                 gggTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGTAGGTCAGCGCCATCCCCTACACATAAATAAAGGTCa                                                              >  2:1151938/1‑100 (MQ=255)
                                 gggTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGCCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCa                                                              >  1:2132413/1‑100 (MQ=255)
                                   gTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGc                                                            <  1:1669004/100‑1 (MQ=255)
                                   gTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGc                                                            <  2:2563530/100‑1 (MQ=255)
                                    tGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGTGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGca                                                           <  1:2973893/100‑1 (MQ=255)
                                       tGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACCCAAATAAAGGGCAGcaa                                                          >  2:2874997/1‑98 (MQ=255)
                                       tGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGcaaca                                                        >  2:1908942/1‑100 (MQ=255)
                                          aGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATg                                                     <  2:516527/100‑1 (MQ=255)
                                          aGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATg                                                     >  2:2337669/1‑100 (MQ=255)
                                           gcgcCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGc                                                    <  2:906205/100‑1 (MQ=255)
                                           gcgcCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGc                                                    >  1:1572719/1‑100 (MQ=255)
                                             gccgccTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGctc                                                  >  2:2223785/1‑100 (MQ=255)
                                             gccgccTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGct                                                   >  1:2848864/1‑99 (MQ=255)
                                               cgccTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTg                                                <  1:2256569/100‑1 (MQ=255)
                                                gccTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTgg                                               <  2:1575854/100‑1 (MQ=255)
                                                gccTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTgg                                               >  1:130376/1‑100 (MQ=255)
                                                 ccTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGc                                              <  2:834580/100‑1 (MQ=255)
                                                  cTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCa                                             >  1:2979027/1‑100 (MQ=255)
                                                  cTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCa                                             >  1:282902/1‑100 (MQ=255)
                                                   tGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAg                                            >  1:753590/1‑100 (MQ=255)
                                                    ggAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCCGGTCAGCGCCATCCCCTGCACCGAAATAAAGGTCAGCAACAATGctcagaca                                             >  2:673932/1‑93 (MQ=255)
                                                      aaCTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTgg                                         <  1:2185813/100‑1 (MQ=255)
                                                      aaCTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTgg                                         <  2:2745109/100‑1 (MQ=255)
                                                        cTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGgag                                       >  1:226619/1‑100 (MQ=255)
                                                         tGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGgaga                                      <  2:1467752/100‑1 (MQ=255)
                                                             aCCCTGCGGCGGCAACCACCTGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCg                                  >  2:1854484/1‑100 (MQ=255)
                                                             aCCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCg                                  >  1:2557299/1‑100 (MQ=255)
                                                              cccTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCGGTGGAGAGCCGc                                 >  1:50842/1‑100 (MQ=255)
                                                              cccTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGc                                 <  2:1978217/100‑1 (MQ=255)
                                                              cccTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGc                                 >  1:440154/1‑100 (MQ=255)
                                                              cccTGCGGCGGCAACCACCAGACCCAACGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGc                                 >  1:1848271/1‑100 (MQ=255)
                                                               ccTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCt                                >  1:1047266/1‑100 (MQ=255)
                                                                cTGCGGCGGCAACCACCAGCCCCAGCGCCGTGCAAGTGCGGGCCATACCCCTCACGTAAATGAAGGCCAGCAACAATGCCCTGGCAGTCGacaggagc                                 >  2:2293069/1‑94 (MQ=255)
                                                                cTGCGGCGGCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTg                               >  1:833953/1‑100 (MQ=255)
                                                                   cggcggCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGttt                            <  2:1432711/100‑1 (MQ=255)
                                                                   cggcggCAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGttt                            <  2:1394926/100‑1 (MQ=255)
                                                                    ggcggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTa                           <  2:2879072/100‑1 (MQ=255)
                                                                    ggcggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTa                           <  2:1342562/100‑1 (MQ=255)
                                                                     gcggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTac                          <  2:690150/100‑1 (MQ=255)
                                                                     gcggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTac                          <  1:447119/100‑1 (MQ=255)
                                                                     gcggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTac                          >  2:1034322/1‑100 (MQ=255)
                                                                      cggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAGCAATGCTCTGGCAGTGGAGAGCCGCTGTTTacc                         >  1:1830002/1‑100 (MQ=255)
                                                                      cggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTacc                         <  2:2345464/100‑1 (MQ=255)
                                                                       gggaaCCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCGGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTacc                         <  2:323958/96‑1 (MQ=255)
                                                                       ggcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTacca                        <  1:2305613/100‑1 (MQ=255)
                                                                        gcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTaccac                       <  1:1255390/100‑1 (MQ=255)
                                                                        gcAACCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTaccac                       >  1:2503878/1‑100 (MQ=255)
                                                                          aaCCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCAccc                     <  1:1387169/100‑1 (MQ=255)
                                                                          aaCCACCAGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCAccc                     >  1:2501452/1‑100 (MQ=255)
                                                                            ccaccaGACCCAGCGCCGTGCAGGTCAGAGCCATCCCCTGCACATAAATATAGGTCGGCAACAAAGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCa                    <  2:3032216/99‑1 (MQ=255)
                                                                             caccaGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGc                  >  1:392080/1‑100 (MQ=255)
                                                                             caccaGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGc                  <  2:217989/100‑1 (MQ=255)
                                                                              accaGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCa                 >  1:40064/1‑100 (MQ=255)
                                                                               ccaGACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCAc                >  1:965902/1‑100 (MQ=255)
                                                                                  gACCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGAt             >  1:479995/1‑100 (MQ=255)
                                                                                   aCCCAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATg            <  1:2578384/100‑1 (MQ=255)
                                                                                      tggcatcGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCa         <  1:347960/94‑1 (MQ=255)
                                                                                      cAGCGCCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCa         <  1:2637361/100‑1 (MQ=255)
                                                                                          gcCGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCAGaaa     <  2:965893/100‑1 (MQ=255)
                                                                                            cGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCAGAAATc   <  1:2119754/100‑1 (MQ=255)
                                                                                            cGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCAGAAATc   <  2:392075/100‑1 (MQ=255)
                                                                                            cGTGCAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCAGAAATc   <  2:2149244/100‑1 (MQ=255)
                                                                                             gtgcAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCAGAAATCa  >  1:112326/1‑100 (MQ=255)
                                                                                                |                                                                                                
GTAAAGACGATGGCGAGGCCAATGAGCACGTATGGGTGCTGTAGCGCCGCCTGGAACTGTAACCCTGCGGCGGCAACCACCAGACCCAGCGCCGTGTAGGTCAGCGCCATCCCCTGCACATAAATAAAGGTCAGCAACAATGCTCTGGCAGTGGAGAGCCGCTGTTTACCACCCAGCACGATGCCAGAAATCA  >  NC_000913/4364285‑4364477

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: