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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001769‑2002064
                                                                                                                                                     |                                                                                                                                                   
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                   >  2:672424/1‑151 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                  >  2:321112/1‑151 (MQ=255)
        agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                           >  1:525372/1‑151 (MQ=255)
         gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGt                                                                                                                                          >  1:134204/1‑151 (MQ=255)
              atttCACGAAAAACAGCCCTGCGTTATATGATTTATGGGAATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGAGGGTTTTTTTGGGGGTTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCg                                                                                                                                      <  1:378043/148‑1 (MQ=255)
                   aCGAAAACCACCCCTGCTTTATAGGATTTCTGGGCTTGATTTCCCGTTTCTGCAGGGTTTTAAATCGGACGATTGGTGGGAGAAATGAGGGTTTTTTTGGGGGTAACGGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAg                                                                                                                                <  1:447076/151‑1 (MQ=255)
                        aaaCAGCCCTGTGTTTTATGATTTATCGGCTTGATTTCCGGTTTCTGCGGTTTTTTTAATCGGACATTTGGTGGGTAAAAGGGGGGTTTTTTTGGGGTTAACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                           <  1:201036/151‑1 (MQ=255)
                          aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                          >  1:196649/1‑150 (MQ=255)
                             ccccTGCTTTATATAGTTTATCGGCATTTTTATCCTTTTCTGCGGGTTTTTTAATCGAACAATTAGTGGTTAAAAGGAGGGTTTTTTTGGGGGTTACAGTTAAATCCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                      <  1:17102/150‑1 (MQ=255)
                                 tGCTTTATAGGAGTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTTTTTGGGGGTTACAGGAAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                   <  1:521701/150‑1 (MQ=255)
                                  gCGTTATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGTTTTTAAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                   <  2:698012/149‑1 (MQ=255)
                                    gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                               >  1:86629/1‑151 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTGGTGGGTAAAAGGGGGGGTTTTTTGGGGGTTACAGGTAATTCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                               <  1:899282/150‑1 (MQ=255)
                                      tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                             >  2:1169298/1‑151 (MQ=255)
                                          tGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGTTTATTTGGGGTTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                          <  1:159310/150‑1 (MQ=255)
                                             gTTATCGGCATGATTATCTGTTTGTGCGGGGTTTTTAATCGGAGGATTGGTGGGTGAAAGGGGGGGTTTTTTGGGGGTGACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                       <  1:290180/150‑1 (MQ=255)
                                              ttACCGGCATAATTATCCTTTTCTGCGGGTTTTTTAATCGGAAAATTAGGGGGTAAAATGGGGGTTTATTGGGGGTTTAAAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                         <  2:677215/147‑1 (MQ=255)
                                                 tCGGAATGATTATCCGTTTCGGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGTTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                   <  2:153584/150‑1 (MQ=255)
                                                    gaaTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                               <  2:104618/149‑1 (MQ=255)
                                                      aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                              >  1:707596/1‑150 (MQ=255)
                                                             tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      >  2:883522/1‑151 (MQ=255)
                                                             tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      >  2:606196/1‑151 (MQ=255)
                                                              ctttttttgggggttttttAGTCGGACGTTTAGGGGGTAAAAGGGGGGTTTTTTTGGGGGTTAAAGGTAATTTCAGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                        <  1:683629/134‑1 (MQ=255)
                                                                gTTTCTGGAGGTTTTTAAACCGGAAGATTAGGGGGAAAAAGGAGGGTTTTTTTGGGGTTAAAAGGTAAATCCCGGGAAAAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                       <  1:786763/147‑1 (MQ=255)
                                                                 tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGACGCGGCGTaa                                                                                   >  2:1066919/1‑150 (MQ=255)
                                                                    cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               >  1:441388/1‑151 (MQ=255)
                                                                      gCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTTTTTGGGGGTAACAGTAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               <  1:746249/149‑1 (MQ=255)
                                                                      gCAGGGTTTTTAATGGGACGATTAGGGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               <  1:397717/149‑1 (MQ=255)
                                                                      gCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                             >  2:202920/1‑151 (MQ=255)
                                                                        aGGGTTTTAAATCGGACGATTGGGGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                            <  2:302406/150‑1 (MQ=255)
                                                                         gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          >  2:942575/1‑151 (MQ=255)
                                                                              tttAATCGGACGTTTGGTGGGGGAAATGAGGGGTTTTTGGGGGTTTACAGTAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                      <  2:489947/150‑1 (MQ=38)
                                                                                 aaTCGAACGATTAGGGGGTGAAAGGAGGTTTTTTTTGGGGTTAACAGGTAAATTCAAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                  <  1:426535/151‑1 (MQ=37)
                                                                                    cgaagtattggggggTAAAATGAGGGTTTTTTTGGGGGTGAAAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                               <  1:618756/139‑1 (MQ=25)
                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                >  1:1028755/1‑150 (MQ=255)
                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                          >  2:1059523/1‑151 (MQ=255)
                                                                                           ttggggggTGAAAGGAGGGTTTATTTGGGGGTTACAGATAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                         <  1:7757/145‑1 (MQ=38)
                                                                                            tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                       >  1:545097/1‑151 (MQ=255)
                                                                                              gTGGGTGAAAGGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                       <  1:969478/149‑1 (MQ=38)
                                                                                                gggaaaaaggggggtTTTTTTGGGGGTAACAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                        <  1:225543/135‑1 (MQ=25)
                                                                                                  gTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                 >  2:578998/1‑151 (MQ=255)
                                                                                                           gggtTTTTTTGGGGTTTAAAGGAAAATTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                           <  1:604040/148‑1 (MQ=25)
                                                                                                           gggtTTATTTGGGGGTTACAGGAAAATGCCGGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  2:522088/151‑1 (MQ=255)
                                                                                                           ggggTTTTTTGGGGGTTACAGAAAAATCCGAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  1:41364/151‑1 (MQ=25)
                                                                                                           ggggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        >  2:700545/1‑151 (MQ=255)
                                                                                                                tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                   >  1:602637/1‑151 (MQ=255)
                                                                                                                  tttGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                   <  1:68081/149‑1 (MQ=255)
                                                                                                                  ttggggggTTACAGGTAAATCCAAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                   <  2:349829/146‑1 (MQ=255)
                                                                                                                   ttGGGGTTAACAGTTAAATTCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                 <  1:242561/150‑1 (MQ=37)
                                                                                                                       gggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                            >  1:207876/1‑151 (MQ=255)
                                                                                                                        ggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           >  2:662534/1‑151 (MQ=255)
                                                                                                                        ggTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           <  2:115368/151‑1 (MQ=255)
                                                                                                                            aCAGTTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                        <  1:8339/150‑1 (MQ=255)
                                                                                                                              aGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGAGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                     >  2:1172875/1‑151 (MQ=255)
                                                                                                                                  aaaTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                  <  2:1018763/150‑1 (MQ=255)
                                                                                                                                   aaTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa                <  1:321097/151‑1 (MQ=255)
                                                                                                                                           gCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg        >  2:981040/1‑151 (MQ=255)
                                                                                                                                             aaaaaaaccCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAgcgcagcgca      <  1:95071/149‑1 (MQ=255)
                                                                                                                                             aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc       >  1:314323/1‑150 (MQ=255)
                                                                                                                                                aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa   <  1:606176/151‑1 (MQ=255)
                                                                                                                                                aaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc    >  1:197964/1‑150 (MQ=255)
                                                                                                                                                 aaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg  <  1:116347/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  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: