Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,043,403 A→G H260H (CAT→CAC yihG ← putative endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,4030AG99.7% 2043.0 / ‑4.3 622H260H (CAT→CACyihGputative endonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (291/328);  total (292/330)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 642 total aligned reads displayed.

TGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAG  >  NC_000913/4043313‑4043492
                                                                                          |                                                                                         
tGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTg                                                                                         >  1:5945091/1‑93 (MQ=255)
tGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTg                                                                                         <  2:1679148/93‑1 (MQ=255)
 gAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACgtgt                                                                                        >  2:5904788/1‑93 (MQ=255)
  aGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTa                                                                                       <  2:4386635/93‑1 (MQ=255)
   gTTCAGCCATTGCTGAAAATGTCGCTTAAAGCTGTTGTCATTGATGTAGTCGCCATGTACCTCATCGGCAATAGGCTGTAAATCCACGTGTAc                                                                                      <  2:2117461/93‑1 (MQ=255)
      cAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGAc                                                                                   <  2:2505538/93‑1 (MQ=255)
        gCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACaa                                                                                 <  2:9299245/93‑1 (MQ=255)
        gCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACaa                                                                                 <  2:3074230/93‑1 (MQ=255)
          cATTGCTGAAAGTGGCGCTTAAAGCGTTTGTCATTGATGTAGTTGCCATGTAACTCATCGGCAATAGGCTGTCAATCCACGTGTACGACCATc                                                                               >  2:5170195/1‑93 (MQ=255)
            ttGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcg                                                                             >  1:9401537/1‑93 (MQ=255)
            ttGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcg                                                                             >  2:7008793/1‑93 (MQ=255)
            ttGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcg                                                                             >  2:8253216/1‑93 (MQ=255)
             tGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCcgc                                                                            <  1:2435721/93‑1 (MQ=255)
                tGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTc                                                                         <  1:2248556/93‑1 (MQ=255)
                tGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTc                                                                         <  1:5406830/93‑1 (MQ=255)
                 gAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCa                                                                        <  1:3465790/93‑1 (MQ=255)
                   aaaTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGt                                                                      <  1:6070391/93‑1 (MQ=255)
                   aaaTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGt                                                                      >  1:6442981/1‑93 (MQ=255)
                    aaTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGtt                                                                     >  2:8663588/1‑93 (MQ=255)
                    aaTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGtt                                                                     >  2:8921376/1‑93 (MQ=255)
                     aTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGttt                                                                    <  1:9551268/93‑1 (MQ=255)
                     aTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGttt                                                                    <  1:6775435/93‑1 (MQ=255)
                     aTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGttt                                                                    <  1:7222782/93‑1 (MQ=255)
                     aTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGttt                                                                    <  1:2836902/93‑1 (MQ=255)
                        gcgcTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTAcc                                                                 >  2:8067430/1‑93 (MQ=255)
                         cgcTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCg                                                                >  1:1363191/1‑93 (MQ=255)
                          gcTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGc                                                               >  1:3307484/1‑93 (MQ=255)
                          gcTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGc                                                               <  1:7026842/93‑1 (MQ=255)
                           cTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCt                                                              <  2:4529031/93‑1 (MQ=255)
                           cTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCt                                                              >  1:10018738/1‑93 (MQ=255)
                             tAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTa                                                            <  1:987218/93‑1 (MQ=255)
                              aaaGCTTTTGTCATTGATGTAGTAGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTaa                                                           <  2:5615420/93‑1 (MQ=255)
                               aaGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAAc                                                          >  2:2218484/1‑93 (MQ=255)
                                aGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACa                                                         <  1:5538480/93‑1 (MQ=255)
                                 gCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACat                                                        >  2:1800329/1‑93 (MQ=255)
                                  cTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACata                                                       <  1:1132796/93‑1 (MQ=255)
                                  cTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACata                                                       <  1:1920757/93‑1 (MQ=255)
                                   ttttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACatat                                                      <  1:2965562/93‑1 (MQ=255)
                                   ttttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACatat                                                      <  1:4417157/93‑1 (MQ=255)
                                    tttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATc                                                     >  2:4181404/1‑93 (MQ=255)
                                      tGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCga                                                   <  1:3625708/93‑1 (MQ=255)
                                       gTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaa                                                  <  1:1718872/93‑1 (MQ=255)
                                       gTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaa                                                  >  2:6904864/1‑93 (MQ=255)
                                        tCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaag                                                 <  1:5262813/93‑1 (MQ=255)
                                        tCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaag                                                 >  2:5730765/1‑93 (MQ=255)
                                        tCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGAGTCAGTTTACCGCTTAACATATCgaag                                                 <  1:6257017/93‑1 (MQ=255)
                                        tCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCAGCGTCAGTTTACCGCTTAACATATCgaag                                                 <  1:7303755/93‑1 (MQ=255)
                                         cATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaaga                                                <  2:5461951/93‑1 (MQ=255)
                                          aTTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaagaa                                               <  2:7355842/93‑1 (MQ=255)
                                           ttGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAAc                                              <  2:3459511/93‑1 (MQ=255)
                                           ttGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAAc                                              <  2:1278474/93‑1 (MQ=255)
                                            tGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACg                                             >  1:4836255/1‑93 (MQ=255)
                                             gATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACgg                                            >  2:5261753/1‑93 (MQ=255)
                                             gATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACgg                                            <  1:4973896/93‑1 (MQ=255)
                                              aTGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGc                                           >  1:7830049/1‑93 (MQ=255)
                                               tGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCt                                          <  2:3641250/93‑1 (MQ=255)
                                               tGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCt                                          >  2:6845164/1‑93 (MQ=255)
                                                gTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTg                                         <  1:5807691/93‑1 (MQ=255)
                                                 tAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGa                                        >  1:4193562/1‑93 (MQ=255)
                                                 tAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGa                                        >  1:7549823/1‑93 (MQ=255)
                                                  aGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGAc                                       <  2:7977581/93‑1 (MQ=255)
                                                  aGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGAc                                       <  2:7786338/93‑1 (MQ=255)
                                                   gTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACg                                      <  1:6904651/93‑1 (MQ=255)
                                                    tCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACgg                                     >  2:3224030/1‑93 (MQ=255)
                                                    tCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACgg                                     <  2:261051/93‑1 (MQ=255)
                                                      gCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGtt                                   >  1:1493979/1‑93 (MQ=255)
                                                       ccATGTAACTCATCAGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTa                                  >  2:6164933/1‑93 (MQ=255)
                                                        cATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAt                                 >  2:7153102/1‑93 (MQ=255)
                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTg                               >  2:7499910/1‑93 (MQ=255)
                                                            tAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTc                             <  2:3973448/93‑1 (MQ=255)
                                                            tAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTc                             <  1:8894846/93‑1 (MQ=255)
                                                              aCTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCg                           >  2:6404489/1‑93 (MQ=255)
                                                              aCTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCg                           <  2:6659023/93‑1 (MQ=255)
                                                                tCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGa                         <  2:1555550/93‑1 (MQ=255)
                                                                 cATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGAt                        <  2:9662107/93‑1 (MQ=255)
                                                                 cATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGAt                        >  1:570021/1‑93 (MQ=255)
                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                    <  2:490618/93‑1 (MQ=255)
                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                    <  2:9593818/93‑1 (MQ=255)
                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                    >  1:3830502/1‑93 (MQ=255)
                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                    <  1:7352379/93‑1 (MQ=255)
                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                    >  1:4672725/1‑93 (MQ=255)
                                                                         aTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGc                >  1:850636/1‑93 (MQ=255)
                                                                         aTAGGCTGTAAAGCCCCGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGc                <  1:9166604/93‑1 (MQ=255)
                                                                          tAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCg               <  1:6571482/93‑1 (MQ=255)
                                                                           aGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGt              <  2:5077887/93‑1 (MQ=255)
                                                                            ggCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTa             <  2:10033379/93‑1 (MQ=255)
                                                                            ggCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTa             >  2:6265975/1‑93 (MQ=255)
                                                                            ggCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTa             <  1:3204446/93‑1 (MQ=255)
                                                                              cTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGTATAACACAGCGTAAc           <  1:2607583/93‑1 (MQ=255)
                                                                              cTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAAc           <  2:8478366/93‑1 (MQ=255)
                                                                               tGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACa          <  1:124752/93‑1 (MQ=255)
                                                                                gTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACAt         >  1:8036281/1‑93 (MQ=255)
                                                                                  aaaTCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTa       <  2:2914566/93‑1 (MQ=255)
                                                                                   actCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTaa      <  2:107213/91‑1 (MQ=255)
                                                                                   aaTCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTaa      <  2:2291017/93‑1 (MQ=255)
                                                                                   aaTCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTaa      <  1:8619844/93‑1 (MQ=255)
                                                                                   aaTCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTaa      <  2:3781481/93‑1 (MQ=255)
                                                                                      ccACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCa   <  1:3948753/93‑1 (MQ=255)
                                                                                       cacgTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAg  <  1:8288/93‑1 (MQ=255)
                                                                                       cacgTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGGTATTGTCCGGATAACACCGCGTAACATTAAGCAg  >  1:1397573/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAG  >  NC_000913/4043313‑4043492

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