Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 561,551 C→T 5.7% Q619* (CAA→TAA)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,5510CT 5.7% 228.2 / 2.0 69Q619* (CAA→TAA) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base T (2/2);  total (34/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.61e-01

TGTAGGGATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGG  >  NC_000913/561409‑561694
                                                                                                                                              |                                                                                                                                               
tGTAGGGATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATa                                                                                                                                           >  2:353805/1‑149 (MQ=255)
 gTAGGGATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAg                                                                                                                                          >  3:199261/1‑149 (MQ=255)
      gATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAgctggc                                                                                                                                     >  3:223223/1‑149 (MQ=255)
      gATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAgctggc                                                                                                                                     >  3:156536/1‑149 (MQ=255)
        tCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAgctggctg                                                                                                                                   >  4:205904/1‑149 (MQ=255)
            cGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAgc                                                                                                                                         >  3:13110/1‑139 (MQ=255)
            cGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAgc                                                                                                                                         <  4:13110/139‑1 (MQ=255)
            cGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAg                                                                                                                               <  3:145511/149‑1 (MQ=255)
              aTAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGgt                                                                                                                             >  2:124409/1‑149 (MQ=255)
                 acgaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTa                                                                                                                          <  4:181527/149‑1 (MQ=255)
                   gaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAgg                                                                                                                        <  4:379749/149‑1 (MQ=255)
                   gaacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAgg                                                                                                                        >  3:379749/1‑149 (MQ=255)
                    aacgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGt                                                                                                                       >  3:157604/1‑149 (MQ=255)
                      cgTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTgg                                                                                                                     <  3:407622/149‑1 (MQ=255)
                        tATTGTCGGACTTAATGTTTCAGTGCCTTTCAACGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGa                                                                                                                   >  2:194333/1‑149 (MQ=255)
                         aTTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGAc                                                                                                                  <  2:309436/149‑1 (MQ=255)
                           tGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCt                                                                                                                <  4:381142/149‑1 (MQ=255)
                                gACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTg                                                                                                           <  4:156536/149‑1 (MQ=255)
                                gACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTg                                                                                                           <  4:223223/149‑1 (MQ=255)
                                  cTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGTAAATTAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGa                                                                                                         >  4:322615/1‑149 (MQ=255)
                                  cTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGTAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGa                                                                                                         >  4:2559/1‑149 (MQ=255)
                                   ttAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGaa                                                                                                        <  1:314202/149‑1 (MQ=255)
                                   ttAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGaa                                                                                                        <  3:205904/149‑1 (MQ=255)
                                    tAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAg                                                                                                       <  4:337662/149‑1 (MQ=255)
                                    tAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAg                                                                                                       <  1:194333/149‑1 (MQ=255)
                                      aTGTTTCCGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCAGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGc                                                                                                     <  1:159168/149‑1 (MQ=255)
                                                    ttCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTa                                                                                       <  3:386786/149‑1 (MQ=255)
                                                      cAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGTAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATc                                                                                     <  3:322615/149‑1 (MQ=255)
                                                      cAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGTAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATc                                                                                     <  3:2559/149‑1 (MQ=255)
                                                       aaTGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCa                                                                                    >  3:444379/1‑149 (MQ=255)
                                                        aTGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCAc                                                                                   <  1:130248/149‑1 (MQ=255)
                                                            tttGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTaa                                                                               <  4:157604/149‑1 (MQ=255)
                                                               gACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGcc                                                                            <  4:444379/149‑1 (MQ=255)
                                                                aCCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCa                                                                           <  1:353805/149‑1 (MQ=255)
                                                                  cAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGcc                                                                                                    >  3:255294/1‑122 (MQ=255)
                                                                  cAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGcc                                                                                                    <  4:255294/122‑1 (MQ=255)
                                                                   aaaCGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAg                                                                                          <  4:396923/131‑1 (MQ=255)
                                                                   aaaCGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAg                                                                                          >  3:396923/1‑131 (MQ=255)
                                                                    aaCGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGt                                                                       <  1:248968/149‑1 (MQ=255)
                                                                       gTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAggg                                                                        <  1:441242/145‑1 (MQ=255)
                                                                       gTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAggg                                                                        >  2:441242/1‑145 (MQ=255)
                                                                          gCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACATTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATAcc                                                                 >  3:114652/1‑149 (MQ=255)
                                                                                  cGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGCAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCtt                                                                                                               >  3:249949/1‑95 (MQ=255)
                                                                                  cGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGCAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCtt                                                                                                               <  4:249949/95‑1 (MQ=255)
                                                                                   gggAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGt                                                                                                                       <  3:218253/86‑1 (MQ=255)
                                                                                   gggAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGt                                                                                                                       >  4:218253/1‑86 (MQ=255)
                                                                                             gcTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGgtgt                                                                                                                           <  3:369864/72‑1 (MQ=255)
                                                                                             gcTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGgtgt                                                                                                                           >  4:369864/1‑72 (MQ=255)
                                                                                             gcTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGgtgt                                                                                                                           >  2:87946/1‑72 (MQ=255)
                                                                                             gcTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGgtgt                                                                                                                           <  1:87946/72‑1 (MQ=255)
                                                                                                    cgcgCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGa                                                                     >  1:82255/1‑119 (MQ=255)
                                                                                                    cgcgCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGa                                                                     <  2:82255/119‑1 (MQ=255)
                                                                                                    cgcgCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAgcgc                                       >  4:191254/1‑149 (MQ=255)
                                                                                                         ttATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGt                                                                                                                       <  3:349584/64‑1 (MQ=255)
                                                                                                         ttATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGt                                                                                                                       >  4:349584/1‑64 (MQ=255)
                                                                                                         ttATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGAt                                                                    <  4:297827/115‑1 (MQ=255)
                                                                                                         ttATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGAt                                                                    >  3:297827/1‑115 (MQ=255)
                                                                                                         ttATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACgg                                  >  4:359079/1‑149 (MQ=255)
                                                                                                               ctcctTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACt                            >  2:305102/1‑149 (MQ=255)
                                                                                                                 cctTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTgg                          <  4:199261/149‑1 (MQ=255)
                                                                                                                  ctTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGc                         >  2:148991/1‑149 (MQ=255)
                                                                                                                   tttAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAAc                                                                                               <  3:411775/78‑1 (MQ=255)
                                                                                                                   tttAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAAc                                                                                               >  4:411775/1‑78 (MQ=255)
                                                                                                                   tttAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCa                        >  2:70463/1‑149 (MQ=255)
                                                                                                                   tttAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCa                        >  1:405692/1‑149 (MQ=255)
                                                                                                                   tttAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCa                        >  4:350564/1‑149 (MQ=255)
                                                                                                                   tttAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCa                        >  2:41250/1‑149 (MQ=255)
                                                                                                                            aaCCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTa               >  3:116767/1‑149 (MQ=255)
                                                                                                                                         aCGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg  <  1:124409/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TGTAGGGATCCCCGATAACGAACGTATTGTCGGACTTAATGTTTCAGTGCCTTTCAATGTTTTGACCAAACGTCGCTACACCCGGGAAAATGCGCTCGACCGCGCTTATGCCTCCTTTAACGCCAACCGTAACAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGG  >  NC_000913/561409‑561694

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