Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,019,460 G→A V359I (GTC→ATC)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,019,4600GA98.5% 196.2 / NA 65V359I (GTC→ATC) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/34);  minor base T (1/0);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACA  >  NC_000913/2019316‑2019600
                                                                                                                                                |                                                                                                                                            
aGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCAtctc                                                                                                                                          <  1:1446543/149‑1 (MQ=255)
 gTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATCCGCTTCCGGTTCCCATctct                                                                                                                                         >  2:955738/1‑149 (MQ=255)
  tGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCtt                                                                                                                                        <  1:446542/149‑1 (MQ=255)
      aTTCAGTTAGGGCAATGCAACATCAGTGGAGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTACCGGTTCCCATCTCTTTACa                                                                                                                                    >  1:1334940/1‑149 (MQ=255)
           gTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGCCGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGc                                                                                                                               >  2:270813/1‑149 (MQ=255)
            ttAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCCACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCg                                                                                                                              >  1:835910/1‑149 (MQ=255)
              cgctcAAAGCAACATCAGTGGCGAAAGCGTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCggg                                                                                                                            <  1:997825/145‑3 (MQ=255)
              aGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtg                                                                                                                            >  1:901665/1‑149 (MQ=255)
               gggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtgt                                                                                                                           <  1:1102488/149‑1 (MQ=255)
                ggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGGGTa                                                                                                                          >  2:943329/1‑149 (MQ=255)
                  cAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAAc                                                                                                                        >  2:746474/1‑149 (MQ=255)
                   aaaGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCACCAACAGCAAACCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGAGGATACGCTTCCGGTTCCCATCTGTTTACAAGGGAGTTTAACa                                                                                                                       <  2:1164412/149‑1 (MQ=255)
                    aaGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGGTCCCATCTCTTTACAAGGGCGTGTAAccg                                                                                                                      >  2:1361572/1‑147 (MQ=255)
                         aCATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGCCGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAAGAGGCAAc                                                                                                                 >  2:329497/1‑149 (MQ=255)
                         aCATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAAc                                                                                                                 >  2:631962/1‑149 (MQ=255)
                               gTGGCGAAAGCTTTAGTGGGCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGTCGTGTAACAGGCAACAgcggc                                                                                                           <  1:120536/149‑1 (MQ=255)
                                tGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAgcggcg                                                                                                          >  1:1259695/1‑149 (MQ=255)
                                  gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                        <  2:64175/149‑1 (MQ=255)
                                  gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                        <  2:547860/149‑1 (MQ=255)
                                          tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                                                                                >  2:1212127/1‑149 (MQ=255)
                                           ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                               <  1:257181/149‑1 (MQ=255)
                                            tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCg                                                                                              >  2:415630/1‑149 (MQ=255)
                                             aGTGGTCAGCAGCATGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGTCGACGATACGCTTGCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCAGCGTTGATATTTTCGc                                                                                             <  1:659560/149‑1 (MQ=255)
                                                 gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa                                                                                         <  1:746474/149‑1 (MQ=255)
                                                 gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGGGTAACAGGCAACCGCGGCGTTGAGATTTTCGCCGaa                                                                                         >  1:1385502/1‑149 (MQ=255)
                                                  tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                                                                                        <  1:584096/149‑1 (MQ=255)
                                                  tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGAGATTTTCGCCTAAc                                                                                        >  1:115425/1‑149 (MQ=255)
                                                  gcagctgctgGCCGCTTCCCCGCAACAGCAAAGCCAACGCACAGCAAACCATGTACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGGTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                                                                                        <  2:1006334/148‑1 (MQ=255)
                                                    agcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCTTCTCTTTACAAGGGCGTGGAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACgg                                                                                      >  2:610797/1‑148 (MQ=255)
                                                      cagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATA‑TTTCGCCTAACGTCag                                                                                   <  1:1108528/149‑1 (MQ=255)
                                                            gCCGCTTCCCAGCAACAGCATAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGCCGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCCACAGCGTCGTGGATATTTTCGCCTTACGTCAGGGGTa                                                                              >  1:1137846/1‑149 (MQ=255)
                                                            gCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa                                                                              >  1:1066860/1‑149 (MQ=255)
                                                            gCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa                                                                              >  2:786627/1‑149 (MQ=255)
                                                              cGCTTCCCAGCAACAGCAAAGCCAACGCCCCGCAAACCATGACCCTCTGGCGGGGGAAGACGCCGATACGCTTCCGGTTCCCATCTCGTGACAAGGGCGTGTAACAGGCAACAGCTGCGTGGATATTTTCGCCGGACGTCAGAGGTAGc                                                                            >  1:816762/1‑149 (MQ=255)
                                                                         aacagcaaAGCCAACGCTCAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTACCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGATGATATTTTCGCCTAACGAAAGAGGTAGCACCGTAATCcg                                                                 >  1:522944/1‑149 (MQ=255)
                                                                          acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCcgc                                                                >  2:1423259/1‑149 (MQ=255)
                                                                          acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCcgc                                                                <  2:835910/149‑1 (MQ=255)
                                                                             gcaaAGCCAACGCACAGCAACCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTc                                                             <  1:795472/149‑1 (MQ=255)
                                                                             gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTc                                                             >  1:961697/1‑149 (MQ=255)
                                                                              caaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCACCAGCGGCGTTGATATTTTCGCCTACCGTCAGAGGTAGCACCGTAATCCGCGTCt                                                            >  1:116149/1‑149 (MQ=255)
                                                                                         cacaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCtt                                                 <  1:415630/149‑1 (MQ=255)
                                                                                          acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTGCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCttt                                                <  1:594458/149‑1 (MQ=255)
                                                                                           caGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttg                                               <  1:631962/149‑1 (MQ=255)
                                                                                             gCAAACCATGATCCTATGGCGGGGGAAGAAGACGATTCGCTTCCGGTTACCATCTCTTTACAAGGGCGTGTATCAGGAAAAAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCTCCGTAATCCGCGTCTTTTCCCCGCTttgtt                                             <  1:236772/149‑1 (MQ=255)
                                                                                             gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAAAAGCGGCGGTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                             <  1:270040/149‑1 (MQ=255)
                                                                                              cAAACCATGAACCGCTGGCGGGGGAAGACGACGATACGCGTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgttg                                            <  1:1361572/149‑1 (MQ=255)
                                                                                                 aCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTgcgc                                         <  2:115425/149‑1 (MQ=255)
                                                                                                   cATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTc                                       >  1:922921/1‑149 (MQ=255)
                                                                                                    aTGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCa                                      >  1:1178254/1‑149 (MQ=255)
                                                                                                      gAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAg                                    >  1:960979/1‑149 (MQ=255)
                                                                                                           tctGGCGGGGGAAGACGACGATACGCTTCCTGTTCCCATCTCTTTACAAGGGCGTGAAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCa                               <  2:437217/149‑1 (MQ=255)
                                                                                                           tctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTAGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCa                               >  1:1106895/1‑149 (MQ=255)
                                                                                                           tctGGAGGGGGAAGCCGACGATACGCTTCCGGTTCCCATCTCTGTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATTGTGTCGCCTAACGGCAGAGGTAGGACCGTAATCCGCGTCTGTGCCCGGCTTTGTTGCGCTCAAGgcgcg                               >  2:823184/1‑148 (MQ=255)
                                                                                                            ctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAg                              <  2:1181738/149‑1 (MQ=255)
                                                                                                            ctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAg                              <  2:1178254/149‑1 (MQ=255)
                                                                                                                 gggggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTGGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATaa                         <  1:656481/149‑1 (MQ=255)
                                                                                                                   gggAAGACGACGATACTCTTCAGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAtt                       <  1:349330/149‑1 (MQ=255)
                                                                                                                   gggAAGAAGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCCCCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAtt                       <  1:343791/149‑1 (MQ=255)
                                                                                                                     gAAGACGACGATACGCTTCCGGATCCCATCTCTTTCCAAGGGCGTGGACCCGGCAACCGCGGCGTGGATATTTTCGCCTAACGTCAGAGGTAGCACCGTCATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAg                     <  2:876393/149‑1 (MQ=255)
                                                                                                                      aaGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGc                    <  1:1371060/149‑1 (MQ=255)
                                                                                                                      aaGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGc                    >  1:231541/1‑149 (MQ=255)
                                                                                                                       aGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGcc                   <  2:1259695/149‑1 (MQ=255)
                                                                                                                               aTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCa           <  2:901665/149‑1 (MQ=255)
                                                                                                                                    cTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGc      >  1:1017678/1‑149 (MQ=255)
                                                                                                                                        cGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGaca  <  1:943329/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
AGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACA  >  NC_000913/2019316‑2019600

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