Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,561,105 A→G V60A (GTG→GCG)  yjiK ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,561,1050AG95.7% 660.5 / 14.1 231V60A (GTG→GCG) yjiKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/8);  new base G (71/150);  total (73/158)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

Only 100 of 240 total aligned reads displayed.

ATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATC  >  NC_000913/4561015‑4561194
                                                                                          |                                                                                         
aTTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCt                                                                                         >  2:768438/1‑93 (MQ=255)
aTTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCt                                                                                         >  1:4322631/1‑93 (MQ=255)
aTTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCt                                                                                         <  2:4211636/93‑1 (MQ=255)
aTTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCt                                                                                         >  1:1638702/1‑93 (MQ=255)
aTTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCt                                                                                         >  2:1538584/1‑93 (MQ=255)
   tCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGcc                                                                                      >  1:2203879/1‑93 (MQ=255)
    cTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCCGCca                                                                                     >  1:2532478/1‑93 (MQ=255)
     tttGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCcac                                                                                    >  1:1345949/1‑93 (MQ=255)
       tGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACaa                                                                                  >  1:4441607/1‑93 (MQ=255)
        gCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAAt                                                                                 <  1:671807/93‑1 (MQ=255)
        gCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAAt                                                                                 <  2:295343/93‑1 (MQ=255)
           gTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGc                                                                              <  1:208564/93‑1 (MQ=255)
           gTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGc                                                                              <  1:3654768/93‑1 (MQ=255)
            tCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCa                                                                             >  2:1463140/1‑93 (MQ=255)
               aTGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCaaaa                                                                          >  1:3502069/1‑93 (MQ=255)
               aTGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCaaaa                                                                          >  1:2468829/1‑93 (MQ=255)
               aTGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCaaaa                                                                          >  1:2583168/1‑93 (MQ=255)
                tggtggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCaaaaa                                                                         >  2:3931719/1‑93 (MQ=255)
                 ggtggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAg                                                                        <  2:2625759/93‑1 (MQ=255)
                 ggtggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAg                                                                        >  2:348815/1‑93 (MQ=255)
                   tggCGTGTTATTTCTGAAATGATGCTGCATGATTGCTTTTACGGGCACAAGATTGCACGAAGAAAGTACTCGCTGCCACAATGGCAAAAAGaa                                                                      <  2:3813930/93‑1 (MQ=255)
                   tggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGaa                                                                      <  1:1535864/93‑1 (MQ=255)
                    ggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAAt                                                                     >  1:4080358/1‑93 (MQ=255)
                    ggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAAt                                                                     >  1:99295/1‑93 (MQ=255)
                     gCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATc                                                                    <  2:821240/93‑1 (MQ=255)
                      cGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATca                                                                   >  2:3197877/1‑93 (MQ=255)
                      cGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATca                                                                   <  2:648119/93‑1 (MQ=255)
                      cGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATca                                                                   <  1:2827134/93‑1 (MQ=255)
                       gTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATcac                                                                  >  2:2681232/1‑93 (MQ=255)
                        tGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATcaca                                                                 >  2:3410213/1‑93 (MQ=255)
                         ggTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACaa                                                                >  2:4429102/1‑93 (MQ=255)
                         ggTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACaa                                                                <  2:169708/93‑1 (MQ=255)
                          gTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACaat                                                               >  1:3887136/1‑93 (MQ=255)
                            aaTTCTGAAATGTTGCTGCATGATTGCTTTTACGTACGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAGTCACaataa                                                             <  2:712143/93‑1 (MQ=255)
                                cTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGga                                                         <  1:527146/93‑1 (MQ=255)
                                cTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGga                                                         >  1:3605191/1‑93 (MQ=255)
                                 tGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGgag                                                        >  1:3721688/1‑93 (MQ=255)
                                  gAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAACAAGAATCACAATAATGgaga                                                       >  1:797449/1‑93 (MQ=255)
                                  gAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGgaga                                                       <  2:931552/93‑1 (MQ=255)
                                  gAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGgaga                                                       <  2:4667552/93‑1 (MQ=255)
                                  gAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGgaga                                                       <  2:1839909/93‑1 (MQ=255)
                                    aaTGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATg                                                     >  1:2354420/1‑93 (MQ=255)
                                    aaTGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATg                                                     >  1:4178654/1‑93 (MQ=255)
                                     atgatgCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGc                                                    >  1:2855793/1‑93 (MQ=255)
                                      tgatgCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCg                                                   >  2:3588766/1‑93 (MQ=255)
                                      tgatgCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCg                                                   >  2:3961614/1‑93 (MQ=255)
                                      tgatgCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCg                                                   <  2:2809208/93‑1 (MQ=255)
                                      tgatgCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCg                                                   <  1:4128427/93‑1 (MQ=255)
                                        atgCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGtt                                                 <  1:4483696/93‑1 (MQ=255)
                                          gctgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTa                                               <  1:304565/93‑1 (MQ=255)
                                           ctgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTAc                                              <  1:4525073/93‑1 (MQ=255)
                                           ctgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTAc                                              <  2:2434986/93‑1 (MQ=255)
                                           ctgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTAc                                              >  2:1906376/1‑93 (MQ=255)
                                            tgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACt                                             <  1:481838/93‑1 (MQ=255)
                                            tgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACt                                             >  2:2518877/1‑93 (MQ=255)
                                             gcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGCGATGCGTTTACtt                                            >  1:2372676/1‑93 (MQ=255)
                                             gcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACtt                                            <  2:742821/93‑1 (MQ=255)
                                               aTGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTaa                                          >  1:3996235/1‑93 (MQ=255)
                                               aTGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTaa                                          <  2:1169687/93‑1 (MQ=255)
                                                tGATTGCTTTTCCGGGCGCAAGATTGCACGCAGAAAGCACACGCTGCCACAATGGCAACAAGAATCACAATAATGGAGATGCGTTTACTTAAc                                         >  2:3136531/1‑93 (MQ=255)
                                                  aTTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGa                                       >  1:4590394/1‑93 (MQ=255)
                                                  aTTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGa                                       >  1:1731162/1‑93 (MQ=255)
                                                  aTTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGa                                       >  1:3783146/1‑93 (MQ=255)
                                                  aTTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGa                                       <  1:3367864/93‑1 (MQ=255)
                                                  aTTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATAACAATAATGGAGATGCGTTTACTTAACGa                                       >  2:2893214/1‑93 (MQ=255)
                                                    tGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGaaa                                     <  1:1491184/93‑1 (MQ=255)
                                                       ttttACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATAc                                  <  1:1261378/93‑1 (MQ=255)
                                                       ttttACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATAc                                  <  2:4115187/93‑1 (MQ=255)
                                                       ttttACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATAc                                  <  1:2002915/93‑1 (MQ=255)
                                                        tttACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACt                                 <  1:3000381/93‑1 (MQ=255)
                                                          tACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACttt                               <  2:541685/93‑1 (MQ=255)
                                                          tACGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACttt                               <  2:1173260/93‑1 (MQ=255)
                                                           aCGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACtttt                              <  1:3424827/93‑1 (MQ=255)
                                                           aCGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACtttt                              <  2:2177492/93‑1 (MQ=255)
                                                            cGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACttttt                             <  1:1915455/93‑1 (MQ=255)
                                                            cGGGCGCAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACttttt                             <  1:2763560/93‑1 (MQ=255)
                                                            cGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACttttt                             <  1:1468027/93‑1 (MQ=255)
                                                                cgcAAGATTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCa                         <  2:4369601/93‑1 (MQ=255)
                                                                    aGATTGCACGAAGAAATTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGt                     <  1:839776/93‑1 (MQ=255)
                                                                      aTTGCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTaa                   <  1:1725521/93‑1 (MQ=255)
                                                                        tGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATc                 <  1:1212414/93‑1 (MQ=255)
                                                                         gCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATcc                <  2:955337/93‑1 (MQ=255)
                                                                         gCACGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATcc                <  1:3926747/93‑1 (MQ=255)
                                                                           aCGAAGAAAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCt              <  2:2020884/93‑1 (MQ=255)
                                                                              aagaaAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTa           <  2:3893722/93‑1 (MQ=255)
                                                                               agaaAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTaa          <  2:3663012/93‑1 (MQ=255)
                                                                               agaaAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTaa          <  1:3340002/93‑1 (MQ=255)
                                                                               agaaAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTaa          <  2:3231612/93‑1 (MQ=255)
                                                                                gaaAGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAg         <  1:1009190/93‑1 (MQ=255)
                                                                                 aaaGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAgg        <  1:574007/93‑1 (MQ=255)
                                                                                  aaGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGc       <  2:484389/93‑1 (MQ=255)
                                                                                  aaGTACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGc       <  2:2383086/93‑1 (MQ=255)
                                                                                     tACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCaaa    <  2:434263/93‑1 (MQ=255)
                                                                                     tACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCaaa    <  2:2122306/93‑1 (MQ=255)
                                                                                     tACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCaaa    <  1:688224/93‑1 (MQ=255)
                                                                                     tACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCaaa    <  1:334448/93‑1 (MQ=255)
                                                                                     tACACGCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCaaa    <  1:1756058/93‑1 (MQ=255)
                                                                                     tACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCaaa    <  2:4033885/93‑1 (MQ=255)
                                                                                       cacgcTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATc  <  2:1985673/93‑1 (MQ=255)
                                                                                       cacgcTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATc  <  2:1594251/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATC  >  NC_000913/4561015‑4561194

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