Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,425,743 Δ1 bp 5.8% pseudogene (114/453 nt) ydaY → pseudogene, Rac prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,425,7430C. 5.8% 230.5 / 7.4 69pseudogene (114/453 nt)ydaYpseudogene, Rac prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base . (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAACGAGGGATCGCGCTGCTTAAG  >  NC_000913/1425595‑1425884
                                                                                                                                                    |                                                                                                                                             
ggTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAACTGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTcc                                                                                                                                               <  3:90666/149‑1 (MQ=255)
   cATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAg                                                                                                                                            >  1:1347/1‑149 (MQ=255)
    aTTGTTCTGGCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAgg                                                                                                                                           <  5:128955/149‑1 (MQ=255)
      tGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTg                                                                                                                                         >  8:168325/1‑149 (MQ=255)
            gtctgGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCt                                                                                                                                   >  7:17374/1‑149 (MQ=255)
                 gCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAg                                                                                                                                            <  5:173499/135‑1 (MQ=255)
                 gCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAg                                                                                                                                            >  6:173499/1‑135 (MQ=255)
                  ccTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtgat                                                                                                                             >  8:156266/1‑149 (MQ=255)
                   cTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtgat                                                                                                                             >  1:102174/1‑148 (MQ=255)
                    tGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtgatga                                                                                                                           >  7:249137/1‑149 (MQ=255)
                      aCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCTAAGGTTTCCGAGGTGGATTCTGTTGATGAtt                                                                                                                         >  8:315656/1‑149 (MQ=255)
                      aCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGAtt                                                                                                                         <  1:16830/149‑1 (MQ=255)
                             aGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATAAATCCGTCGCCACATTGAGAAAATTGATGCCAGTTTGCACGTCGGAACGAATGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCAt                                                                                                                  <  6:278197/149‑1 (MQ=255)
                                 gAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCCCATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGa                                                                                                               >  5:192305/1‑148 (MQ=255)
                                       gcGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGcc                                                                                                        >  1:67409/1‑149 (MQ=255)
                                       gcGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGcc                                                                                                        >  5:100558/1‑149 (MQ=255)
                                        cGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGAAGATTTACTCATTGATAATGCCg                                                                                                       >  5:213161/1‑149 (MQ=255)
                                         gATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGc                                                                                                      <  8:186652/149‑1 (MQ=255)
                                             aCGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGt                                                                                                  >  2:249808/1‑149 (MQ=255)
                                              cGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGtt                                                                                                 <  7:168325/149‑1 (MQ=255)
                                               gATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTg                                                                                                                                         <  2:151254/108‑1 (MQ=255)
                                               gATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTg                                                                                                                                         >  1:151254/1‑108 (MQ=255)
                                                   aGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATct                                                                                             <  2:329333/148‑1 (MQ=255)
                                                   aGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATct                                                                                             <  2:82313/148‑1 (MQ=255)
                                                     tATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATctgct                                                                                          >  8:93130/1‑149 (MQ=255)
                                                      aTCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATctgct                                                                                          >  3:213408/1‑148 (MQ=255)
                                                         gCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGaaa                                                                                      <  8:17374/149‑1 (MQ=255)
                                                              cGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTg                                                                                 >  3:30235/1‑149 (MQ=255)
                                                                     cATTAATCCGTCGCCACATTGAGAAAAAGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAggg                                                                          <  1:120588/149‑1 (MQ=255)
                                                                      aTTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGGTTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTTGAAAgggg                                                                         >  7:79070/1‑149 (MQ=255)
                                                                       ttAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtgat                                                                                                                             <  6:211949/96‑1 (MQ=255)
                                                                       ttAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtgat                                                                                                                             >  5:211949/1‑96 (MQ=255)
                                                                        tAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTTGAAAGGGGtt                                                                       >  5:164765/1‑149 (MQ=255)
                                                                          aTCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTgg                                                                     <  4:213408/149‑1 (MQ=255)
                                                                           tCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGt                                                                    <  6:192305/149‑1 (MQ=255)
                                                                              gTCGCCACATTGAGAAAATGGATGCCAGTTAGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGTTGGATTCTGGTGATGATTTACTCATTGATAATTCAGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGaa                                                                 >  3:11076/1‑149 (MQ=255)
                                                                               tCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAAc                                                                >  4:139609/1‑149 (MQ=255)
                                                                                   cacaTTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGt                                                            <  6:52644/149‑1 (MQ=255)
                                                                                     caTTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGtt                                                           >  2:122219/1‑148 (MQ=255)
                                                                                       ttGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAAt                                                        <  7:156266/149‑1 (MQ=255)
                                                                                       ttGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAAt                                                        >  8:310295/1‑149 (MQ=255)
                                                                                         gagaCAATGGATGCCAGTGTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGGTTCCGAGGTGGATTCTGTTGATGAGTTACTCATTGATACTGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATg                                                       <  5:59017/148‑1 (MQ=255)
                                                                                         gagaAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCTAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATg                                                       <  7:315656/148‑1 (MQ=255)
                                                                                         gagaAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATg                                                       <  7:93130/148‑1 (MQ=255)
                                                                                         gagaAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATg                                                       <  2:67409/148‑1 (MQ=255)
                                                                                          agaAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGgt                                                     <  8:226418/149‑1 (MQ=255)
                                                                                             aaaTGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGtt                                                  <  6:164765/149‑1 (MQ=255)
                                                                                               aTGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGa                                                <  2:1347/149‑1 (MQ=255)
                                                                                                   aTGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGtt                                            <  1:249808/149‑1 (MQ=255)
                                                                                                      ccAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCa                                         <  4:23487/149‑1 (MQ=255)
                                                                                                         gTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTg                                      <  8:249137/149‑1 (MQ=255)
                                                                                                             gCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGaa                                   >  4:5534/1‑148 (MQ=255)
                                                                                                              cACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAat                                  <  6:100558/148‑1 (MQ=255)
                                                                                                               aCGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAatat                                <  3:139609/149‑1 (MQ=255)
                                                                                                                cGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAAtata                               <  3:5534/149‑1 (MQ=255)
                                                                                                                  tCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTTGTTAATGTTGTTGAGGTTGCACTGGAATATACg                             >  6:312774/1‑149 (MQ=255)
                                                                                                                         gAAGGAGTTTGATATTTCATAGGTTTC‑‑AGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGa                                                <  1:258499/121‑1 (MQ=255)
                                                                                                                         gAAGGAGTTTGATATTTCAAAGGTTTC‑‑AGGTGTATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTTGTGAACTGGTTAATGGTGTTGa                                                >  2:253961/1‑121 (MQ=255)
                                                                                                                         gAAGGAGTTTGATATTTCAAAGGTTTC‑‑AGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGa                                                <  1:253961/121‑1 (MQ=255)
                                                                                                                         gAAGGAGTTTGATATTTCAAAGGTTTC‑‑AGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGa                                                >  2:258499/1‑121 (MQ=255)
                                                                                                                          aaGGAGTTTGATATTTCAAAGGTGTCCGAGGTGGATTCTGTTGATGATTTACTCATTGAGAATGCCGCTCGTTATCTGCTGAACGACTGGAAAGGGGTTGGGGAAGTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAAc                      >  2:214149/1‑148 (MQ=255)
                                                                                                                           aggaggTTGATATTTCAAAGGTTTCCGCGGTGGATTAGGTTGATGATTTACTCATTGATACTGCCACTCGTTATCTGCTGCAAGACTGGAACGGGGTTGGTGAACTGGTTAATGGTGATGAGGTTGCACTGGAATATACGGCAGAACGa                    <  5:35205/149‑1 (MQ=255)
                                                                                                                                   gATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGGGAACTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAACGAGGGATcg             <  5:312774/148‑1 (MQ=255)
                                                                                                                                    atatTTCAAAGGTTTCCGAGGTGTATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTCGAAAGGGGTTGTTGAACTGGTTAATCGTGTTGAGGTTGCACTGGAATATACGTCAGAACGAGGGATcgcg           >  5:182924/1‑149 (MQ=255)
                                                                                                                                     tatTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAACTAGGGATCgcgc          >  5:120120/1‑149 (MQ=255)
                                                                                                                                        ttCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAACGAGGGATCGCgctgc       >  4:118109/1‑149 (MQ=255)
                                                                                                                                          cAAAGGTTTCCGAGGTGGATTCTGTTCATGATTTACTCAATGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGTTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAACGAGGGATCGCgctgat      >  5:83345/1‑146 (MQ=255)
                                                                                                                                           aaaGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAACGAGGGATCGCGCTGCtt     <  1:122219/148‑1 (MQ=255)
                                                                                                                                           aaaGGTTTCCGAGGTGGATTCTGTAGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGTAAAGGGGTTGTTGAACTCGTTAATCGTCTTGAGGTTGCACTGGAATATACGGCAGAACGAGGGATCGCGCTGCtt     >  5:92510/1‑148 (MQ=255)
                                                                                                                                                tttCCGAGGTTGATTCTGTTGATGATTTAATCATTGATAATTCAGCTCGTTATCTGCTGAATGACTGTAAAGGGTTTGGTGAACTGGTTAATGGTGTTGATGTTGAACTTGAATATACTGCAGAAAGAGGGATCGAGCTTATTAAg  >  5:266974/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTATCTGCTGAAAGACTGGAAAGGGGTTGGTGAACTGGTTAATGGTGTTGAGGTTGCACTGGAATATACGGCAGAACGAGGGATCGCGCTGCTTAAG  >  NC_000913/1425595‑1425884

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