Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC50.9% 68.2 / 94.5 58intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (25/3);  new base C (0/29);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 3.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-08

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001769‑2002064
                                                                                                                                                     |                                                                                                                                                   
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑c                                                                                                                                                   >  2:672424/1‑151 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑cc                                                                                                                                                  >  2:321112/1‑151 (MQ=255)
        agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCgg                                                                                                                                           >  1:525372/1‑151 (MQ=255)
         gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGt                                                                                                                                          >  1:134204/1‑151 (MQ=255)
              atttCACGAAAAACAGCCCTGCGTTATATGATTTATGGGAATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGAGGGTTTTTTTGGGGGTTACAGTTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCg                                                                                                                                      <  1:378043/148‑1 (MQ=255)
                   aCGAAAACCACCCCTGCTTTATAGGATTTCTGGGCTTGATTTCCCGTTTCTGCAGGGTTTTAAATCGGACGATTGGTGGGAGAAATGAGGGTTTTTTTGGGGGTAACGGGAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAg                                                                                                                                <  1:447076/151‑1 (MQ=255)
                        aaaCAGCCCTGTGTTTTATGATTTATCGGCTTGATTTCCGGTTTCTGCGGTTTTTTTAATCGGACATTTGGTGGGTAAAAGGGGGGTTTTTTTGGGGTTAACAGGAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                           <  1:201036/151‑1 (MQ=255)
                          aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                          >  1:196649/1‑150 (MQ=255)
                             ccccTGCTTTATATAGTTTATCGGCATTTTTATCCTTTTCTGCGGGTTTTTTAATCGAACAATTAGTGGTTAAAAGGAGGGTTTTTTTGGGGGTTACAGTTAAATCCCAGCAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                      <  1:17102/150‑1 (MQ=255)
                                 tGCTTTATAGGAGTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTTTTTGGGGGTTACAGGAAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                   <  1:521701/150‑1 (MQ=255)
                                  gCGTTATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGTTTTTAAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                   <  2:698012/149‑1 (MQ=255)
                                    gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                               >  1:86629/1‑151 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTGGTGGGTAAAAGGGGGGGTTTTTTGGGGGTTACAGGTAATTCCAAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                               <  1:899282/150‑1 (MQ=255)
                                      tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                             >  2:1169298/1‑151 (MQ=255)
                                          tGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGTTTATTTGGGGTTTACAGGAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                          <  1:159310/150‑1 (MQ=255)
                                             gTTATCGGCATGATTATCTGTTTGTGCGGGGTTTTTAATCGGAGGATTGGTGGGTGAAAGGGGGGGTTTTTTGGGGGTGACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                       <  1:290180/150‑1 (MQ=255)
                                              ttACCGGCATAATTATCCTTTTCTGCGGGTTTTTTAATCGGAAAATTAGGGGGTAAAATGGGGGTTTATTGGGGGTTTAAAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                         <  2:677215/147‑1 (MQ=255)
                                                 tCGGAATGATTATCCGTTTCGGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGTTTACAGGAAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                   <  2:153584/150‑1 (MQ=255)
                                                    gaaTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                               <  2:104618/149‑1 (MQ=255)
                                                      aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                              >  1:707596/1‑150 (MQ=255)
                                                             tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      >  2:883522/1‑151 (MQ=255)
                                                             tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      >  2:606196/1‑151 (MQ=255)
                                                              ctttttttgggggttttttAGTCGGACGTTTAGGGGGTAAAAGGGGGGTTTTTTTGGGGGTTAAAGGT‑AATTTCAGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                        <  1:683629/134‑1 (MQ=255)
                                                                gTTTCTGGAGGTTTTTAAACCGGAAGATTAGGGGGAAAAAGGAGGGTTTTTTTGGGGTTAAAAGGTAAATCCCGGGAAAAAAAAAG‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                       <  1:786763/147‑1 (MQ=255)
                                                                 tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGACGCGGCGTaa                                                                                   >  2:1066919/1‑150 (MQ=255)
                                                                    cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               >  1:441388/1‑151 (MQ=255)
                                                                      gCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTTTTTGGGGGTAACAGTAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               <  1:746249/149‑1 (MQ=255)
                                                                      gCAGGGTTTTTAATGGGACGATTAGGGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               <  1:397717/149‑1 (MQ=255)
                                                                      gCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                             >  2:202920/1‑151 (MQ=255)
                                                                        aGGGTTTTAAATCGGACGATTGGGGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                            <  2:302406/150‑1 (MQ=255)
                                                                         gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          >  2:942575/1‑151 (MQ=255)
                                                                              tttAATCGGACGTTTGGTGGGGGAAATGAGGGGTTTTTGGGGGTTTACAGTAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                      <  2:489947/150‑1 (MQ=38)
                                                                                 aaTCGAACGATTAGGGGGTGAAAGGAGGTTTTTTTTGGGGTTAACAGGTAAATTCAAGGCAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                  <  1:426535/151‑1 (MQ=37)
                                                                                    cgaagtattggggggTAAAATGAGGGTTTTTTTGGGGGTGAAAGGAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                               <  1:618756/139‑1 (MQ=25)
                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                >  1:1028755/1‑150 (MQ=255)
                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                          >  2:1059523/1‑151 (MQ=255)
                                                                                           ttggggggTGAAAGGAGGGTTTATTTGGGGGTTACAGATAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                         <  1:7757/145‑1 (MQ=38)
                                                                                            tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                       >  1:545097/1‑151 (MQ=255)
                                                                                              gTGGGTGAAAGGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                       <  1:969478/149‑1 (MQ=38)
                                                                                                gggaaaaaggggggtTTTTTTGGGGGTAACAGTAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                        <  1:225543/135‑1 (MQ=25)
                                                                                                  gTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                 >  2:578998/1‑151 (MQ=255)
                                                                                                           gggtTTTTTTGGGGTTTAAAGGAAAATTCAAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                           <  1:604040/148‑1 (MQ=25)
                                                                                                           gggtTTATTTGGGGGTTACAGGAAAATGCCGGGCAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  2:522088/151‑1 (MQ=255)
                                                                                                           ggggTTTTTTGGGGGTTACAGAAAAATCCGAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  1:41364/151‑1 (MQ=25)
                                                                                                           ggggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        >  2:700545/1‑151 (MQ=255)
                                                                                                                tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                   >  1:602637/1‑151 (MQ=255)
                                                                                                                  tttGGGGGTTACAGGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                   <  1:68081/149‑1 (MQ=255)
                                                                                                                  ttggggggTTACAGGTAAATCCAAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                   <  2:349829/146‑1 (MQ=255)
                                                                                                                   ttGGGGTTAACAGTTAAATTCCAGCAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                 <  1:242561/150‑1 (MQ=37)
                                                                                                                       gggTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                            >  1:207876/1‑151 (MQ=255)
                                                                                                                        ggTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           >  2:662534/1‑151 (MQ=255)
                                                                                                                        ggTTAAAGGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           <  2:115368/151‑1 (MQ=255)
                                                                                                                            aCAGTTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                        <  1:8339/150‑1 (MQ=255)
                                                                                                                              aGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGAGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                     >  2:1172875/1‑151 (MQ=255)
                                                                                                                                  aaaTTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                  <  2:1018763/150‑1 (MQ=255)
                                                                                                                                   aaTTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa                <  1:321097/151‑1 (MQ=255)
                                                                                                                                           gCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg        >  2:981040/1‑151 (MQ=255)
                                                                                                                                             aaaaaaacc‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAgcgcagcgca      <  1:95071/149‑1 (MQ=255)
                                                                                                                                             aGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc       >  1:314323/1‑150 (MQ=255)
                                                                                                                                                aaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa   <  1:606176/151‑1 (MQ=255)
                                                                                                                                                aaaaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc    >  1:197964/1‑150 (MQ=255)
                                                                                                                                                 aaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg  <  1:116347/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001769‑2002064

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