Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,664,972 T→C 5.4% F156S (TTT→TCT)  dmsD → twin‑argninine leader‑binding protein for DmsA and TorA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,664,9720TC 5.4% 230.7 / 3.1 74F156S (TTT→TCT) dmsDtwin‑argninine leader‑binding protein for DmsA and TorA
Reads supporting (aligned to +/- strand):  ref base T (37/33);  new base C (2/2);  total (39/35)
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.46e-01

GGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGAT  >  NC_000913/1664831‑1665118
                                                                                                                                             |                                                                                                                                                  
ggATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTg                                                                                                                                              >  2:287064/1‑148 (MQ=255)
       cgAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAaca                                                                                                                                       >  7:291660/1‑148 (MQ=255)
        gagaAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAacac                                                                                                                                      <  7:394723/148‑1 (MQ=255)
              ggCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTc                                                                                                                                >  2:447571/1‑148 (MQ=255)
                 aTTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTg                                                                                                                             <  4:383956/148‑1 (MQ=255)
                  ttCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGa                                                                                                                            >  1:230675/1‑148 (MQ=255)
                   tCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCTCTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGAt                                                                                                                           >  3:308780/1‑148 (MQ=255)
                    cAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATg                                                                                                                          >  5:57533/1‑148 (MQ=255)
                         tGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGtttt                                                                                                                      <  6:351705/147‑1 (MQ=255)
                           aaaTGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTAt                                                                                                                   <  1:447571/148‑1 (MQ=255)
                                aaGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGGCGCCAGACAGAATGCGAAGAATTGATGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGaa                                                                                                               <  1:25467/147‑1 (MQ=255)
                                    aaaaCGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAg                                                                                                           >  5:208427/1‑147 (MQ=255)
                                      aaCGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGttt                                                                                                                       >  8:98109/1‑133 (MQ=255)
                                      aaCGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGttt                                                                                                                       <  7:98109/133‑1 (MQ=255)
                                      aaCGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAg                                                                                                           >  5:280920/1‑145 (MQ=255)
                                      aaCGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAg                                                                                                           >  2:376516/1‑145 (MQ=255)
                                       aCGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGa                                                                                                       <  2:230675/148‑1 (MQ=255)
                                       aCGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGa                                                                                                       <  6:208427/148‑1 (MQ=255)
                                           aCCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAacac                                                                                                   >  5:180245/1‑148 (MQ=255)
                                              ggAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCa                                                                                                                                          <  1:16399/106‑1 (MQ=255)
                                              ggAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCa                                                                                                                                          >  2:16399/1‑106 (MQ=255)
                                              ggAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTCTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCt                                                                                                >  5:421164/1‑148 (MQ=255)
                                              ggAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTCTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCt                                                                                                >  5:427609/1‑148 (MQ=255)
                                                aaGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCttt                                                                                              <  6:57533/148‑1 (MQ=255)
                                                  gATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCt                                                                                            <  2:222845/148‑1 (MQ=255)
                                                          tGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGc                                                                                     <  6:280920/147‑1 (MQ=255)
                                                               ctctgctgTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTgg                                                                                <  4:308780/145‑1 (MQ=255)
                                                                   gctgTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGt                                                                                                                                               <  8:384288/80‑1 (MQ=255)
                                                                   gctgTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGt                                                                                                                                               >  7:384288/1‑80 (MQ=255)
                                                                   gctgTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGaa                                                                           >  5:341430/1‑148 (MQ=255)
                                                                     ggTTGATGGCTGGGGTGTTTGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCTCGTTTTTCCGTGGTCCACACGTTTTCTTGATGTTTTTATCGATAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACt                                                                         <  6:47789/147‑1 (MQ=255)
                                                                        tGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGc                                                                      >  3:87863/1‑148 (MQ=255)
                                                                                    ggTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCt                                                          >  8:96785/1‑148 (MQ=255)
                                                                                     gTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTg                                                         >  8:367005/1‑148 (MQ=255)
                                                                                     gTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTg                                                         <  5:378402/148‑1 (MQ=255)
                                                                                      ttGGCGGAGAATGGTCTCCATACAGAATGCGAAGAATTGCTGTCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATTTTTTTATCGAAAAAGCAGAACACCCTTTCTACAGTGCACTGGGTGAACTGGCTCGCTTAACGCTgg                                                        <  1:287064/148‑1 (MQ=255)
                                                                                       tGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGc                                                                                     <  6:309352/118‑1 (MQ=255)
                                                                                       tGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGc                                                                                     >  5:309352/1‑118 (MQ=255)
                                                                                       tGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCa                                                                                    <  1:250016/119‑1 (MQ=255)
                                                                                       tGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCa                                                                                    >  2:250016/1‑119 (MQ=255)
                                                                                            gagaATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAAt                                                  <  7:367005/148‑1 (MQ=255)
                                                                                             agaATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATg                                                 <  6:180245/148‑1 (MQ=255)
                                                                                               aaTGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTATTTCCGTCGTCAACACCTTTTCTTCATGTTTTTATCGAAAAAGCAGAACCCCCATTCTAACGTGCACTGGGTAAACTGTCTCGCTTAAAGCAGGCCCAAAgg                                                >  1:203586/1‑147 (MQ=255)
                                                                                                  ggTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACt                                                                                  >  3:388678/1‑110 (MQ=255)
                                                                                                  ggTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACt                                                                                  <  4:388678/110‑1 (MQ=255)
                                                                                                      gCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATcaca                                        <  4:87863/148‑1 (MQ=255)
                                                                                                       ccAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTCGCTCGCTTAACGCTGGCGCAATGGCAATCACaa                                       >  2:139583/1‑148 (MQ=255)
                                                                                                           acagaATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGc                                                                                     <  8:412062/98‑1 (MQ=255)
                                                                                                           acagaATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGc                                                                                     >  7:412062/1‑98 (MQ=255)
                                                                                                             agaCTGCGATGAATTGCTGGCATGGCACCTTTTTCCGTGGTGAACTCGTTTTCTTGCTTTTTTTATCGTAAAAGCCGAACACCGTTTCTACCGATCACTGGGTGAACTGGCGCGCTTATCGCTGGCGCAATGGCAATCACAACTGtt                                  <  8:201216/147‑1 (MQ=255)
                                                                                                                aTGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAAtt                              >  2:432054/1‑148 (MQ=255)
                                                                                                                aTGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAAtt                              >  4:88609/1‑148 (MQ=255)
                                                                                                                  gCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGc                                               >  8:24182/1‑129 (MQ=255)
                                                                                                                  gCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGc                                               <  7:24182/129‑1 (MQ=255)
                                                                                                                     aagaaTTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTAt                                                                                                                   <  2:207759/58‑1 (MQ=255)
                                                                                                                     aagaaTTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTAt                                                                                                                   >  1:207759/1‑58 (MQ=255)
                                                                                                                         aTTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTg                                    <  6:86978/133‑1 (MQ=255)
                                                                                                                         aTTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTg                                    >  5:86978/1‑133 (MQ=255)
                                                                                                                         aTTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTcgcg                     <  6:341430/148‑1 (MQ=255)
                                                                                                                             cTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGtt                  >  6:16364/1‑147 (MQ=255)
                                                                                                                                gcatggcaCCTTTCTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAAc              <  6:427609/148‑1 (MQ=255)
                                                                                                                                gcatggcaCCTTTCTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAAc              <  6:421164/148‑1 (MQ=255)
                                                                                                                                    ggcaCCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTaa                                                              <  7:174433/96‑1 (MQ=255)
                                                                                                                                    ggcaCCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTaa                                                              >  8:174433/1‑96 (MQ=255)
                                                                                                                                    ggcaCCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCa                                                    <  5:259261/106‑1 (MQ=255)
                                                                                                                                    ggcaCCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCa                                                    >  6:259261/1‑106 (MQ=255)
                                                                                                                                    ggcaCCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCt          >  1:283373/1‑148 (MQ=255)
                                                                                                                                      caCCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGt        >  3:87917/1‑148 (MQ=255)
                                                                                                                                       aCCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGtt       <  5:16364/148‑1 (MQ=255)
                                                                                                                                        ccTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTcgc                      <  4:2742/132‑1 (MQ=255)
                                                                                                                                        ccTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTcgc                      >  3:2742/1‑132 (MQ=255)
                                                                                                                                        ccTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGttt      >  8:194029/1‑148 (MQ=255)
                                                                                                                                        ccTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGttt      >  4:55339/1‑148 (MQ=255)
                                                                                                                                         cTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACATCTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTc     >  4:354660/1‑148 (MQ=255)
                                                                                                                                            tttCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCg    <  1:432054/146‑1 (MQ=255)
                                                                                                                                             ttCCGTTTTCAACAAGTTTTCTTGATGTTTTTATTGAAAAAGCAGAACAGACTTTCTACCGTGCACTGGGTGAACTGGCTAGCTTAACGCTGTCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGAt  <  1:139583/147‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGAT  >  NC_000913/1664831‑1665118

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: ‑