Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 999,964 T→C 5.2% V250A (GTA→GCA)  elfC → putative outer membrane fimbrial subunit export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913999,9640TC 5.2% 245.5 / 3.3 77V250A (GTA→GCA) elfCputative outer membrane fimbrial subunit export usher protein
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base C (2/2);  total (37/40)
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.63e-01

AATTTAAACAGTGGCGTTAATTTAGGCCCATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCC  >  NC_000913/999819‑1000104
                                                                                                                                                 |                                                                                                                                            
aaTTTAAACAGTGGCGTTAATTTAGGCCCATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTTGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCATCTGAAAGGCGAACTGACCGTa                                                                                                                                             <  3:65228/147‑1 (MQ=255)
    tAAACAGTGGCGTTAATTTAGGCCCATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGtg                                                                                                                                         <  3:360537/147‑1 (MQ=255)
         aGTGGCGTTAATTTAGGCCCATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGGTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGAtt                                                                                                                                   >  7:120172/1‑148 (MQ=255)
            ggCGTTAATTTAGGCCCATGGAGATTCCGCAACAATTAAACATGGAGTCGTAGTAGTGGCCAAAACGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATtata                                                                                                                                >  1:22179/1‑148 (MQ=255)
              cGTTAATTTAGGCCCATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACt                                                                                                                              <  2:387603/148‑1 (MQ=255)
                            cATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGAtttttt                                                                                                                >  6:195826/1‑148 (MQ=255)
                            cATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGAtttttt                                                                                                                <  5:195826/148‑1 (MQ=255)
                             aTGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTc                                                                                                               >  8:127664/1‑148 (MQ=255)
                               ggAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGa                                                                                                             >  1:163726/1‑148 (MQ=255)
                                 agaTTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATa                                                                                                           <  5:331300/148‑1 (MQ=255)
                                 agaTTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATa                                                                                                           <  8:119816/148‑1 (MQ=255)
                                    ttGCGCACCAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGCCCGTAGGTGATGATTATACTGCTGGCGAGTTTTTCGATAgtg                                                                                                        >  4:197103/1‑148 (MQ=255)
                                            caaTTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCttt                                                                                                <  2:22179/148‑1 (MQ=255)
                                               ttCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCgt                                                                                             <  2:163726/148‑1 (MQ=255)
                                                      tGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGc                                                                                      <  7:127664/148‑1 (MQ=255)
                                                      tGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGc                                                                                      <  4:166902/148‑1 (MQ=255)
                                                       ggAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGtgtg                                                                                       >  7:229642/1‑146 (MQ=255)
                                                        gAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGAt                                                                                                            <  8:252969/124‑1 (MQ=255)
                                                        gAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGAt                                                                                                            >  7:252969/1‑124 (MQ=255)
                                                            cGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTgg                                                                                >  3:21261/1‑148 (MQ=255)
                                                            cGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTgg                                                                                >  5:77950/1‑148 (MQ=255)
                                                            cGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTgg                                                                                >  5:379668/1‑148 (MQ=255)
                                                            cGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTgg                                                                                >  4:83740/1‑148 (MQ=255)
                                                               agtagtGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGAtttttt                                                                                                                >  2:79968/1‑113 (MQ=255)
                                                               agtagtGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGAtttttt                                                                                                                <  1:79968/113‑1 (MQ=255)
                                                                  agtGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTc                                                                                               >  2:265006/1‑127 (MQ=255)
                                                                  agtGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTc                                                                                               <  1:265006/127‑1 (MQ=255)
                                                                  agtGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGt                                                                             >  3:214173/1‑145 (MQ=255)
                                                                  agtGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGt                                                                             >  7:117405/1‑145 (MQ=255)
                                                                   gtGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAg                                                                          <  4:214173/147‑1 (MQ=255)
                                                                    tGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCACAt                                                                        >  7:351068/1‑148 (MQ=255)
                                                                     ggCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGAt                                                                        >  8:138756/1‑147 (MQ=255)
                                                                          aaCCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGacaac                                                                  <  4:21261/148‑1 (MQ=255)
                                                                          aaCCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGacaac                                                                  <  7:138756/148‑1 (MQ=255)
                                                                          aaCCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGacaac                                                                  <  3:83740/148‑1 (MQ=255)
                                                                              gCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGCAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACAtgc                                                              >  6:105121/1‑148 (MQ=255)
                                                                              gCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGCAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACAtgc                                                              >  4:67667/1‑148 (MQ=255)
                                                                              gCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGCAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACAtgc                                                              <  5:105121/148‑1 (MQ=255)
                                                                              gCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGCAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACAtgc                                                              <  3:67667/148‑1 (MQ=255)
                                                                                 gAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACAtgctgc                                                           >  8:206956/1‑148 (MQ=255)
                                                                                     ggAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATtat                                                                                                                                 <  3:87336/74‑1 (MQ=255)
                                                                                     ggAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATtat                                                                                                                                 >  4:87336/1‑74 (MQ=255)
                                                                                      gaagacTCTCAGGAGCTATTTGCAACGGGCGGTTATTCCTCTGATAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACTTCATGCTGCCAGAc                                                      <  4:258574/148‑1 (MQ=255)
                                                                                      gaagaaTCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGa                                                       >  2:272641/1‑147 (MQ=255)
                                                                                      gaagaaTCTCAGCAGCTATTTGCAACAGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGa                                                       <  7:206956/147‑1 (MQ=255)
                                                                                       aagaaTCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACa                                                     <  7:40351/148‑1 (MQ=255)
                                                                                        agaaTCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTa                                                                                                                                             >  4:28917/1‑59 (MQ=255)
                                                                                        agaaTCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTa                                                                                                                                             <  3:28917/59‑1 (MQ=255)
                                                                                                agcagcTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGaa                                              <  8:117405/146‑1 (MQ=255)
                                                                                                   agcTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGt                                         <  6:77950/148‑1 (MQ=255)
                                                                                                      tATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTg                                      >  1:240065/1‑148 (MQ=255)
                                                                                                         ttGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCa                                                                                     <  2:390494/98‑1 (MQ=255)
                                                                                                         ttGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCa                                                                                     >  1:390494/1‑98 (MQ=255)
                                                                                                          tGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATg                                                                       >  2:142451/1‑111 (MQ=255)
                                                                                                          tGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATg                                                                       <  1:142451/111‑1 (MQ=255)
                                                                                                           gCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGcc                                  <  8:229642/147‑1 (MQ=255)
                                                                                                                  gCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGc                          <  1:272641/148‑1 (MQ=255)
                                                                                                                      ttATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGa                                                                      >  6:79736/1‑100 (MQ=255)
                                                                                                                      ttATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGa                                                                      <  5:79736/100‑1 (MQ=255)
                                                                                                                      ttATTCCTCTGAAAAGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTgg                      >  7:42345/1‑148 (MQ=255)
                                                                                                                      gtATGCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGGCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGATTGCGCCTGTGGTGCGTgg                      <  3:197103/147‑1 (MQ=255)
                                                                                                                       tATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATg                                                                       >  7:94606/1‑98 (MQ=255)
                                                                                                                       tATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATg                                                                       <  8:94606/98‑1 (MQ=255)
                                                                                                                       tATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTgg                      >  5:39895/1‑147 (MQ=255)
                                                                                                                        aTTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGt                     <  6:379668/147‑1 (MQ=255)
                                                                                                                         ttCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAgt                                                   <  2:8325/116‑2 (MQ=255)
                                                                                                                         ttCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAgt                                                   >  1:8325/1‑115 (MQ=255)
                                                                                                                           ccTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTg                       <  6:383862/142‑1 (MQ=255)
                                                                                                                           ccTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTg                       >  5:383862/1‑142 (MQ=255)
                                                                                                                             tctGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGc                >  3:299799/1‑147 (MQ=255)
                                                                                                                              ctGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGcc               >  7:370392/1‑147 (MQ=255)
                                                                                                                                  aaGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAg          >  6:309921/1‑148 (MQ=255)
                                                                                                                                    ggCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCa        <  7:344641/148‑1 (MQ=255)
                                                                                                                                    ggCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCa        <  5:309921/148‑1 (MQ=255)
                                                                                                                                    ggCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCAAAAAGc         >  7:314882/1‑147 (MQ=255)
                                                                                                                                    ggCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGGGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCa        <  8:114634/148‑1 (MQ=255)
                                                                                                                                          cTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGccc  <  8:120172/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
AATTTAAACAGTGGCGTTAATTTAGGCCCATGGAGATTGCGCAACAATTCAACATGGAGTCGTAGTAGTGGCCAAACCGCAGAATGGAAGAATCTCAGCAGCTATTTGCAACGGGCGGTTATTCCTCTGAAAGGCGAACTGACCGTAGGTGATGATTATACTGCAGGCGATTTTTTCGATAGTGTCAGCTTTCGTGGTGTGCAGCTGGCGTCAGATGACAACATGCTGCCAGACAGCCTGAAAGGGTTTGCGCCTGTGGTGCGTGGTATCGCCAAAAGCAATGCCC  >  NC_000913/999819‑1000104

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: