Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,072,613 G→T 5.6% P84H (CCC→CAC)  rutB ← ureidoacrylate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,072,6130GT 5.6% 280.9 / 3.7 88P84H (CCC→CAC) rutBureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base T (3/2);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCT  >  NC_000913/1072471‑1072756
                                                                                                                                              |                                                                                                                                               
cagcaCAATATCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGc                                                                                                                                           >  6:107068/1‑149 (MQ=255)
    acaATATCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGc                                                                                                                                       >  2:131870/1‑149 (MQ=255)
    acaATATCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGc                                                                                                                                       >  1:278950/1‑149 (MQ=255)
       atatCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTc                                                                                                                                    <  1:326026/149‑1 (MQ=255)
          tCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCa                                                                                                                                         >  5:225967/1‑141 (MQ=255)
          tCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCa                                                                                                                                         <  6:225967/141‑1 (MQ=255)
            gCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCGGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACAt                                                                                                                               <  5:107068/149‑1 (MQ=255)
             ccAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATa                                                                                                                              <  7:304002/149‑1 (MQ=255)
             ccAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATa                                                                                                                              <  7:157494/149‑1 (MQ=255)
               aGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACt                                                                                                                            >  5:324275/1‑149 (MQ=255)
                  cTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGtt                                                                                                                         <  2:357717/149‑1 (MQ=255)
                     cGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCa                                                                                                                       <  7:307337/148‑1 (MQ=255)
                       gCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGa                                                                                                                                  <  5:315240/135‑1 (MQ=255)
                       gCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGa                                                                                                                                  >  6:315240/1‑135 (MQ=255)
                             gTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCa                                                                                                              >  6:245136/1‑149 (MQ=255)
                                    caccaGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGaa                                                                                                       <  2:264634/149‑1 (MQ=255)
                                     accaGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGaaa                                                                                                      >  3:228338/1‑149 (MQ=255)
                                     accaGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGGTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACGGTTCATCCCAGCCATTTTGaac                                                                                                      >  3:179858/1‑148 (MQ=255)
                                     accaGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGAGGCTGCTTACGCATGGTTTTCAGCGCGTTCGATTTATGAAAATTCGGTGAGCCGGGACCGCCAGCCTCGACATACTGTTCATCCACCCATTTTGaacc                                                                                                     >  8:344463/1‑147 (MQ=255)
                                         gTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAg                                                                                                  <  1:131870/149‑1 (MQ=255)
                                          ttGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCGga                                                                                                 >  3:333521/1‑149 (MQ=255)
                                            gATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATAc                                                                                                                             >  7:338939/1‑119 (MQ=255)
                                            gATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATAc                                                                                                                             <  8:338939/119‑1 (MQ=255)
                                             aTAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAgatga                                                                                              >  3:159687/1‑149 (MQ=255)
                                                aTCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATc                                                                                            <  5:322467/148‑1 (MQ=255)
                                                 tCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAg                                                                                          >  1:312875/1‑149 (MQ=255)
                                                 tCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAg                                                                                          >  8:282677/1‑149 (MQ=255)
                                                   ccAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGc                                                                                         <  5:301226/148‑1 (MQ=255)
                                                    cAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCAt                                                                                       >  6:291317/1‑149 (MQ=255)
                                                    cAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCAt                                                                                       <  3:296793/149‑1 (MQ=255)
                                                    cAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCAt                                                                                       <  7:251529/149‑1 (MQ=255)
                                                     aGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATc                                                                                      <  2:340332/149‑1 (MQ=255)
                                                          ccTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGc                                                                                 >  5:219881/1‑149 (MQ=255)
                                                                gCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATa                                                                                                                              <  4:36101/98‑1 (MQ=255)
                                                                  cAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTcgcg                                                                         >  7:206826/1‑149 (MQ=255)
                                                                  cAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTcgc                                                                          >  1:4441/1‑148 (MQ=255)
                                                                            cccTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGTGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAc                                                               >  2:201797/1‑149 (MQ=255)
                                                                            cccTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAc                                                               >  1:28071/1‑149 (MQ=255)
                                                                            cccTGCAGCTGCGGCTGATTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGTGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAc                                                               >  4:78968/1‑149 (MQ=255)
                                                                              cTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTTCCGCTCGCGCTGCGGTCAcgg                                                             <  4:228338/149‑1 (MQ=255)
                                                                                    cTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGttt                                                       >  5:327503/1‑149 (MQ=255)
                                                                                    cTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGttt                                                       <  5:134889/149‑1 (MQ=255)
                                                                                         gctgcgTACGCATGGTGTTCAGGGCGGTCGATTTCTGAAAATTCGGGGAGCCGGGTCCGCCAGCCTCGACATACTGTGCATCCCAGCCATTTTGAAACCAGATCATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATg                                                  <  2:278140/149‑1 (MQ=255)
                                                                                          ctgctTACGCCTGGTTTTCAGGGCGTTCGATTTATGAAAATACGGGGACCCGGGTCCGCCAGCCTCGCCATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGt                                                 <  6:219881/149‑1 (MQ=255)
                                                                                            gctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAcggcg                                                           <  4:239012/137‑1 (MQ=255)
                                                                                            gctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAcggcg                                                           >  3:239012/1‑137 (MQ=255)
                                                                                            gctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTg                                               >  5:280639/1‑149 (MQ=255)
                                                                                             ctTACGCCTGGTTATCAGGGCGTTCGATTTCTGAAACTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGtt                                                <  6:238864/147‑1 (MQ=255)
                                                                                             ctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTgg                                              >  4:339750/1‑149 (MQ=255)
                                                                                             ctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTgg                                              <  5:169571/149‑1 (MQ=255)
                                                                                             ctTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCACTGCGGTCACGGCGGTTTGAATGTTgg                                              <  4:159687/149‑1 (MQ=255)
                                                                                              ttACGCATGGTTTTCAGGTCGTTAGATTTATGAAAATTCGGTGAGCCGTGTCAGACAGCCTAGACATAATGTTAATCAAAGCCATTTTGAAACCAGATGATCATCATACCTGCCGCTCGCTCTGCGCTCACGGCGCTTTGAATGTTgg                                              >  4:113481/1‑148 (MQ=255)
                                                                                                  gCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATg                                         >  4:198412/1‑149 (MQ=255)
                                                                                                     tGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGAcc                                      >  1:243746/1‑149 (MQ=255)
                                                                                                       gTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATcc                                                                                     >  1:351547/1‑100 (MQ=255)
                                                                                                       gTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATcc                                                                                     <  2:351547/100‑1 (MQ=255)
                                                                                                       gTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCgg                                    >  3:121114/1‑149 (MQ=255)
                                                                                                        ttttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCggg                                   >  2:138733/1‑149 (MQ=255)
                                                                                                        ttttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCggg                                   >  6:82494/1‑149 (MQ=255)
                                                                                                          ttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGc                                  <  2:271927/148‑1 (MQ=255)
                                                                                                            cAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGa                                <  2:4441/148‑1 (MQ=255)
                                                                                                            cAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGa                                <  8:206826/148‑1 (MQ=255)
                                                                                                              gggCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTg                             >  3:151390/1‑149 (MQ=255)
                                                                                                               ggCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTgg                            >  1:333549/1‑149 (MQ=255)
                                                                                                                  gTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGtt                          >  2:241032/1‑148 (MQ=255)
                                                                                                                  gTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTg                         <  5:291317/149‑1 (MQ=255)
                                                                                                                   ttCGATTTATTAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTga                        <  3:113481/149‑1 (MQ=255)
                                                                                                                   ttCGATTTATGAAAATTCGGTGAGCCGTGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTga                        <  1:201797/149‑1 (MQ=255)
                                                                                                                   ttCGATTTATGAAAATTCGGTGAGCCGTGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTga                        <  3:78968/149‑1 (MQ=255)
                                                                                                                   ttCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTga                        <  1:138733/149‑1 (MQ=255)
                                                                                                                   ttCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTga                        <  1:241032/149‑1 (MQ=255)
                                                                                                                      gATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTgaga                      >  6:51349/1‑148 (MQ=255)
                                                                                                                       aTTTATGAAAATTCGGTGAGCCGGGTCCGCCCGCCTCGACCTACGGTTCATCCCAGCCATTTTGACACCCGATGATCCGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGAcc                    >  5:263600/1‑148 (MQ=255)
                                                                                                                         ttATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGa                                <  7:210963/135‑1 (MQ=255)
                                                                                                                         ttATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGa                                >  8:210963/1‑135 (MQ=255)
                                                                                                                             gAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGTGGTTTGAATGTTGGCAATGACCGGGCGa                                >  1:35614/1‑131 (MQ=255)
                                                                                                                             gAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGTGGTTTGAATGTTGGCAATGACCGGGCGa                                <  2:35614/131‑1 (MQ=255)
                                                                                                                             gAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTgg                                              <  7:20766/117‑1 (MQ=255)
                                                                                                                             gAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTgg                                              >  8:20766/1‑117 (MQ=255)
                                                                                                                              aaaaTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAccc             >  5:235667/1‑149 (MQ=255)
                                                                                                                               aaaTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAAccc            <  2:312875/149‑1 (MQ=255)
                                                                                                                                 aTTCGGTGAGCCGGGTCCGACAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCgg          <  2:243746/149‑1 (MQ=255)
                                                                                                                                    cGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACAt                   >  7:301681/1‑137 (MQ=255)
                                                                                                                                    cGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACAt                   <  8:301681/137‑1 (MQ=255)
                                                                                                                                      gTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCgaga     <  5:245136/149‑1 (MQ=255)
                                                                                                                                       tGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGc         >  7:63978/1‑144 (MQ=255)
                                                                                                                                       tGAGCCGGGGCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGAt    <  6:324275/149‑1 (MQ=255)
                                                                                                                                           ccGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCt  >  6:201038/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCT  >  NC_000913/1072471‑1072756

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