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 ← peroxyureidoacrylate/ureidoacrylate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,072,3600TA 5.3% 482.9 / 9.8 149E168D (GAA→GATrutBperoxyureidoacrylate/ureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (61/80);  new base A (5/3);  total (66/83)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 151 total aligned reads displayed.

AGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCCAC  >  NC_000913/1072212‑1072509
                                                                                                                                                    |                                                                                                                                                     
aGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAtt                                                                                                                                                      >  1:819691/1‑150 (MQ=255)
 gAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTc                                                                                                                                                     <  1:2411328/150‑1 (MQ=255)
    tGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCAGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGAc                                                                                                                                                  <  2:2001226/150‑1 (MQ=255)
     gTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACg                                                                                                                                                 >  2:875051/1‑150 (MQ=255)
      tttCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGc                                                                                                                                                <  2:2472727/150‑1 (MQ=255)
       ttCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCa                                                                                                                                               <  2:1790664/150‑1 (MQ=255)
       ttCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCa                                                                                                                                               <  2:770179/150‑1 (MQ=255)
        tCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAg                                                                                                                                              <  2:2278075/150‑1 (MQ=255)
         cGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGa                                                                                                                                             <  2:2361594/150‑1 (MQ=255)
         cGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGa                                                                                                                                             <  1:1251948/150‑1 (MQ=255)
          gACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGAc                                                                                                                                            >  1:1461836/1‑150 (MQ=255)
           ccGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAATCGACGctacgc                                                                                                                                           <  1:118034/149‑6 (MQ=255)
             gTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGtt                                                                                                                                         <  1:1979958/150‑1 (MQ=255)
                gCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGt                                                                                                                                      >  2:1544232/1‑150 (MQ=255)
                  tGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAg                                                                                                                                    <  1:2354415/150‑1 (MQ=255)
                   gACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGc                                                                                                                                   <  1:139684/150‑1 (MQ=255)
                   gACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGAATCGACGctacgc                                                                                                                                           >  2:1436375/1‑137 (MQ=255)
                    aCCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCg                                                                                                                                  <  1:1863163/150‑1 (MQ=255)
                      ccAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGAt                                                                                                                                >  2:1403353/1‑150 (MQ=255)
                      ccAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGAt                                                                                                                                <  1:1162364/150‑1 (MQ=255)
                           cAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCgg                                                                                                                           <  2:252067/150‑1 (MQ=255)
                           cAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCgg                                                                                                                           >  2:1046448/1‑150 (MQ=255)
                               aaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaa                                                                                                                       <  1:1107067/150‑1 (MQ=255)
                                 aGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGaaaa                                                                                                                     <  2:642579/150‑1 (MQ=255)
                                  ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                    >  2:1350628/1‑150 (MQ=255)
                                  ggTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAc                                                                                                                    >  1:1396185/1‑150 (MQ=255)
                                   gTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAcc                                                                                                                   <  2:427473/150‑1 (MQ=255)
                                   gTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAAcc                                                                                                                   <  1:1123667/150‑1 (MQ=255)
                                           atTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                                                                                           <  1:370327/150‑1 (MQ=255)
                                           atTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCg                                                                                                           >  2:2350459/1‑150 (MQ=255)
                                            ttGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGt                                                                                                          <  2:2134422/150‑1 (MQ=255)
                                            ttGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGt                                                                                                          <  1:1329662/150‑1 (MQ=255)
                                               aacaacGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTAtt                                                                                                       >  1:1491688/1‑150 (MQ=255)
                                                      cgGCTTTCTGCGCAAATTTCGGCCGCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGc                                                                                                <  2:446709/150‑1 (MQ=255)
                                                      cgGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTTGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGc                                                                                                <  1:1815716/150‑1 (MQ=255)
                                                            tCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCa                                                                                          <  2:1273160/150‑1 (MQ=255)
                                                                cgcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAAt                                                                                      >  1:42562/1‑150 (MQ=255)
                                                                 gcAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATg                                                                                     >  1:268159/1‑150 (MQ=255)
                                                                  cAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGc                                                                                    <  1:509133/150‑1 (MQ=255)
                                                                  cAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGc                                                                                    >  1:1773504/1‑150 (MQ=255)
                                                                    aaTTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTg                                                                                  <  2:538325/150‑1 (MQ=255)
                                                                      tttCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTc                                                                                <  1:1318169/150‑1 (MQ=255)
                                                                          ggCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACACCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgc                                                                            >  2:346669/1‑150 (MQ=255)
                                                                           gCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcg                                                                           <  1:1815250/150‑1 (MQ=255)
                                                                           gCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcg                                                                           >  1:623225/1‑150 (MQ=255)
                                                                            ccccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcgg                                                                          >  2:238288/1‑150 (MQ=255)
                                                                            ccccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcgg                                                                          <  2:1461836/150‑1 (MQ=255)
                                                                              ccGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAgcggcg                                                                        <  1:995426/150‑1 (MQ=255)
                                                                                 ccTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTAt                                                                     <  1:1937662/150‑1 (MQ=255)
                                                                                   tGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGCGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTg                                                                   <  2:1655781/150‑1 (MQ=255)
                                                                                    ggTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTACAGAAAAAAGCCGTCGACTAGCGTCGATTCGACGCAGACGTTCGTAGAGATGCCCGTGAAAACCAGATCGCGTATACCGCGGCTGCGCAAAATGCTGTACAGCGGCGTATTga                                                                  >  2:17522/1‑150 (MQ=255)
                                                                                      tGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaag                                                                <  1:279611/150‑1 (MQ=255)
                                                                                      tGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaag                                                                <  1:1643947/150‑1 (MQ=255)
                                                                                      tGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaag                                                                <  1:1544232/150‑1 (MQ=255)
                                                                                       gAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTgaaga                                                               <  2:358415/150‑1 (MQ=255)
                                                                                         gTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGaaa                                                             >  1:79366/1‑150 (MQ=255)
                                                                                           tGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAAcc                                                           <  2:798951/150‑1 (MQ=255)
                                                                                             cGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGc                                                         <  2:2015128/150‑1 (MQ=255)
                                                                                               tCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTg                                                       >  1:2422796/1‑150 (MQ=255)
                                                                                               tCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTg                                                       >  1:1045203/1‑150 (MQ=255)
                                                                                                cTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGt                                                      >  2:1919994/1‑150 (MQ=255)
                                                                                                 ttCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTa                                                     <  1:2158047/150‑1 (MQ=255)
                                                                                                  tCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAg                                                    >  2:1231646/1‑150 (MQ=255)
                                                                                                   cAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAgc                                                   >  2:194301/1‑150 (MQ=255)
                                                                                                    aaGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAgcg                                                  <  1:1046448/150‑1 (MQ=255)
                                                                                                        accacGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGAGGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGAGCAAAATGCTGTCGAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGc                                              <  2:1859788/150‑1 (MQ=255)
                                                                                                         ccacGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCt                                             >  2:1090193/1‑150 (MQ=255)
                                                                                                          cacGCCGAAATACTCCAGAAAAAAGCCGGCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCTCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTATCGCGTCtt                                            <  2:1172348/150‑1 (MQ=255)
                                                                                                            cGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCg                                          >  1:2087791/1‑150 (MQ=255)
                                                                                                              ccGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgc                                        <  2:1785618/150‑1 (MQ=255)
                                                                                                               cGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgcc                                       >  1:347826/1‑149 (MQ=255)
                                                                                                                gAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgcag                                      >  1:153157/1‑150 (MQ=255)
                                                                                                                gAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgcag                                      <  1:289729/150‑1 (MQ=255)
                                                                                                                 aaaTACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgcagc                                     <  1:1202580/150‑1 (MQ=255)
                                                                                                                 aaaTACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGgcagc                                     <  2:608804/150‑1 (MQ=255)
                                                                                                                   aTACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGcac                                   >  1:1890655/1‑150 (MQ=255)
                                                                                                                    tACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTg                                                                                  >  1:2254424/1‑102 (MQ=255)
                                                                                                                      cTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAat                                >  2:765052/1‑150 (MQ=255)
                                                                                                                      cTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTAAGAAACCCCGGTAGCGCGGCGTACGCAGCACCata                               >  1:535566/1‑150 (MQ=255)
                                                                                                                        ccAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAatat                              <  2:693231/150‑1 (MQ=255)
                                                                                                                        ccAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAatat                              >  1:652046/1‑150 (MQ=255)
                                                                                                                           gAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGc                           <  2:819691/150‑1 (MQ=255)
                                                                                                                             aaaaaGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCa                         <  2:623225/150‑1 (MQ=255)
                                                                                                                                  gCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTg                    >  1:2205459/1‑150 (MQ=255)
                                                                                                                                      gcgcgTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGGGTTTTGAAGACACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCtgctgc                <  1:1858133/149‑1 (MQ=255)
                                                                                                                                         cgTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAgg                       <  1:953645/140‑1 (MQ=255)
                                                                                                                                         cgTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCAcc             <  1:121029/150‑1 (MQ=255)
                                                                                                                                          tctgCGTCGAATCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCa            >  2:1406669/4‑150 (MQ=255)
                                                                                                                                          tctgCGTCGAATCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCa            >  1:558490/4‑150 (MQ=255)
                                                                                                                                           tAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAg           <  1:875051/150‑1 (MQ=255)
                                                                                                                                           tAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAg           >  2:204023/1‑150 (MQ=255)
                                                                                                                                           ctgCGTCGAATCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAg           >  1:1469975/3‑150 (MQ=255)
                                                                                                                                            tgCGTCGAATCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGt          >  2:721649/2‑150 (MQ=255)
                                                                                                                                             gCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGtt         >  2:2042029/1‑150 (MQ=255)
                                                                                                                                                tCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCAt      <  2:1881045/150‑1 (MQ=255)
                                                                                                                                                gcGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCAt      <  2:1602011/149‑1 (MQ=255)
                                                                                                                                                  gATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATcc    >  1:1490380/1‑150 (MQ=255)
                                                                                                                                                  gATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATcc    >  1:478154/1‑150 (MQ=255)
                                                                                                                                                  gATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATcc    >  2:2194897/1‑150 (MQ=255)
                                                                                                                                                    ttCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATccac  <  1:1251955/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGAATGTTTCGACGTCGCTGACCCAGCCAAAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCCAC  >  NC_000913/1072212‑1072509

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: