Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,682 T→A 6.8% V381V (GTT→GTA) ‡ waaA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,6820TA 6.8% 272.4 / 7.1 88V381V (GTT→GTA) ‡waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base T (41/40);  new base A (2/4);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGC  >  NC_000913/3809536‑3809828
                                                                                                                                                  |                                                                                                                                                  
gggCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTc                                                                                                                                                  >  4:60077/1‑149 (MQ=255)
 ggCATAATCCGCTGGATGCTGCCGCACAGGCTATTCCGGTATTTATGGGGCCGCATTCTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTc                                                                                                                                                  <  2:256241/148‑1 (MQ=255)
  gCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCct                                                                                                                                                >  4:369963/1‑149 (MQ=255)
    aTAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCctct                                                                                                                                              <  8:288185/149‑1 (MQ=255)
       aTCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                           <  5:99744/149‑1 (MQ=255)
         tcGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                           <  3:2741/146‑1 (MQ=255)
          cGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgc                                                                                                                                         <  8:222839/148‑4 (MQ=255)
           gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACCTTTGCGCGAGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAATACGCTTGCAAAAGAGGTTTCCTCTTTACTCa                                                                                                                                       <  1:387507/149‑1 (MQ=255)
           gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTc                                                                                                                                                  <  8:48100/138‑1 (MQ=255)
           gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTc                                                                                                                                                  >  7:48100/1‑138 (MQ=255)
           gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTa                                                                                                                                           >  4:2741/1‑145 (MQ=255)
           gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgcaa                                                                                                                                       <  5:396161/149‑6 (MQ=255)
           gCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgcaa                                                                                                                                       <  7:43582/149‑6 (MQ=255)
           gCTGGAAGCTGCCGCACACGCGATTCCGGTATGGATGGGGCCGCATACTTTTACCTTTAAAGACATTTGCGCGCGGCGGGAGCAGGCAAGCGGGCTTATTACCGTTACCGCTGCCCCTCCGCTTGCAAAAGAGGTGTCCTCTTTACTCa                                                                                                                                       <  1:393174/149‑1 (MQ=255)
            cTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCAc                                                                                                                                      <  2:240855/149‑1 (MQ=255)
            cTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCAc                                                                                                                                      <  3:195838/149‑1 (MQ=255)
                aaGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGAc                                                                                                                                  >  4:222953/1‑149 (MQ=255)
                  gCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGc                                                                                                                                <  3:61115/149‑1 (MQ=255)
                      ccGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGAt                                                                                                                            <  1:39945/149‑1 (MQ=255)
                        gCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCct                                                                                                                                                >  6:345241/1‑127 (MQ=255)
                        gCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCct                                                                                                                                                <  5:345241/127‑1 (MQ=255)
                        gCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTa                                                                                                                          <  5:377779/149‑1 (MQ=255)
                         cacaCGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTAc                                                                                                                         <  3:222953/149‑1 (MQ=255)
                          acacGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTAcc                                                                                                                        >  8:225998/1‑149 (MQ=255)
                            acGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGa                                                                                                                                   <  4:250484/136‑1 (MQ=255)
                            acGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGa                                                                                                                                   >  3:250484/1‑136 (MQ=255)
                              gcaaTTCCGGTAATGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGAGCGGCTGGAGCAGGCAACCGGGCAGATTACCCTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAg                                                                                                                    >  4:270696/4‑149 (MQ=255)
                                tATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTCCTCACCGACGCCGATTACCGTAGtt                                                                                                                  >  1:65293/1‑149 (MQ=255)
                                  ttCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgcaag                                                                                                                                      >  4:220223/1‑121 (MQ=255)
                                  ttCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGAAACCTCttttgcaag                                                                                                                                      <  3:220223/127‑7 (MQ=255)
                                     cGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa                                                                                                              <  3:60077/148‑1 (MQ=255)
                                     cGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTAt                                                                                                             >  5:346862/1‑149 (MQ=255)
                                      ggTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCTAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATg                                                                                                            <  3:126646/149‑1 (MQ=255)
                                            gATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGt                                                                                                       <  6:349613/148‑1 (MQ=255)
                                             aTGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCa                                                                                                     >  5:59976/1‑149 (MQ=255)
                                               ggggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATg                                                                                                   >  2:318574/1‑149 (MQ=255)
                                                 ggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGt                                                                                                                   <  3:253750/131‑1 (MQ=255)
                                                 ggCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGt                                                                                                                   >  4:253750/1‑131 (MQ=255)
                                                    cGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc                                                                                                                >  7:3097/1‑131 (MQ=255)
                                                    cGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTc                                                                                                                <  8:3097/131‑1 (MQ=255)
                                                     gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGTTGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGtt                                                                                              >  4:152743/1‑148 (MQ=255)
                                                     gCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGACACTACGCTTGCAAAAGATGTTTCCTCTTTACACACAGACGCCTATTACCGTAGTTTCTATGGCCGTCATGCCGtt                                                                                              >  5:163231/1‑148 (MQ=255)
                                                        tACTTTTAACTTTAAAGTCATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAg                                                                                          <  2:175472/149‑1 (MQ=255)
                                                         actataaaCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGAAGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTAACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAg                                                                                          >  2:74380/7‑148 (MQ=255)
                                                           ttttAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTAc                                                                                       >  1:350338/1‑149 (MQ=255)
                                                               aaCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTAcc                                                                                                                        >  4:34110/1‑112 (MQ=255)
                                                               aaCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGt                                                                                               >  6:234503/1‑137 (MQ=255)
                                                               aaCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGt                                                                                               <  5:234503/137‑1 (MQ=255)
                                                               aaCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGt                                                                                    <  1:318574/148‑1 (MQ=255)
                                                               aaCTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTa                                                                                   >  1:340792/1‑149 (MQ=255)
                                                                   ttAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCa                                                                                >  1:46196/1‑148 (MQ=255)
                                                                          cATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAgg                                                                        <  3:369963/149‑1 (MQ=255)
                                                                           aTTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAggg                                                                       >  3:15480/1‑149 (MQ=255)
                                                                                cgcgcgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGgcgcg                                                                   <  7:225998/148‑1 (MQ=255)
                                                                                 gcgcgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGgcgcgc                                                                  >  5:212533/1‑148 (MQ=255)
                                                                                 gcgcgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGgcgcgc                                                                  >  5:224456/1‑148 (MQ=255)
                                                                                 gcgcgGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGgcgcgc                                                                  <  6:212533/148‑1 (MQ=255)
                                                                                             cAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCt                                                     >  8:394930/1‑149 (MQ=255)
                                                                                              aggacaGCGTGCTTATTACCGTTAACGATGCCACTACGCTTGCAATAGAGGTTTCCTCTTTACTCACCGACCCCGATTACCGTAGTTTCTATGGCCGTCATTCAGTTTAAGTAGTTTATCATAACCAGGGCGCGCTACAGCGTCTGCtt                                                    <  5:34164/144‑1 (MQ=255)
                                                                                               ggCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCtt                                                    >  4:202706/1‑148 (MQ=255)
                                                                                                  aaGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAAc                                                >  8:382565/1‑149 (MQ=255)
                                                                                                     cGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActgc                                             <  4:15480/149‑1 (MQ=255)
                                                                                                     cGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActgc                                             <  3:202706/149‑1 (MQ=255)
                                                                                                      gggCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActgct                                            <  6:346862/149‑1 (MQ=255)
                                                                                                       ggCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAActgctg                                           <  3:152743/149‑1 (MQ=255)
                                                                                                        gCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGGCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTgg                                          <  6:59976/149‑1 (MQ=255)
                                                                                                            aTTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAAc                                       <  2:341008/148‑1 (MQ=255)
                                                                                                            aTTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAAc                                       <  2:350338/148‑1 (MQ=255)
                                                                                                               accgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCtt                                    <  3:270696/148‑1 (MQ=255)
                                                                                                                 cgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTAcc                                 >  6:20832/1‑149 (MQ=255)
                                                                                                                 cgttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTAcc                                 >  4:175563/1‑149 (MQ=255)
                                                                                                                   ttaccgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTg                               >  2:223431/1‑149 (MQ=255)
                                                                                                                     accgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGc                              <  3:204571/148‑1 (MQ=255)
                                                                                                                       cgATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCAc                           <  2:340792/149‑1 (MQ=255)
                                                                                                                             cACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCACCTGCTGGACCCTTACCTGCCACCGaaac                     >  6:10522/1‑148 (MQ=255)
                                                                                                                              aCTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAAc                    >  8:149460/1‑149 (MQ=255)
                                                                                                                              aCTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAAc                    >  4:140336/1‑149 (MQ=255)
                                                                                                                                     ttGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGaa                                        >  6:328608/1‑122 (MQ=255)
                                                                                                                                     ttGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGaa                                        <  5:328608/122‑1 (MQ=255)
                                                                                                                                      gagtaaaGAGGAAACCTCTTTACTCACCGCCGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAg            >  2:357013/5‑149 (MQ=255)
                                                                                                                                         aaaaGAGGTTTCCTCTTTACTCCCCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGATCCTTACCTGCCACCGAAAACGCATTGAGgtt         >  3:90142/1‑149 (MQ=255)
                                                                                                                                         aaaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgtt         >  6:385538/1‑149 (MQ=255)
                                                                                                                                          aaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttg        >  1:209391/1‑149 (MQ=255)
                                                                                                                                          aaaGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGgttg        >  3:187991/1‑149 (MQ=255)
                                                                                                                                             tagGTTTCCTGTTTACTCACCGTCTCCGATTACCGTAGTTTCTATGGCGGTCATGCCGTTGAAGTACTGTATCAAAAGCAGGGCTCGCTACATCGTCTGCTTCAACTGCTTGAACCTTACCTGCCACCGAAATCGCATTGAGGTTGTTa     <  6:163231/148‑1 (MQ=255)
                                                                                                                                              agGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGa             <  1:366119/140‑1 (MQ=255)
                                                                                                                                              agGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGa             >  2:366119/1‑140 (MQ=255)
                                                                                                                                                gTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGc  <  1:223431/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGCCGTTGAAGTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGC  >  NC_000913/3809536‑3809828

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