Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,122,580 (A)6→7 intergenic (‑84/+50) yceB ← / ← grxB hypothetical protein/glutaredoxin 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,122,5741.A98.9% 320.6 / ‑3.4 89intergenic (‑78/+56)yceB/grxBhypothetical protein/glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (52/36);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 105 total aligned reads displayed.

TAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTT  >  NC_000913/1122482‑1122661
                                                                                             |                                                                                         
tAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCcagcag                                                                                            <  1:994414/93‑1 (MQ=255)
tAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCcagcag                                                                                            >  2:181970/1‑93 (MQ=255)
 aaaaaCTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGa                                                                                          >  2:1975212/1‑93 (MQ=255)
  aaaaCTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaa                                                                                        >  1:3166107/1‑93 (MQ=255)
  aaaaCTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaa                                                                                        >  2:3386659/1‑93 (MQ=255)
  aaaaCTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaa                                                                                        <  2:4648485/93‑1 (MQ=255)
   atgcTTGTTCATGTTAGTCCGCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaaa                                                                                       <  2:3339096/90‑1 (MQ=255)
       ttGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaaaaaaa                                                                                   >  2:1632473/1‑92 (MQ=255)
        tGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaaaaaaac                                                                                  <  2:4557409/93‑3 (MQ=255)
         gTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCGTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaaaaaaact                                                                                 >  1:3388394/1‑90 (MQ=255)
         gTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaaaaaaact                                                                                 <  1:770948/93‑4 (MQ=255)
          ttCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGaaaaaaactg                                                                                >  2:2083265/1‑89 (MQ=255)
             aTGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGaaa                                                                              >  2:2269006/1‑93 (MQ=255)
             aTGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGaaa                                                                              >  1:1894761/1‑93 (MQ=255)
              tGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGaaaa                                                                             <  1:1505755/93‑1 (MQ=255)
              tGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGaaaa                                                                             >  1:4069485/1‑93 (MQ=255)
              tGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGaaaa                                                                             >  1:2306108/1‑93 (MQ=255)
                 tAGTCCGCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAAcgg                                                                          <  1:1650197/93‑1 (MQ=255)
                  aGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAAcggc                                                                         <  1:585446/93‑1 (MQ=255)
                   gTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAAcggcg                                                                        >  2:852881/1‑93 (MQ=255)
                      cTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCa                                                                     >  1:3002322/1‑93 (MQ=255)
                       tCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCaa                                                                    >  2:1506350/1‑93 (MQ=255)
                       tCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCaa                                                                    >  2:710167/1‑93 (MQ=255)
                       tCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCaa                                                                    >  2:2798392/1‑93 (MQ=255)
                       tCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCaa                                                                    <  2:3818416/93‑1 (MQ=255)
                       tCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCaa                                                                    >  2:1732305/1‑93 (MQ=255)
                       tCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAGAAAAACTGAAAACGGCGGCaa                                                                    >  1:2043344/1‑93 (MQ=255)
                         caCATGAAATCTGCGTCAATAGCCTGACACAACGCAGCATGTGTCACCAGCGATAACCTCGCCAGCAGAAAAAAACGTAAATCGGCGGCAAcc                                                                  <  2:1975296/92‑1 (MQ=255)
                             tGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCcgcg                                                              >  1:531060/1‑93 (MQ=255)
                             tGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCcgcg                                                              >  1:77313/1‑93 (MQ=255)
                                 aTCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCCaaaa                                                          >  1:133557/1‑91 (MQ=255)
                                 aTCTGCGTAAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAATCTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATa                                                          <  1:3355697/93‑1 (MQ=255)
                                    tGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAg                                                       >  2:4589413/1‑93 (MQ=255)
                                    tGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAg                                                       <  1:2331453/93‑1 (MQ=255)
                                     gCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAgg                                                      >  1:1116801/1‑93 (MQ=255)
                                      cGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAAACCGCGAATACAGGc                                                     >  2:1903233/1‑93 (MQ=255)
                                        tCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTg                                                   <  1:1456510/93‑1 (MQ=255)
                                        tCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTg                                                   >  1:1062925/1‑93 (MQ=255)
                                          aaTATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGGGGACACCCGCGAATACAGGGTGcc                                                 >  1:22722/1‑93 (MQ=255)
                                           atatCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCcg                                                >  1:2827720/1‑93 (MQ=255)
                                           atatCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACCGCGGCAACCCGCGAATACACGCTGCcg                                                >  1:2142875/1‑93 (MQ=255)
                                            tatCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGccac                                               >  2:2238846/1‑91 (MQ=255)
                                             atCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCcgcg                                              >  1:3425973/1‑93 (MQ=255)
                                              tCCTGACACAACGCAGCATGTGTCACCAGCGATACACTCGCCAGCAGAAAAAAACTGAAAACGGCGGTAACCCGCGAATACTGGCTGCcgggc                                             >  2:2872326/1‑90 (MQ=255)
                                              tCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCgcgg                                             >  2:714368/1‑93 (MQ=255)
                                                cTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCgcggcg                                           >  1:2672231/1‑93 (MQ=255)
                                                 tGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGcggcgg                                          <  2:1703611/93‑1 (MQ=255)
                                                 tGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGcggcgg                                          <  2:1390388/93‑1 (MQ=255)
                                                 tGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGcggcgg                                          >  2:138175/1‑93 (MQ=255)
                                                   acacAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAAACCGCGAATACAGGCTGCCGCGGCgggg                                        >  1:2127525/1‑92 (MQ=255)
                                                     acaACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCa                                      >  1:1214367/1‑93 (MQ=255)
                                                     acaACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCa                                      >  1:3030084/1‑93 (MQ=255)
                                                     acaACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCa                                      >  2:1398224/1‑93 (MQ=255)
                                                       aaCGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAgg                                    <  2:1584047/93‑1 (MQ=255)
                                                       aaCGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAgg                                    >  2:1589418/1‑93 (MQ=255)
                                                        aCGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGa                                   <  1:3403093/93‑1 (MQ=255)
                                                         cGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGAt                                 <  2:2733905/93‑1 (MQ=255)
                                                         cGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGAt                                 <  2:3880848/93‑1 (MQ=255)
                                                         agcagcaTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGAt                                 <  1:2458213/92‑1 (MQ=255)
                                                          gcagcaTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGAtt                                >  1:712383/1‑93 (MQ=255)
                                                            agcaTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTaa                              <  1:46593/93‑1 (MQ=255)
                                                             gcaTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTGAGGATTaaa                             >  1:1523340/1‑93 (MQ=255)
                                                             gcaTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTaaa                             <  2:1061875/93‑1 (MQ=255)
                                                              caTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAAt                            <  2:1248286/93‑1 (MQ=255)
                                                              caTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAAt                            <  1:1506350/93‑1 (MQ=255)
                                                               aTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATc                           <  1:4164933/93‑1 (MQ=255)
                                                               aTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATc                           <  2:2750740/93‑1 (MQ=255)
                                                               aTGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGGATAAATc                           >  2:4249549/1‑93 (MQ=255)
                                                                  tgtCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGGTTAAATCGcc                        >  1:465839/1‑93 (MQ=255)
                                                                   gtCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCa                       <  2:1896701/93‑1 (MQ=255)
                                                                   gtCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGGTTAAATCGCCa                       >  1:1796519/1‑93 (MQ=255)
                                                                    tCACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCAt                      <  1:250592/93‑1 (MQ=255)
                                                                     cACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCAtt                     >  2:2255876/1‑93 (MQ=255)
                                                                     cACCAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCAtt                     <  2:3356974/93‑1 (MQ=255)
                                                                       ccAGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATtga                   <  2:1126340/93‑1 (MQ=255)
                                                                         aGCGATAAACTCGCCAGCAGAAAAAAACTGAGAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATtgatg                 <  1:2937637/93‑1 (MQ=255)
                                                                         aGCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGATTTAATCGCCATtgatg                 >  2:374257/1‑93 (MQ=255)
                                                                          gCGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATtgatga                <  1:777914/93‑1 (MQ=255)
                                                                           cGATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTgatgat               <  2:52425/93‑1 (MQ=255)
                                                                            gATAAACTCGCGAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATa              <  1:1361352/93‑1 (MQ=255)
                                                                            gATAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGATTAAATCGCCATTGATGATa              >  2:4596475/1‑93 (MQ=255)
                                                                             aTAAACTCGCCAGCAGAAAAAAACTGAAACCGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATaa             <  2:4373489/93‑1 (MQ=255)
                                                                             aTAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATaa             <  1:1685923/93‑1 (MQ=255)
                                                                             aTAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATaa             <  1:3182534/93‑1 (MQ=255)
                                                                              tAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCCTTGATGATAAc            >  1:2184037/1‑93 (MQ=255)
                                                                              tAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATaaa            >  2:1327198/1‑92 (MQ=255)
                                                                              tAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAAc            <  1:3454026/93‑1 (MQ=255)
                                                                              tAAACTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAAc            >  1:4279241/1‑93 (MQ=255)
                                                                               aaaCTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGGTTAAATCGCCATTGATGATAACa           >  2:3330977/1‑93 (MQ=255)
                                                                                aaCTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGAGTTAAACCCCCATGATTGATAACa           >  2:2500342/1‑93 (MQ=255)
                                                                                  cTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAAt        <  1:2773066/93‑1 (MQ=255)
                                                                                  cTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGATTTAATCGCCATTGGTGATAACAAAt        >  1:763863/1‑93 (MQ=255)
                                                                                  cTCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGATTAAATCGCCATTGATGATAACAAAt        >  1:3343023/1‑93 (MQ=255)
                                                                                   tCGCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGGTTTAAACGCCATTGATGATAACAAAtt       >  1:3249154/1‑93 (MQ=255)
                                                                                     gCCAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGa     <  2:2308773/93‑1 (MQ=255)
                                                                                      ccAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGGTTAAATCGCCATTGATGATAACAAATTTAt    >  2:1283266/1‑93 (MQ=255)
                                                                                      ccAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGAt    <  1:4473591/93‑1 (MQ=255)
                                                                                      ccAGCAGAAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGAt    >  2:2220239/1‑93 (MQ=255)
                                                                                       cagcagaAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGAtt   <  2:3874272/93‑1 (MQ=255)
                                                                                        agcagaAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGGCAGGATTAAATCGCCATTGATGATAACAAATTGAttt  >  2:812321/1‑93 (MQ=255)
                                                                                             |                                                                                         
TAAAAACTTGTTCATGGTAGTCCTCGACATGAAATCTGCGTCAATATCCTGACACAACGCAGCATGTGTCACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTT  >  NC_000913/1122482‑1122661

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