Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,964,425 G→C 7.5% A209G (GCC→GGC)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,4250GC 7.5% 364.4 / 17.4 107A209G (GCC→GGC) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (48/51);  new base C (4/4);  total (52/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

Only 100 of 109 total aligned reads displayed.

GCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTG  >  NC_000913/1964281‑1964561
                                                                                                                                                |                                                                                                                                        
gcCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATc                                                                                                                                      >  1:194410/1‑149 (MQ=255)
     aaTTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccgg                                                                                                                                 <  2:194410/149‑2 (MQ=255)
            gCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACaaa                                                                                                                          >  2:311905/1‑149 (MQ=255)
                ccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACaaa                                                                                                                          <  4:1102457/145‑1 (MQ=255)
                   gACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAacttactt                                                                                                                   <  1:730914/149‑1 (MQ=255)
                        tCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGcccc                                                                                                              <  3:225717/149‑1 (MQ=255)
                           ccAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTc                                                                                                           <  3:660822/149‑1 (MQ=255)
                             aaTGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCa                                                                                                         >  3:851571/1‑149 (MQ=255)
                              aTGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCAt                                                                                                        <  4:748927/149‑1 (MQ=255)
                               tGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCAtt                                                                                                       <  3:150819/149‑1 (MQ=255)
                               cgatctaTAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGTCGCGAACAAACTTACTTGCCCCGTCCAtt                                                                                                       <  3:875980/143‑1 (MQ=255)
                                    aaTAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATgg                                                                                                                                          <  4:1464541/109‑1 (MQ=255)
                                    aaTAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATgg                                                                                                                                          >  3:1464541/1‑109 (MQ=255)
                                    aaTAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGcc                                                                                                  >  1:1040818/1‑149 (MQ=255)
                                      tAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACAATGATGAGAATCCCGGCTATGGCATCGCCGCGAACAAACTTACTTGCGCCGTCCATTGAgatcg                                                                                                >  4:875980/1‑145 (MQ=255)
                                           tataACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTcc                                                                                                          >  1:1390951/1‑134 (MQ=255)
                                           tataACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTcc                                                                                                          <  2:1390951/134‑1 (MQ=255)
                                             taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                                                         <  4:851571/149‑1 (MQ=255)
                                             taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                                                         <  2:1040818/149‑1 (MQ=255)
                                              aaCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                                                        <  1:62617/149‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                                                         >  1:265251/1‑113 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                                                         <  2:265251/113‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTc                                                                                                           >  1:1417714/1‑127 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTc                                                                                                           <  1:459879/127‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTc                                                                                                           <  2:1417714/127‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGa                                                                                                     >  3:1401701/1‑133 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGa                                                                                                     <  4:1401701/133‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                                                        >  2:173858/1‑146 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                                                        >  4:73218/1‑146 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                                                        <  3:73218/146‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                                                        >  2:652280/1‑146 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                                                        <  1:652280/146‑1 (MQ=255)
                                                 tttCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                                                        <  1:173858/146‑1 (MQ=255)
                                                    ccgccgCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTcc                                                                                  <  4:835286/149‑1 (MQ=255)
                                                           gtgtCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGa                                                                                                     <  1:770274/123‑1 (MQ=255)
                                                           gtgtCCAATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGACCAAACTTACTTGCCCCGTCCATTGa                                                                                                     >  2:770274/1‑123 (MQ=255)
                                                               cccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTc                                                                       <  3:513375/149‑1 (MQ=255)
                                                                ccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccgg                                                                                                                                 >  4:1218287/1‑89 (MQ=255)
                                                                ccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccgg                                                                                                                                 <  3:1218287/90‑2 (MQ=255)
                                                                 cATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCgg                                                                     >  3:1187903/1‑149 (MQ=255)
                                                                  aTGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGa                                                                    <  3:955539/149‑1 (MQ=255)
                                                                  aTGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGa                                                                    >  2:1348876/1‑149 (MQ=255)
                                                                  aTGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGa                                                                    <  1:1348876/149‑1 (MQ=255)
                                                                     cTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                                                         <  3:453512/125‑1 (MQ=255)
                                                                     cTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                                                         >  4:453512/1‑125 (MQ=255)
                                                                      tCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGt                                                                             <  4:554584/136‑1 (MQ=255)
                                                                      tCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGt                                                                             >  3:554584/1‑136 (MQ=255)
                                                                         tGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGt                                                             >  4:1101389/1‑149 (MQ=255)
                                                                            cATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATcgccgc                                                                                                                                 >  4:506021/1‑78 (MQ=255)
                                                                            cATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATcgccgc                                                                                                                                 <  3:506021/78‑1 (MQ=255)
                                                                            cATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCa                                                                                                         >  3:1092344/1‑102 (MQ=255)
                                                                            cATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCa                                                                                                         <  4:1092344/102‑1 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTc                                                       >  1:1279359/1‑149 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTc                                                       >  4:576547/1‑149 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTc                                                       <  3:1204122/149‑1 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTc                                                       <  3:1101389/149‑1 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATCGCCGGGa                                                                                                                               <  3:1006217/77‑1 (MQ=255)
                                                                               gTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATCGCCGGGa                                                                                                                               >  4:1006217/1‑77 (MQ=255)
                                                                                  gcagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                                                         >  4:206728/1‑80 (MQ=255)
                                                                                  gcagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                                                         <  3:206728/80‑1 (MQ=255)
                                                                                  gcagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCAtt                                                                                                       >  3:1146150/1‑98 (MQ=255)
                                                                                  gcagcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCAtt                                                                                                       <  4:1146150/98‑1 (MQ=255)
                                                                                    agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTaaaaa                                                                                           <  1:1234556/108‑1 (MQ=255)
                                                                                    agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTaaaaa                                                                                           >  2:1234556/1‑108 (MQ=255)
                                                                                    agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc                                                    <  2:192042/147‑1 (MQ=255)
                                                                                    agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc                                                    >  1:192042/1‑147 (MQ=255)
                                                                                    agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTc                                                  >  1:126144/1‑149 (MQ=255)
                                                                                        cGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCAt                                                                                                        >  2:800545/1‑91 (MQ=255)
                                                                                        cGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCAt                                                                                                        <  1:800545/91‑1 (MQ=255)
                                                                                            gACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAc                                          <  3:576547/149‑1 (MQ=255)
                                                                                            gACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAc                                          <  2:126144/149‑1 (MQ=255)
                                                                                                 gCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAc                                          >  2:684831/1‑144 (MQ=255)
                                                                                                 gCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAc                                          <  1:684831/144‑1 (MQ=255)
                                                                                                     cgcgccgACAATGTTAATGACCATGATGATGAGCCCGGCGATGGCATCGCCGCGAACAAACGTACTTGCCCCATCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTCTTCGCCTCATCTTCACcaatcaagc                                 <  4:463857/147‑3 (MQ=255)
                                                                                                      ccgccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaatcaat                                  >  1:1180360/1‑147 (MQ=255)
                                                                                                        gccgACAATGTTAATGACCATGATGAGGGTCCCGGCGATTGCCTCGCCACGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCTGCTTCCTGAGTCACTTCGGAGCGCCGTTTTTTCGCCTCATCTTCACCAATCACTccag                              >  3:463857/1‑147 (MQ=255)
                                                                                                           gACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtt                                                            >  1:604042/1‑116 (MQ=255)
                                                                                                           gACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtt                                                            <  2:604042/116‑1 (MQ=255)
                                                                                                            aCAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTg                                                                                                      <  1:797141/73‑1 (MQ=255)
                                                                                                            aCAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTg                                                                                                      >  2:797141/1‑73 (MQ=255)
                                                                                                              aaTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTc                                                       >  3:89563/1‑118 (MQ=255)
                                                                                                              aaTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTaa                        <  2:1279359/149‑1 (MQ=255)
                                                                                                              aaTGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGACTCACTTCGGAGCGGCGTTTTTTc                                                       <  4:89563/118‑1 (MQ=255)
                                                                                                                tGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAg                                                                                                    <  4:1458171/71‑1 (MQ=255)
                                                                                                                tGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAg                                                                                                    >  3:1458171/1‑71 (MQ=255)
                                                                                                                  ttAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAc                                          >  3:749655/1‑127 (MQ=255)
                                                                                                                  ttAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAc                                          <  4:749655/127‑1 (MQ=255)
                                                                                                                    aaTGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCaatcaatc                                 <  4:441164/134‑1 (MQ=255)
                                                                                                                    aaTGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCaatcaatc                                 >  3:441164/1‑134 (MQ=255)
                                                                                                                     aTGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTc                                            >  1:477847/1‑122 (MQ=255)
                                                                                                                     aTGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTc                                            <  2:477847/122‑1 (MQ=255)
                                                                                                                     aTGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGc                 >  3:1075420/1‑149 (MQ=255)
                                                                                                                            tgatgaGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATc          >  4:391386/1‑149 (MQ=255)
                                                                                                                             gatgaGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGa                                                                    <  4:703709/90‑1 (MQ=255)
                                                                                                                             gatgaGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGa                                                                    >  3:703709/1‑90 (MQ=255)
                                                                                                                                   gATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtttttt                                                        >  3:215565/1‑96 (MQ=255)
                                                                                                                                   gATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtttttt                                                        <  4:215565/96‑1 (MQ=255)
                                                                                                                                    aTCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTg  >  3:91955/1‑149 (MQ=255)
                                                                                                                                      cgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaatcaat                                  >  3:534223/3‑115 (MQ=255)
                                                                                                                                      cgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACcaatcaat                                  <  4:534223/113‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
GCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTG  >  NC_000913/1964281‑1964561

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