Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,022,388 T→A 5.9% Y174N (TAC→AAC)  fliP → flagellar biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,022,3880TA 5.9% 220.9 / 3.7 68Y174N (TAC→AAC) fliPflagellar biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base A (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GGAGGCGCTGGAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACT  >  NC_000913/2022240‑2022536
                                                                                                                                                    |                                                                                                                                                    
ggaggCGCTGGAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCt                                                                                                                                                      <  5:217873/149‑1 (MQ=255)
   ggCGCTGGAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACg                                                                                                                                                   >  2:227600/1‑149 (MQ=255)
     cgcTGGAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTg                                                                                                                                                 <  3:136687/149‑1 (MQ=255)
      gcTGGAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGa                                                                                                                                                >  4:91118/1‑149 (MQ=255)
       cTGGAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGAc                                                                                                                                               >  5:114410/1‑149 (MQ=255)
          gAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAg                                                                                                                                            <  1:335897/149‑1 (MQ=255)
           aaaaaGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGc                                                                                                                                           <  8:262120/149‑1 (MQ=255)
                ggggCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGtt                                                                                                                                      >  6:168719/1‑149 (MQ=255)
                       aGCCGCTGCGTGACTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAAcc                                                                                                                               <  6:48922/149‑1 (MQ=255)
                               cGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTTACCAGCGAGTTGAAAACCGCATTTCa                                                                                                                       <  7:171633/149‑1 (MQ=255)
                               cGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTc                                                                                                                        <  3:91118/148‑1 (MQ=255)
                                gTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCAttt                                                                                                                         >  1:87136/1‑146 (MQ=255)
                                 tGAGTTTATGCTGCGTCAGACCCGTGATGCAGATTTAGGGTTGTTTGCCAGAATGGCGAATACCGGCACGTTGCATGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGGGACCAGCGAGTTGAAAACCGCATTTCAGa                                                                                                                     <  2:70671/149‑1 (MQ=255)
                                   aGTTTATGCTGCGTCGGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATa                                                                                                                   <  6:106613/149‑1 (MQ=255)
                                   aGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGa                                                                                                                     >  4:147907/1‑147 (MQ=255)
                                         tgctgcGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCAACGTGACCAGCGAGTTGAAAGCCGCAt                                                                                                                           <  7:321700/135‑1 (MQ=255)
                                         tgctgcGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCAACGTGACCAGCGAGTTGAAAGCCGCAt                                                                                                                           >  8:321700/1‑135 (MQ=255)
                                         tgctgcGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCAACGTGACCAGCGAGTTGAAAACCGCAt                                                                                                                           >  6:104324/1‑135 (MQ=255)
                                         tgctgcGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCAACGTGACCAGCGAGTTGAAAACCGCAt                                                                                                                           <  5:104324/135‑1 (MQ=255)
                                           ctgcGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCa                                                                                                            <  6:114410/148‑1 (MQ=255)
                                           ctgcGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCa                                                                                                            <  8:34193/148‑1 (MQ=255)
                                               gTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGCCCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGAtt                                                                                                       >  1:304965/1‑149 (MQ=255)
                                                 cAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGAtttt                                                                                                     >  4:114262/1‑149 (MQ=255)
                                                    aCCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCAt                                                                                                  <  7:70373/149‑1 (MQ=255)
                                                     cccGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATc                                                                                                 <  8:248175/149‑1 (MQ=255)
                                                        gTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCt                                                                                              >  5:216101/1‑149 (MQ=255)
                                                        gTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCt                                                                                              >  4:110133/1‑149 (MQ=255)
                                                                gATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGAt                                                                                      <  2:87136/149‑1 (MQ=255)
                                                                  tttAGGGTTGTTTGCCACACTGGAGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGAAGCGCATTTTGCTCCCGGCCTACGTGACGATCGTGTTGAAAACGGCATTTCAGATAGGCTTCACGAGTTTCATCCCTTTTTTGATTa                                                                                    <  6:123494/149‑1 (MQ=255)
                                                                     aGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATc                                                                                  >  8:386967/1‑148 (MQ=255)
                                                                     aGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATc                                                                                  <  1:227600/148‑1 (MQ=255)
                                                                     aGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATc                                                                                  <  1:251025/148‑1 (MQ=255)
                                                                               gCCAGACTGGCGAATACCGGCCCGTTGCAGGGCCCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGa                                                                        >  3:279049/1‑148 (MQ=255)
                                                                                  aGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGc                                                                    >  4:131293/1‑149 (MQ=255)
                                                                                      tGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGACGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACATGGTGATAGCCAGc                                                                >  7:217900/1‑149 (MQ=255)
                                                                                        gCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCgt                                                              >  4:80202/1‑149 (MQ=255)
                                                                                          gAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCgtgt                                                            >  4:222336/1‑149 (MQ=255)
                                                                                           aaTACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGtt                                                           >  3:64802/1‑149 (MQ=255)
                                                                                             tACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGa                                                         >  4:33390/1‑149 (MQ=255)
                                                                                               ccGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGGGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGAt                                                        <  1:52736/148‑1 (MQ=255)
                                                                                                cGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCt                                                                             <  5:215369/126‑1 (MQ=255)
                                                                                                cGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCt                                                                             >  6:215369/1‑126 (MQ=255)
                                                                                                cGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATgg                                                      <  4:64833/149‑1 (MQ=255)
                                                                                                cGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATgg                                                      >  7:18153/1‑149 (MQ=255)
                                                                                                cGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATgg                                                      >  5:320090/1‑149 (MQ=255)
                                                                                                cGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATgg                                                      >  3:335810/1‑149 (MQ=255)
                                                                                                         gCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGGGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTgggg                                              <  2:58770/148‑1 (MQ=255)
                                                                                                                ccTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGgatgatgat                                      <  8:22576/149‑1 (MQ=255)
                                                                                                                ccTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGgatgatgat                                      <  3:114262/149‑1 (MQ=255)
                                                                                                                  tGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCgt                                                              <  6:249335/123‑1 (MQ=255)
                                                                                                                  tGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCgt                                                              >  5:249335/1‑123 (MQ=255)
                                                                                                                      gCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTcc                                <  3:80202/149‑1 (MQ=255)
                                                                                                                              gATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGcca                        >  1:19873/1‑149 (MQ=255)
                                                                                                                              gATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGcca                        >  1:26395/1‑149 (MQ=255)
                                                                                                                              gATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGcca                        >  5:140469/1‑149 (MQ=255)
                                                                                                                              gATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTATTGATTATCGACCTGGt                                                                          >  3:101046/1‑99 (MQ=255)
                                                                                                                              gATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTATTGATTATCGACCTGGt                                                                          <  4:101046/99‑1 (MQ=255)
                                                                                                                                   gcATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGc                          >  3:175197/1‑142 (MQ=255)
                                                                                                                                   gcATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGc                          <  4:175197/142‑1 (MQ=255)
                                                                                                                                          gCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGCCCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGc            >  5:86585/1‑149 (MQ=255)
                                                                                                                                           cTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGGTCCCCCAGCCACCATTGCTCTGcc           <  7:386967/149‑1 (MQ=255)
                                                                                                                                            tCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGccc          >  8:2575/1‑149 (MQ=255)
                                                                                                                                            tCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCAACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGccc          >  1:1035/1‑149 (MQ=255)
                                                                                                                                            tCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCGTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGccc          >  5:333862/1‑149 (MQ=255)
                                                                                                                                              ccGGCCTACGTGACCAGCGAGTTGAAACCCGCATTTCAGATAGGCATCACGAATTTCATCCCTTATTTGATTATCGACATCGTGAAATCCAGCGTGTTGATGGCATTGGTGATGATCATGGTTCCCCCAGCAACAAATGCTCTGCCCtt        >  8:176130/1‑149 (MQ=255)
                                                                                                                                              ccGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCtt        <  4:335810/149‑1 (MQ=255)
                                                                                                                                               cGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCtt        >  2:36938/1‑148 (MQ=255)
                                                                                                                                                  ccTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTaaa    <  5:168719/149‑1 (MQ=255)
                                                                                                                                                  ccTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTaaa    <  6:320090/149‑1 (MQ=255)
                                                                                                                                                    tACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGCCGTGGTGATAGCCAGAGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCCCAATTGCTCTGCCCTTGAAACt  >  2:162321/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAGGCGCTGGAAAAAGGGGCGCAGCCGCTGCGTGAGTTTATGCTGCGTCAGACCCGTGAGGCAGATTTAGGGTTGTTTGCCAGACTGGCGAATACCGGCCCGTTGCAGGGACCTGAAGCCGTGCCGATGCGCATTTTGCTCCCGGCCTACGTGACCAGCGAGTTGAAAACCGCATTTCAGATAGGCTTCACGATTTTCATCCCTTTTTTGATTATCGACCTGGTGATAGCCAGCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACT  >  NC_000913/2022240‑2022536

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: ‑