Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,664,796 G→T 100% V97V (GTG→GTT 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,7960GT100.0% 335.4 / ‑6.7 92V97V (GTG→GTTdmsDtwin‑argninine leader‑binding protein for DmsA and TorA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (44/47);  total (45/47)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCA  >  NC_000913/1664650‑1664936
                                                                                                                                                  |                                                                                                                                            
ccAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCt                                                                                                                                            <  1:853793/149‑1 (MQ=255)
 cAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTg                                                                                                                                           >  3:348191/1‑149 (MQ=255)
  aGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAAtctgttctgt                                                                                                                                          <  2:7182/149‑1 (MQ=255)
      gCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTgg                                                                                                                                      <  4:837331/149‑1 (MQ=255)
        gTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCg                                                                                                                                    >  1:704053/1‑149 (MQ=255)
         tCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGa                                                                                                                                   >  2:286686/1‑149 (MQ=255)
          cATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGAt                                                                                                                                  >  3:543367/1‑149 (MQ=255)
               gcgcCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAAc                                                                                                                             <  2:659070/149‑1 (MQ=255)
               gcgcCGCCGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAAc                                                                                                                             <  2:658774/149‑1 (MQ=255)
                   cgcTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTg                                                                                                                         >  1:722673/1‑149 (MQ=255)
                    gcTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTgg                                                                                                                        <  2:600837/149‑1 (MQ=255)
                     cTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGc                                                                                                                       <  2:530287/149‑1 (MQ=255)
                       ggTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCAc                                                                                                                     <  4:816142/149‑1 (MQ=255)
                        gTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACt                                                                                                                    <  1:707041/149‑1 (MQ=255)
                             ggCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTc                                                                                                               >  4:204801/1‑149 (MQ=255)
                             ggCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTc                                                                                                               >  2:865884/1‑149 (MQ=255)
                                ttttCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGt                                                                                                            >  3:85590/1‑149 (MQ=255)
                                 tttCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTg                                                                                                           <  3:420158/149‑1 (MQ=255)
                                   tCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGa                                                                                                         <  4:757915/149‑1 (MQ=255)
                                     aGACACAGGGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGCCGGTCTGGCTGGATCGCGAATCAGCTCTGTTTGGCGATTGAACATAGGCACTTCGTCAGTGGGTg                                                                                                       >  2:33584/1‑149 (MQ=255)
                                      gACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATgc                                                                                                      <  2:425821/149‑1 (MQ=255)
                                        cacaGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATgcgc                                                                                                    >  3:329555/1‑149 (MQ=255)
                                        cacaGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATgcgc                                                                                                    <  3:204801/149‑1 (MQ=255)
                                         acaGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGcgcg                                                                                                   >  1:270400/1‑149 (MQ=255)
                                             tgtgAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGaa                                                                                               >  3:83497/1‑149 (MQ=255)
                                             tgtgAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGaa                                                                                               >  3:321983/1‑149 (MQ=255)
                                              gtgAAGAGACACACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAACCTGTTCTGTTTGGCACTTCAACATTGGCACTTCGTCAGTGGATGCGCGAGaaa                                                                                              <  2:411304/149‑1 (MQ=255)
                                                gAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAgg                                                                                            >  1:351340/1‑149 (MQ=255)
                                                 aaGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGCGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGc                                                                                           >  1:532964/1‑149 (MQ=255)
                                                   gagaCTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACCTCGTCAGTGGATGCGCGAGAAAGGCAt                                                                                         >  4:328296/1‑149 (MQ=255)
                                                       cTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAg                                                                                     <  3:36439/149‑1 (MQ=255)
                                                        tCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGt                                                                                    >  4:921794/1‑149 (MQ=255)
                                                        tCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGt                                                                                    >  1:925883/1‑149 (MQ=255)
                                                         cACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGtt                                                                                   <  2:407359/149‑1 (MQ=255)
                                                          aCGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGttt                                                                                  <  2:422759/149‑1 (MQ=255)
                                                          aCGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGttt                                                                                  <  2:776769/149‑1 (MQ=255)
                                                                aGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATg                                                                            >  4:25660/1‑149 (MQ=255)
                                                                aGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATg                                                                            >  3:537869/1‑149 (MQ=255)
                                                                 ggCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGa                                                                           <  4:85590/149‑1 (MQ=255)
                                                                  gCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGaa                                                                          <  1:96336/149‑1 (MQ=255)
                                                                  gCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGaa                                                                          <  2:925883/149‑1 (MQ=255)
                                                                  gCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGaa                                                                          <  1:592669/149‑1 (MQ=255)
                                                                  gCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGaa                                                                          <  1:286686/149‑1 (MQ=255)
                                                                        cAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCaaaa                                                                    <  4:321983/149‑1 (MQ=255)
                                                                        cAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCaaaa                                                                    <  4:376442/149‑1 (MQ=255)
                                                                         aGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAAc                                                                   <  4:863928/149‑1 (MQ=255)
                                                                          gCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACg                                                                  >  3:259169/1‑149 (MQ=255)
                                                                            gtttgtttGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGaa                                                                >  1:105471/1‑149 (MQ=255)
                                                                            gtttgtttGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGaa                                                                >  2:96434/1‑149 (MQ=255)
                                                                            gtttgtttGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGaa                                                                <  3:921794/149‑1 (MQ=255)
                                                                             tttgtttgTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTg                                                                                                           <  3:30729/105‑1 (MQ=255)
                                                                             tttgtttgTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTg                                                                                                           >  4:30729/1‑105 (MQ=255)
                                                                             tttgtttgTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGGAGCAAAACGAAc                                                               <  1:865884/149‑1 (MQ=255)
                                                                             tttgtttgTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAAc                                                               >  2:350698/1‑149 (MQ=255)
                                                                             tttgtttgTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCACGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAAc                                                               >  4:139735/1‑149 (MQ=255)
                                                                                gtttgtCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAAt                                                                             >  4:827574/1‑132 (MQ=255)
                                                                                gtttgtCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAAt                                                                             <  3:827574/132‑1 (MQ=255)
                                                                                gtttgtCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCgg                                                            >  4:520527/1‑149 (MQ=255)
                                                                                     tCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGAt                                                       <  2:704053/149‑1 (MQ=255)
                                                                                     tCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTCGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGAt                                                       <  4:734787/149‑1 (MQ=255)
                                                                                       ggCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGaa                                                                          <  4:130110/128‑1 (MQ=255)
                                                                                       ggCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGaa                                                                          >  3:130110/1‑128 (MQ=255)
                                                                                        gCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGaa                                                          <  1:862239/143‑1 (MQ=255)
                                                                                        gCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGaa                                                          >  2:862239/1‑143 (MQ=255)
                                                                                        gCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCAt                                                    >  3:270324/1‑149 (MQ=255)
                                                                                         cccGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCAtt                                                   >  4:573677/1‑149 (MQ=255)
                                                                                          ccGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTg                                                                                                                         <  2:780536/78‑1 (MQ=255)
                                                                                          ccGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTg                                                                                                                         >  1:780536/1‑78 (MQ=255)
                                                                                          ccGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCAttt                                                  >  3:414953/1‑149 (MQ=255)
                                                                                          ccGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCAttt                                                  <  4:207870/149‑1 (MQ=255)
                                                                                           cGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCAtttt                                                 <  3:395557/149‑1 (MQ=255)
                                                                                            gTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTg                                                <  1:814930/149‑1 (MQ=255)
                                                                                            gTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTCTGGTGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTg                                                <  2:667087/149‑1 (MQ=255)
                                                                                              gggCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGt                                              >  3:834898/1‑149 (MQ=255)
                                                                                                  aCTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCAc                                          <  4:834898/149‑1 (MQ=255)
                                                                                                    tGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCActg                                        >  1:131383/1‑149 (MQ=255)
                                                                                                    tGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCActg                                        >  1:430431/1‑149 (MQ=255)
                                                                                                      ccGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCActgct                                      <  2:131383/149‑1 (MQ=255)
                                                                                                         tctcCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGtt                                   <  2:430431/149‑1 (MQ=255)
                                                                                                            ccgccgTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGAt                                <  3:25660/149‑1 (MQ=255)
                                                                                                              gccgTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATgg                              >  4:816383/1‑149 (MQ=255)
                                                                                                                    ggggTCGGTCTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGt                        >  4:366666/1‑149 (MQ=255)
                                                                                                                          ggtcTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTgg                  >  3:532083/1‑149 (MQ=255)
                                                                                                                           gtcTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGc                 <  3:573677/149‑1 (MQ=255)
                                                                                                                            tcTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGa                                 >  1:913081/1‑132 (MQ=255)
                                                                                                                            tcTGGCTGGATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGa                                 <  2:913081/132‑1 (MQ=255)
                                                                                                                                 ctggATCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGaa           <  4:543367/149‑1 (MQ=255)
                                                                                                                                   ggATCGCGAATCTGTGCTTTTTGGAGAGTCAACATTGGCACTCCGTCGGTGGATGCGCGAGAAAGGCATTCCGTTTCAAATGAAACAAAACGAACCGGAAGATCAGTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATg         >  4:374920/1‑149 (MQ=255)
                                                                                                                                      tCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGGTGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTc      <  4:532083/149‑1 (MQ=255)
                                                                                                                                      tCGCGAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTc      <  4:414953/149‑1 (MQ=255)
                                                                                                                                        gcgAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGc    >  3:752021/1‑149 (MQ=255)
                                                                                                                                        gcgAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGc    >  4:136970/1‑149 (MQ=255)
                                                                                                                                          gAATCTGTTCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCca  >  2:788660/1‑149 (MQ=255)
                                                                                                                                               tgttCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTgg                      <  3:195262/124‑1 (MQ=255)
                                                                                                                                               tgttCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTgg                      >  4:195262/1‑124 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCA  >  NC_000913/1664650‑1664936

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