Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 893,111 A→T 5.4% D52V (GAT→GTT)  ybjN → negative regulator of motility; multicopy suppressor of coaA(Ts)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913893,1110AT 5.4% 236.7 / 2.2 74D52V (GAT→GTT) ybjNnegative regulator of motility; multicopy suppressor of coaA(Ts)
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base T (2/2);  total (37/37)
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

TCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGCCAAAGCTGGTGGTT  >  NC_000913/892963‑893259
                                                                                                                                                    |                                                                                                                                                    
tCGCTGGGCGTTCCTGGTCTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGa                                                                                                                                                      <  6:33603/149‑1 (MQ=255)
   cTGGTCGTTCCTGGTCTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATaa                                                                                                                                                   <  8:54075/149‑1 (MQ=255)
      gTCGTTCCTGGTCTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAac                                                                                                                                                  <  8:64999/147‑1 (MQ=255)
                gTCTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGtt                                                                                                                                       <  4:296344/148‑1 (MQ=255)
                 tCTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGtttt                                                                                                                                     <  8:124331/149‑1 (MQ=255)
                 tCTAGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGtttt                                                                                                                                     >  4:218659/1‑149 (MQ=255)
                  cTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGTATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTc                                                                                                                                    >  1:311550/1‑149 (MQ=255)
                   tGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCt                                                                                                                                   >  2:331982/1‑149 (MQ=255)
                           cTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGACAACACGATCCTGTTTTCTGCCATGGc                                                                                                                           <  1:38071/149‑1 (MQ=255)
                            tGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATaa                                                                                                                                                   <  7:274954/124‑1 (MQ=255)
                            tGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATaa                                                                                                                                                   >  8:274954/1‑124 (MQ=255)
                               gTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGTCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCCTGGCGGaa                                                                                                                       <  5:95684/149‑1 (MQ=255)
                                  aaTGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGGTCCTGTTTTCTGCCATGGCGGAAGTc                                                                                                                    <  8:303618/149‑1 (MQ=255)
                                   aTGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTcc                                                                                                                   >  2:83879/1‑149 (MQ=255)
                                    tGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCg                                                                                                                  <  8:70785/149‑1 (MQ=255)
                                      gCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGCCGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGAc                                                                                                                >  6:282696/1‑149 (MQ=255)
                                      gCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGAc                                                                                                                >  6:126206/1‑149 (MQ=255)
                                      gCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCAGCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGa                                                                                                                 <  7:302875/148‑1 (MQ=255)
                                              aCCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGTTAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGcc                                                                                                        >  2:174647/1‑149 (MQ=255)
                                              aCCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGTTAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGcc                                                                                                        >  2:169634/1‑149 (MQ=255)
                                                cTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCg                                                                                                       >  3:254672/1‑148 (MQ=255)
                                                         aGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCt                                                                                                                                          <  8:71116/104‑1 (MQ=255)
                                                           tttttttGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTgg                                                                                           >  3:299920/1‑149 (MQ=255)
                                                                ttGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTggcggc                                                                                       <  8:312411/148‑1 (MQ=255)
                                                                   aaTGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGa                                                                                    <  5:180724/148‑1 (MQ=255)
                                                                   aaTGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGAt                                                                                   <  4:214717/149‑1 (MQ=255)
                                                                         aTAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCt                                                                             <  7:226329/149‑1 (MQ=255)
                                                                         aTAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCt                                                                             <  3:255646/149‑1 (MQ=255)
                                                                          tAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTg                                                                            <  6:139544/149‑1 (MQ=255)
                                                                            aCTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGcc                                                                          <  4:254672/149‑1 (MQ=255)
                                                                             cTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCa                                                                         <  2:264267/149‑1 (MQ=255)
                                                                             cTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCa                                                                         >  8:212909/1‑149 (MQ=255)
                                                                              tGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCAt                                                                        >  7:248961/1‑149 (MQ=255)
                                                                                  aGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAAt                                                                    >  2:61656/1‑149 (MQ=255)
                                                                                  aGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGGGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAAt                                                                    <  5:339204/149‑1 (MQ=255)
                                                                                     ctctGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGc                                                                  >  2:31471/1‑148 (MQ=255)
                                                                                      tctGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGc                                                                  >  8:1962/1‑147 (MQ=255)
                                                                                           aTCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCg                                                           >  8:94862/1‑149 (MQ=255)
                                                                                                ccccATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGa                                                       <  3:218659/148‑1 (MQ=255)
                                                                                                ccccATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGAc                                                      >  3:208192/1‑149 (MQ=255)
                                                                                                 cccATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGAcc                                                     <  4:229810/149‑1 (MQ=255)
                                                                                                  gcATATGCAGAATTTCGACGGTGTCTTTGATTCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTg                                                                                                                                  <  2:249436/70‑1 (MQ=255)
                                                                                                   cATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGt                                                   >  8:312743/1‑149 (MQ=255)
                                                                                                        gCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAg                                              >  3:128644/1‑149 (MQ=255)
                                                                                                          aGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCa                                            >  2:77450/1‑149 (MQ=255)
                                                                                                              attCCACGGTGACTTAGATGCCAAAATCGATCTGATCGATAACACGATACTCTTTTCTGCCATGGCGGAAGTCCGCCCTTCAGCCGTATTGCCGATGGCGGCGGTTTTATCTGCCATAAAAGCCAGTTAGCTGCCCGAGAAAGCATATc                                        >  6:272385/2‑149 (MQ=255)
                                                                                                                tCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCtt                                      >  3:162092/1‑149 (MQ=255)
                                                                                                                  gACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTg                                     <  8:248961/148‑1 (MQ=255)
                                                                                                                  gACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTAGCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGa                                    <  1:83879/149‑1 (MQ=255)
                                                                                                                   aCGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGat                                   <  1:61656/149‑1 (MQ=255)
                                                                                                                     ggTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGatat                                 >  4:224950/1‑149 (MQ=255)
                                                                                                                        gtCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATg                                >  2:178755/1‑147 (MQ=255)
                                                                                                                           tttGATGCCAAAATCGATCTGATCGATGACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGga                           >  5:269238/1‑149 (MQ=255)
                                                                                                                           tttGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGga                           >  1:179823/1‑149 (MQ=255)
                                                                                                                           tttGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGga                           >  2:261033/1‑149 (MQ=255)
                                                                                                                             tGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGgat                          <  7:212909/148‑1 (MQ=255)
                                                                                                                              gATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATg                                <  1:37968/141‑1 (MQ=255)
                                                                                                                              gATGCCAAAATCGATATGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGACGCTGGAGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATg                                >  2:37968/1‑141 (MQ=255)
                                                                                                                                tGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATa                      <  5:126206/149‑1 (MQ=255)
                                                                                                                                  ccAAAATCGATCTGATCGTTAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATaa                     <  1:174647/148‑1 (MQ=255)
                                                                                                                                  ccAAAATCGATCTGATCGTTAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATaa                     <  1:169634/148‑1 (MQ=255)
                                                                                                                                  ccAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAAt                    >  4:324516/1‑149 (MQ=255)
                                                                                                                                   cAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATc                   <  4:299920/149‑1 (MQ=255)
                                                                                                                                     aaaTCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTg                 <  1:31471/149‑1 (MQ=255)
                                                                                                                                        tCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAGTCTGCCa              <  4:208192/149‑1 (MQ=255)
                                                                                                                                        tCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGcc               >  5:41463/1‑148 (MQ=255)
                                                                                                                                           aTCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATATGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGCCAAAg           >  6:137644/1‑149 (MQ=255)
                                                                                                                                           aTCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGCGGCGGATTTATCt                                                                             <  6:7406/82‑1 (MQ=255)
                                                                                                                                           aTCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGCGGCGGATTTATCt                                                                             >  5:7406/1‑82 (MQ=255)
                                                                                                                                              tGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGCCAAAGCtg        >  7:174321/1‑149 (MQ=255)
                                                                                                                                              tGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGCCAAAGCtg        >  5:154997/1‑149 (MQ=255)
                                                                                                                                                aTCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGCCAAAGCtggt      >  2:228855/1‑149 (MQ=255)
                                                                                                                                                 tCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCaa                                                                     <  2:346638/85‑1 (MQ=255)
                                                                                                                                                 tCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCaa                                                                     >  1:346638/1‑85 (MQ=255)
                                                                                                                                                 tCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGCCAAAGCtggtg     <  2:179823/149‑1 (MQ=255)
                                                                                                                                                    aTAACTCGATCCTGTTTTCTGCCATGGCGGAAGGCCTACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCATTTCGCTGACCGTGAAAGCATTTCTTGTTATGCAGGTTGTTAATCTGCCAAAGCTGGTGGtt  <  2:311550/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGCTGGTCGTTCCTGGTCTGGATACGCTGCGTCAATGGCTCGATGACCTGGGGATGAGTTTTTTTGAATGTGATAACTGTCAGGCTCTGCATCTGCCCCATATGCAGAATTTCGACGGTGTCTTTGATGCCAAAATCGATCTGATCGATAACACGATCCTGTTTTCTGCCATGGCGGAAGTCCGACCTTCAGCCGTATTGCCGCTGGCGGCGGATTTATCTGCCATCAATGCCAGTTCGCTGACCGTGAAAGCATTTCTTGATATGCAGGATGATAATCTGCCAAAGCTGGTGGTT  >  NC_000913/892963‑893259

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: