Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,799,671 C→T 5.5% A176T (GCA→ACA)  waaZ ← lipopolysaccharide KdoIII transferase; lipopolysaccharide core biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,799,6710CT 5.5% 243.6 / 2.9 73A176T (GCA→ACA) waaZlipopolysaccharide KdoIII transferase; lipopolysaccharide core biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base T (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

AAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTCTTTATAA  >  NC_000913/3799525‑3799817
                                                                                                                                                  |                                                                                                                                                  
aaaTGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTgca                                                                                                                                                   <  6:241557/148‑1 (MQ=255)
   tGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTgcagca                                                                                                                                                <  3:275763/148‑1 (MQ=255)
     gCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCaa                                                                                                                                               <  8:91894/147‑1 (MQ=255)
           ttttGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCaa                                                                                                                                         <  6:266988/147‑1 (MQ=255)
               gAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  6:65910/147‑1 (MQ=255)
               gAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  8:59259/147‑1 (MQ=255)
               gAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  2:47445/147‑1 (MQ=255)
               gAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  3:119999/147‑1 (MQ=255)
                aaTAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTAt                                                                                                                                   <  3:294335/148‑1 (MQ=255)
                           aCTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  8:282032/135‑1 (MQ=255)
                           aCTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     >  7:282032/1‑135 (MQ=255)
                                      aGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGAc                                                                                                              >  2:294536/1‑147 (MQ=255)
                                            tGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACTAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCATCAAAGGAAATAGTATGGCATGAGCAATAACCAAGACTAATGt                                                                                                        >  8:124026/1‑147 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGGACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACAAAGACTAATGTcc                                                                                                      >  8:161033/1‑148 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAAt                                                                                                                                        >  4:134961/1‑114 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAAt                                                                                                                                        <  3:134961/114‑1 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  6:206972/117‑1 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     >  5:206972/1‑117 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGc                                                                                                                                <  5:29500/122‑1 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGc                                                                                                                                >  6:29500/1‑122 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTcc                                                                                                      >  1:210855/1‑148 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTcc                                                                                                      >  5:131248/1‑148 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTcc                                                                                                      >  3:44572/1‑148 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTcc                                                                                                      >  3:48658/1‑148 (MQ=255)
                                             gggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGAGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACAAAGACTAATGTcc                                                                                                      >  7:131215/1‑148 (MQ=255)
                                              ggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCa                                                                                                                         <  5:154732/128‑1 (MQ=255)
                                              ggATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCa                                                                                                                         >  6:154732/1‑128 (MQ=255)
                                               gATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCtt                                                                                                    >  2:115046/1‑148 (MQ=255)
                                                  ttttATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGc                                                                                                  >  6:294398/1‑147 (MQ=255)
                                                    ttATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCaa                                                                                                <  1:115046/147‑1 (MQ=255)
                                                       ttCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATaa                                                                                                                     >  1:64963/1‑123 (MQ=255)
                                                       ttCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATaa                                                                                                                     <  2:64963/123‑1 (MQ=255)
                                                        tctcATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATc                                                                                            >  8:108091/1‑147 (MQ=255)
                                                          tcatcaTAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCaa                                                                                         <  1:294536/148‑1 (MQ=255)
                                                           catcatAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAAc                                                                                        >  8:264739/1‑148 (MQ=255)
                                                            atcatAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCa                                                                                          <  6:131248/145‑1 (MQ=255)
                                                             tcatAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGTAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAAc                                                                                        >  5:253616/1‑146 (MQ=255)
                                                              catAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAg                                                                                                                <  3:262639/121‑1 (MQ=255)
                                                              catAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAg                                                                                                                >  4:262639/1‑121 (MQ=255)
                                                              catAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCAATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAAc                                                                                        >  7:266846/1‑145 (MQ=255)
                                                                        agaACAGCTACCCGTTAAATCAAGACCAGAACAATTAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCtttt                                                                           <  4:48658/148‑1 (MQ=255)
                                                                         gaACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCtttt                                                                           <  4:44572/147‑1 (MQ=255)
                                                                          aaCAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCtttttt                                                                         <  7:264739/148‑1 (MQ=255)
                                                                            cAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAg                                                                       <  2:210855/148‑1 (MQ=255)
                                                                                            cAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     >  7:187078/1‑70 (MQ=255)
                                                                                            cAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  8:187078/70‑1 (MQ=255)
                                                                                             aaGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCAATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATc                                                      <  8:266846/148‑1 (MQ=255)
                                                                                                   aGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     >  1:314543/1‑63 (MQ=255)
                                                                                                   aGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGt                                                                                                                                     <  2:314543/63‑1 (MQ=255)
                                                                                                   aGAACAAATAATCCGCGCAAACTTAAGTGAATATGCGAATTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGAAGACCTTTTTTAGACAACGGTACAGAGATCAGCAGt                                                >  4:123711/1‑148 (MQ=255)
                                                                                                    gAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTATAGAAAACGGTACAGAGATcagcagat                                               >  4:43345/1‑146 (MQ=255)
                                                                                                      aCAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCtttttt                                                                         <  8:200605/120‑1 (MQ=255)
                                                                                                      aCAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCtttttt                                                                         >  7:200605/1‑120 (MQ=255)
                                                                                                          ataatCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTa                                                                                                            >  7:178505/1‑81 (MQ=255)
                                                                                                          ataatCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTa                                                                                                            <  8:178505/81‑1 (MQ=255)
                                                                                                           taatCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATcagca                                                  >  3:292698/1‑138 (MQ=255)
                                                                                                           taatCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATcagca                                                  <  4:292698/138‑1 (MQ=255)
                                                                                                           taatCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTtgt                                         >  1:223936/1‑147 (MQ=255)
                                                                                                             atCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATTGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGa                                       <  5:294398/147‑1 (MQ=255)
                                                                                                                 gcgcATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTg                                  <  3:43345/148‑1 (MQ=255)
                                                                                                                    cATACTTAAGTGAATATGCGATTTGAATTGTAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGtct                               <  6:253616/148‑1 (MQ=255)
                                                                                                                       aattAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATCATCATTAACCAAGACAAATGTCCTAGCAACATCAAACCAGACGACCATATTAAGACAAAGGTACAGAGATCAGCAGTTCTTTTTGAATCAGTCTCa                             >  2:195692/3‑147 (MQ=255)
                                                                                                                       aCTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCa                             >  8:184901/1‑147 (MQ=255)
                                                                                                                            aGTGAATATGCGATTTGAATTGTAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAAt                                     >  3:132421/1‑134 (MQ=255)
                                                                                                                            aGTGAATATGCGATTTGAATTGTAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAAt                                     <  4:132421/134‑1 (MQ=255)
                                                                                                                              tGAATATGCGATTTGAATTTCAGAAAAGGCAATAGTATGTAATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGAAGACCTTTTTTAGACAACGGTACAGAGATCAGAAGTTCTTTGTGAATCTGTCTCAAAATAtt                      >  2:5404/1‑147 (MQ=255)
                                                                                                                                  tatGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAAtt                 >  2:147382/1‑148 (MQ=255)
                                                                                                                                   atGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAAtt                 >  4:287687/1‑147 (MQ=255)
                                                                                                                                   atGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAAtt                 >  5:21155/1‑147 (MQ=255)
                                                                                                                                    tGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAAttt                >  8:126223/1‑147 (MQ=255)
                                                                                                                                    tGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAAttt                <  7:126223/147‑1 (MQ=255)
                                                                                                                                      cGAATTGATTTTCATCAAAGGCAATTGTATGGCATGAGCAATATCAAAGACTAATTTCCTTGCAAAATCCAACTAGACTACCTTTTTTAGACAACGGTACAGAGCTCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTaa              <  8:131215/147‑1 (MQ=255)
                                                                                                                                         tttGAATTGCAGCAACGGCAATAGTAGGGCATGAGCAATAACCAAGACTACTGTCATTTCAAAATACAACAAGACGAACTTTATTAGAAAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTAAAAATATTAAATTTAAtctt          <  1:5404/148‑1 (MQ=255)
                                                                                                                                         tttGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAAtctt          >  3:57735/1‑148 (MQ=255)
                                                                                                                                             aaTTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATcttctt       >  5:240276/1‑147 (MQ=255)
                                                                                                                                                  cagcaAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTCTTTATaa  <  7:184901/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTCTTTATAA  >  NC_000913/3799525‑3799817

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