Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,072,360 T→A 5.3% E168D (GAA→GAT)  rutB ← ureidoacrylate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,072,3600TA 5.3% 340.3 / 6.7 94E168D (GAA→GAT) rutBureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (44/45);  new base A (2/3);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCC  >  NC_000913/1072213‑1072507
                                                                                                                                                   |                                                                                                                                                   
gAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAtt                                                                                                                                                    >  1:699171/1‑149 (MQ=255)
        cGACGTCGCTCATCCAGCCAAAAAAGGTTTCGAGATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAATCGACGcta                                                                                                                                            <  3:1224176/149‑3 (MQ=255)
             tCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCg                                                                                                                                                  >  2:139222/1‑138 (MQ=255)
             tCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCg                                                                                                                                                  <  1:139222/138‑1 (MQ=255)
                     ccAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGa                                                                                                                               >  1:912608/1‑149 (MQ=255)
                     ccAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGa                                                                                                                               >  4:968384/1‑149 (MQ=255)
                     ccAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGa                                                                                                                               >  4:23455/1‑149 (MQ=255)
                     ccAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGa                                                                                                                               <  4:860168/149‑1 (MQ=255)
                         ccAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGcc                                                                                                                           >  4:592864/1‑149 (MQ=255)
                           aaaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGa                                                                                                                                           <  3:87997/131‑1 (MQ=255)
                           aaaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGa                                                                                                                                           >  4:87997/1‑131 (MQ=255)
                            aaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTc                                                                                                                                                   >  4:1313506/1‑122 (MQ=255)
                            aaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTc                                                                                                                                                   <  3:1313506/122‑1 (MQ=255)
                                  gTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                  <  2:699171/149‑1 (MQ=255)
                                  gTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                  <  3:23455/149‑1 (MQ=255)
                                   tttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGc                                                                                                                            <  4:30031/138‑1 (MQ=255)
                                   tttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGc                                                                                                                            >  3:30031/1‑138 (MQ=255)
                                   tttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAATCGACGcta                                                                                                                                            <  1:1225888/122‑3 (MQ=255)
                                   tttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAATCGACGcta                                                                                                                                            >  2:1225888/1‑120 (MQ=255)
                                    ttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTg                                                                                                                                      <  3:1450327/127‑1 (MQ=255)
                                    ttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTg                                                                                                                                      >  4:1450327/1‑127 (MQ=255)
                                       gATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAg                                                                                                                                            >  3:988497/1‑118 (MQ=255)
                                       gATATCGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAg                                                                                                                                            <  4:988497/118‑1 (MQ=255)
                                        atatTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATg                                                                                                            >  3:772277/1‑149 (MQ=255)
                                           ttGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                                                                                         >  3:37325/1‑149 (MQ=255)
                                           ttGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                                                                                         >  2:872116/1‑149 (MQ=255)
                                           ttGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                                                                                         >  4:1459633/1‑149 (MQ=255)
                                           ttGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                                                                                         >  1:462253/1‑149 (MQ=255)
                                             gAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                                                                                       <  3:1459633/149‑1 (MQ=255)
                                             gAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                                                                                       <  1:872116/149‑1 (MQ=255)
                                             gAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                                                                                       <  4:772277/149‑1 (MQ=255)
                                                caacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTc                                                                                                    <  3:968384/149‑1 (MQ=255)
                                                       gCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCt                                                                                             >  1:391309/1‑149 (MQ=255)
                                                        cTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTg                                                                                            >  3:1130562/1‑149 (MQ=255)
                                                                      ttCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTc                                                                              <  2:462253/149‑1 (MQ=255)
                                                                          gCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGt                                                                                                                                        >  3:1269049/1‑87 (MQ=255)
                                                                          gCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGt                                                                                                                                        <  4:1269049/87‑1 (MQ=255)
                                                                            cccGGCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcgg                                                                        <  4:37325/149‑1 (MQ=255)
                                                                            cccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATg                                                                                   <  3:463618/138‑1 (MQ=255)
                                                                            cccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATg                                                                                   >  4:463618/1‑138 (MQ=255)
                                                                            cccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcgg                                                                        >  4:1364163/1‑149 (MQ=255)
                                                                                   ggTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTg                                                                 <  3:1364163/149‑1 (MQ=255)
                                                                                         ttGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcggcg                                                                      >  2:274054/1‑138 (MQ=255)
                                                                                         ttGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcggcg                                                                      <  1:274054/138‑1 (MQ=255)
                                                                                         ttGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGaaa                                                           <  2:912608/149‑1 (MQ=255)
                                                                                           gCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAAcc                                                         >  2:322588/1‑149 (MQ=255)
                                                                                                    aGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                  >  4:1437037/1‑83 (MQ=255)
                                                                                                    aGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                  <  3:53704/83‑1 (MQ=255)
                                                                                                    aGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                  <  3:1437037/83‑1 (MQ=255)
                                                                                                    aGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                  >  4:53704/1‑83 (MQ=255)
                                                                                                     gCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAgcgc                                               <  4:1130562/149‑1 (MQ=255)
                                                                                                      caccacGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGcgcg                                              <  3:592864/149‑1 (MQ=255)
                                                                                                        ccacGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGc                                            <  2:391309/149‑1 (MQ=255)
                                                                                                         cacGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCt                                           >  4:280176/1‑149 (MQ=255)
                                                                                                                  aTACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGt                                                                               >  3:846180/1‑104 (MQ=255)
                                                                                                                  aTACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGt                                                                               <  4:846180/104‑1 (MQ=255)
                                                                                                                      tCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               >  2:573510/1‑116 (MQ=255)
                                                                                                                      tCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               <  1:573510/116‑1 (MQ=255)
                                                                                                                         aGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                                                                                       >  1:428056/1‑73 (MQ=255)
                                                                                                                         aGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTa                                                                                                       <  2:428056/73‑1 (MQ=255)
                                                                                                                         aGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATc                           >  1:763893/1‑149 (MQ=255)
                                                                                                                         aGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATc                           <  3:280176/149‑1 (MQ=255)
                                                                                                                         aGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATc                           >  3:604582/1‑149 (MQ=255)
                                                                                                                          gAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               >  1:354270/1‑112 (MQ=255)
                                                                                                                          gAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               <  2:354270/112‑1 (MQ=255)
                                                                                                                          gAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgca                                     >  3:1231237/1‑138 (MQ=255)
                                                                                                                          gAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgca                                     <  4:1231237/138‑1 (MQ=255)
                                                                                                                          gAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCg                          >  3:1433837/1‑149 (MQ=255)
                                                                                                                           aaaaaaGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               >  1:181256/1‑111 (MQ=255)
                                                                                                                           aaaaaaGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               <  2:181256/111‑1 (MQ=255)
                                                                                                                             aaaaGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCa                       <  4:1433837/149‑1 (MQ=255)
                                                                                                                              aaaGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAg                      <  2:763893/149‑1 (MQ=255)
                                                                                                                                aGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGc                    <  4:604582/149‑1 (MQ=255)
                                                                                                                                    gTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAgcgc                                               <  1:332055/118‑1 (MQ=255)
                                                                                                                                    gTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAgcgc                                               >  2:332055/1‑118 (MQ=255)
                                                                                                                                     tCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcg                                                                         >  2:1020772/1‑91 (MQ=255)
                                                                                                                                     tCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcg                                                                         <  1:1020772/91‑1 (MQ=255)
                                                                                                                                       gcgTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATg                                                                                   >  1:951232/1‑79 (MQ=255)
                                                                                                                                       gcgTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATg                                                                                   <  2:951232/79‑1 (MQ=255)
                                                                                                                                       gcgTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCa             >  1:1194714/1‑149 (MQ=255)
                                                                                                                                          ctgCGTCGAATCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGc                 <  3:994370/140‑1 (MQ=255)
                                                                                                                                          ctgCGTCGAATCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGc                 >  4:994370/3‑142 (MQ=255)
                                                                                                                                            gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAGAACAAGATGGc                                                                                                          >  3:905247/1‑51 (MQ=255)
                                                                                                                                            gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGc                                                                                                          <  4:905247/51‑1 (MQ=255)
                                                                                                                                            gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               <  3:502494/94‑1 (MQ=255)
                                                                                                                                            gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaa                                                               >  4:502494/1‑94 (MQ=255)
                                                                                                                                            gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCa          >  2:1441230/1‑147 (MQ=255)
                                                                                                                                            gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCa          <  1:1441230/147‑1 (MQ=255)
                                                                                                                                            gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGt        <  1:322588/149‑1 (MQ=255)
                                                                                                                                               tCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTgc                                                                                           >  2:259622/1‑63 (MQ=255)
                                                                                                                                               tCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTgc                                                                                           <  1:259622/63‑1 (MQ=255)
                                                                                                                                               tCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCa     <  2:1194714/149‑1 (MQ=255)
                                                                                                                                                  aTTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTc                                         >  2:191808/1‑110 (MQ=255)
                                                                                                                                                  aTTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTc                                         <  1:191808/110‑1 (MQ=255)
                                                                                                                                                  aTTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATcc  >  1:832472/1‑149 (MQ=255)
                                                                                                                                                   ttCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTc                                         <  3:1168289/109‑1 (MQ=255)
                                                                                                                                                   ttCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTc                                         >  4:1168289/1‑109 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCC  >  NC_000913/1072213‑1072507

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