Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 137,057 G→A intergenic (‑140/‑26) yacC ← / → yacK hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913137,0570GA99.6% 987.4 / ‑5.6 278intergenic (‑140/‑26)yacC/yacKhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (94/182);  total (94/184)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 301 total aligned reads displayed.

AGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGT  >  NC_000913/136967‑137146
                                                                                          |                                                                                         
agGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTTTTTCAcc                                                                                         >  1:2740453/1‑93 (MQ=255)
agGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCAcc                                                                                         <  1:3314728/93‑1 (MQ=255)
 ggCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCt                                                                                        <  1:6097601/93‑1 (MQ=255)
  gCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTg                                                                                       <  1:6088701/93‑1 (MQ=255)
    gacgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCt                                                                                     <  1:127909/93‑1 (MQ=255)
     acgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCtt                                                                                    >  2:4886249/1‑93 (MQ=255)
     acgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCtt                                                                                    >  1:1941351/1‑93 (MQ=255)
     acgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCtt                                                                                    <  2:1652979/93‑1 (MQ=255)
     acgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCtt                                                                                    <  1:3710597/93‑1 (MQ=255)
     acgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCtt                                                                                    <  1:5296468/93‑1 (MQ=255)
     acgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCtt                                                                                    >  1:6722112/1‑93 (MQ=255)
      cgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGTGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTa                                                                                   <  1:1599407/93‑1 (MQ=255)
       gacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTaa                                                                                  <  2:850534/93‑1 (MQ=255)
       gacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTaa                                                                                  >  2:3189519/1‑93 (MQ=255)
        acCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTTACCTGCTTAAg                                                                                 >  1:2361644/1‑93 (MQ=255)
           attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAaaaaa                                                                              >  2:3740579/1‑92 (MQ=255)
           attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAAt                                                                              <  1:2567257/93‑1 (MQ=255)
             tattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATaa                                                                            <  2:70578/93‑1 (MQ=255)
              attGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAg                                                                           <  2:4575352/93‑1 (MQ=255)
               ttGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAgg                                                                          >  2:6657561/1‑93 (MQ=255)
                tGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGa                                                                         <  1:6326082/93‑1 (MQ=255)
                tGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGa                                                                         <  2:1689807/93‑1 (MQ=255)
                 gTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGaa                                                                        <  1:4321680/93‑1 (MQ=255)
                 gTCAATAAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTCTGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGaa                                                                        <  1:5451137/93‑1 (MQ=255)
                  tCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGaaa                                                                       <  1:4297952/93‑1 (MQ=255)
                   cGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAAt                                                                      <  1:628945/93‑1 (MQ=255)
                    gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGGAGGACTATGCTCAACCTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATa                                                                     >  2:6320105/1‑93 (MQ=255)
                    gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATTCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAAATAGGAAATa                                                                     >  2:3362748/1‑93 (MQ=255)
                      ttAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAAc                                                                   >  1:2834551/1‑93 (MQ=255)
                       tAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACt                                                                  <  1:3951102/93‑1 (MQ=255)
                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTa                                                                 <  2:1334878/93‑1 (MQ=255)
                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTa                                                                 <  2:14622/93‑1 (MQ=255)
                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTa                                                                 <  2:3714306/93‑1 (MQ=255)
                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTa                                                                 <  1:4165890/93‑1 (MQ=255)
                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTa                                                                 <  1:4046219/93‑1 (MQ=255)
                         aaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTAt                                                                <  1:1100935/93‑1 (MQ=255)
                          aTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGGTTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATg                                                               >  2:2723723/1‑93 (MQ=255)
                           ttGTGTCTGCGGCTTGACCTTCCCGGAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAAATATgg                                                              >  1:638946/1‑92 (MQ=255)
                             gtgtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCaa                                                            <  2:5606627/93‑1 (MQ=255)
                             gtgtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCaa                                                            <  1:6171900/93‑1 (MQ=255)
                              tgtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAc                                                           <  1:3220260/93‑1 (MQ=255)
                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         <  2:3663564/93‑1 (MQ=255)
                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         <  2:1850776/93‑1 (MQ=255)
                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         <  1:4599643/93‑1 (MQ=255)
                                  tGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcg                                                       <  1:1351238/93‑1 (MQ=255)
                                    cGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCCACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTg                                                     >  2:3898255/1‑93 (MQ=255)
                                     ggCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGa                                                    <  2:4375573/93‑1 (MQ=255)
                                      gCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGAt                                                   >  1:550767/1‑93 (MQ=255)
                                       cTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGAtt                                                  <  1:2000230/93‑1 (MQ=255)
                                         tGACCTTCCCGTGAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTc                                                <  2:5213670/93‑1 (MQ=255)
                                         tGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTc                                                >  2:5107327/1‑93 (MQ=255)
                                         tGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTc                                                <  1:3076847/93‑1 (MQ=255)
                                          gACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCt                                               >  2:2523651/1‑93 (MQ=255)
                                          gACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCt                                               <  1:770083/93‑1 (MQ=255)
                                             cTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaa                                            >  2:5766238/1‑93 (MQ=255)
                                              ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                           >  1:5045513/1‑93 (MQ=255)
                                               tCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaaa                                          <  1:4772835/93‑1 (MQ=255)
                                                cccGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAat                                         <  1:3612804/93‑1 (MQ=255)
                                                 ccGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAata                                        >  1:5843799/1‑93 (MQ=255)
                                                  cGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAatat                                       >  2:2153516/1‑93 (MQ=255)
                                                  cGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAatat                                       <  1:3812895/93‑1 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      <  2:5937639/93‑1 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      <  1:2407611/93‑1 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      <  2:2455028/93‑1 (MQ=255)
                                                    tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                     >  2:3457663/1‑93 (MQ=255)
                                                    tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                     <  2:4400502/93‑1 (MQ=255)
                                                    tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                     <  1:2929660/93‑1 (MQ=255)
                                                     aaGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTcc                                    <  1:2123007/93‑1 (MQ=255)
                                                       ggggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGt                                  >  2:4135636/1‑93 (MQ=255)
                                                        gggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                 <  2:5518720/93‑1 (MQ=255)
                                                         ggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcg                                <  2:3524641/93‑1 (MQ=255)
                                                         ggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcg                                <  1:1432570/93‑1 (MQ=255)
                                                          gAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgc                               <  2:6164518/93‑1 (MQ=255)
                                                           aaGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgcg                              >  1:665118/1‑93 (MQ=255)
                                                            aGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCgcgc                             >  2:5846519/1‑93 (MQ=255)
                                                            aGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCgcgc                             <  2:344249/93‑1 (MQ=255)
                                                             ggACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCt                            >  2:4016528/1‑93 (MQ=255)
                                                               aCTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                          >  2:2185664/1‑93 (MQ=255)
                                                               aCTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                          <  1:3662536/93‑1 (MQ=255)
                                                                cTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTggg                         >  2:2940905/1‑93 (MQ=255)
                                                                  aTGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtg                       <  2:6527231/93‑1 (MQ=255)
                                                                   tGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtgt                      <  2:3794388/93‑1 (MQ=255)
                                                                    gCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGtgtg                     <  2:2058088/93‑1 (MQ=255)
                                                                      tCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGc                   <  1:5111355/93‑1 (MQ=255)
                                                                      tCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGc                   <  1:5247697/93‑1 (MQ=255)
                                                                      tCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGc                   >  1:3575426/1‑93 (MQ=255)
                                                                        aaCGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCtt                 <  2:2690431/93‑1 (MQ=255)
                                                                           gTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCgg              >  1:6485037/1‑93 (MQ=255)
                                                                            tttGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGc             >  2:1103417/1‑93 (MQ=255)
                                                                              tGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCtt           <  2:547022/93‑1 (MQ=255)
                                                                              tGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCtt           <  2:2781153/93‑1 (MQ=255)
                                                                              tGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCtt           <  1:5386984/93‑1 (MQ=255)
                                                                               gATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt          <  1:2245056/93‑1 (MQ=255)
                                                                               gATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt          >  2:2267941/1‑93 (MQ=255)
                                                                                 ttttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGc        <  1:3968483/93‑1 (MQ=255)
                                                                                 ttttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGATGGGTGTGGCTTCCGCTTTGc        >  1:963572/1‑93 (MQ=255)
                                                                                  tttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGcc       >  2:4074725/1‑93 (MQ=255)
                                                                                   ttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCg      >  1:4422415/1‑93 (MQ=255)
                                                                                     gTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCt    <  1:1153977/93‑1 (MQ=255)
                                                                                       ttCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt  >  2:1266547/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGT  >  NC_000913/136967‑137146

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