Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913559,1810GA40.3% 90.0 / 25.6 68E70K (GAA→AAA) sfmCputative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (23/17);  new base A (0/27);  total (23/45)
Fisher's exact test for biased strand distribution p-value = 1.84e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAA  >  NC_000913/559034‑559329
                                                                                                                                                   |                                                                                                                                                    
gcCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTaaaagaa                                                                                                                                                    >  1:133146/1‑150 (MQ=255)
  ccATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTaaaagaaa                                                                                                                                                   >  1:1183839/1‑149 (MQ=255)
         gcCGGAGGGATCGCAATAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGtcattc                                                                                                                                           >  1:970977/1‑150 (MQ=255)
            ggAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATtcatca                                                                                                                                       <  2:734014/151‑1 (MQ=255)
             gAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTcatcat                                                                                                                                      <  1:187946/151‑1 (MQ=255)
              atggATGGCATTAGGTGCCACTAGTATTATTTATCCGGCTGATTCTAACCAGACTGCGGTGTGGATTAGATATTGCCATACCAATGGGCGAATTCTGGTCAATTCGTGGATTGGAACCATCGTCGGTGTAAAATAAATGCCATgcatca                                                                                                                                       <  1:748315/147‑1 (MQ=255)
                  aTCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTacac                                                                                                                                 <  2:790226/151‑1 (MQ=255)
                    cGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACAcc                                                                                                                                <  2:578709/150‑1 (MQ=255)
                    cGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCg                                                                                                                               <  2:539664/151‑1 (MQ=255)
                      cATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGcc                                                                                                                             <  2:610088/151‑1 (MQ=255)
                                cACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACtgtttgtt                                                                                                                   >  1:944922/1‑151 (MQ=255)
                                   gtgTATTATTTACCCCGCTGATGCTAAACAGACGGCGGTATGGATTAAAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGGGGTGAAAAAAAAAAGTCATTCATCATTCCACCGCCACTGTTTGTTAGt                                                                                                                <  1:1061412/149‑1 (MQ=255)
                                    cgTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTg                                                                                                               >  2:281961/1‑151 (MQ=255)
                                     gTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTg                                                                                                               >  2:760822/1‑150 (MQ=255)
                                     gTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTg                                                                                                               >  2:393954/1‑150 (MQ=255)
                                      ttttgtttatCCCTCTGATTCTAAACACACTGCGTTTTGTATTAAAAATAGCCATACCAATGAGCGCTTTTTGTTCAGTTGGTGGATTAAAAACAGCCGGGGTGAAAAAGAAAAGCCCTCCATCATCACACCGCCACTGTTTGTTAGt                                                                                                                <  1:1096502/143‑1 (MQ=255)
                                      tattatTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGa                                                                                                              >  1:63223/1‑150 (MQ=255)
                                        ttttttacCCCGCTGATGCAAAACAGCCGGGGTTATGGATTAGAAATACCCATACCAATGAGCGTTTTCTGGTCAATTCGTGAATTAAAAACAGCAGCGGTTAAAAAAAAAATTCATTCATCATTACACCGCCACtgtttgtt                                                                                                                   <  1:957479/140‑1 (MQ=255)
                                        ttatttacCCCGCTGATGCTAAACGGCCTGCGGTATGGATAAAAAATAGCCATACAAAGGAGCGTTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATCCATCTTTACACCGCCACTGTTTGTTAGTg                                                                                                               <  1:42742/147‑1 (MQ=255)
                                         tatttatCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAAccc                                                                                                          >  1:625237/1‑151 (MQ=255)
                                              atCCCGCTGATGCAAAAAAAATTGCGGTATGAATTAGAAATAGCCATCCCAATGAGCGTTTTCTGGTAAATTCGTGTATTAAAAACGGCAGCGTTATAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCaaaa                                                                                                      <  1:570784/150‑1 (MQ=255)
                                                cccGCTGATGCAAAACATACGGCGGTATGAATAAAAAATAGCCATACAAATGAGCGCTTTCTGGTCATTTCGTGGATTAAAAACAGCAGGGGTGAAAAAAAAAAGTCATTCATAATTACACCGCCACTGTTTGTTAGTGAACCCaaaa                                                                                                      <  1:634343/148‑1 (MQ=255)
                                                cccGCTGAGGCAAAACAGACTGCGGTATGGATTAGAAATAGCCATACAAATGAGCGTTTTCTGGTCAATTCGTGGATTAAAAACAGCAGCGGTGAAAAAAAAAAGTAATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCaaaa                                                                                                      <  1:338582/148‑1 (MQ=255)
                                                  cGCTGATGCTAAACAGACTGCGGTATGAATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGaa                                                                                                 <  2:133392/151‑1 (MQ=255)
                                                         gaaaaaCAGACTGCGTTATGGATAAAAAATAGCCATCCAAATGAGCTCTTCCGGGTAATTTCGTGTATTGAACACAGCAGCGTGGAAAAAAAAAAGAAATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACtt                                                                                          <  1:259265/148‑1 (MQ=255)
                                                         acaaaacagctggCGGTATGGATTAAAAATCGCCATACAAAGGAGCTCTTTCTGGTCAATTCGTGGTTTGAAAACGGCGGGGGTGAAAAAAAAAAATCTTTCTTCATTACCCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACt                                                                                           <  1:699635/138‑1 (MQ=255)
                                                            aaagagatgggtgTATGATTTAAAAATAGCCATACCAATGGGCGTTTTCTGTTCATTTGGTGTATAAAAAACAGCAGCGATAAAAAAAAAAATTCATCCATAATTCCGCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGc                                                                                       <  1:138512/139‑1 (MQ=255)
                                                            aaaCAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGAAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACt                                                                                           <  1:1198546/147‑1 (MQ=255)
                                                            aaaCAGACGGGGGTATGGATTAAAAATAACCATACAAATGAGCGTTTTCGGGTAATTTCGTGGATTAAAAAGAGAAGCGGTTTAAAAGAAAAGTAATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACtt                                                                                          <  2:840219/148‑1 (MQ=255)
                                                              aCAGACTGCGGTGTGGATTAGAAATCGCCATACCAATGAGTGCTTTTTGGTCATTTCGTGGATTAAAAACAGCAGCGGTTTAAAAAAAAAGTCATTCTTCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTg                                                                                        <  2:898516/148‑1 (MQ=255)
                                                              aCAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGt                                                                                     >  1:1037620/1‑151 (MQ=255)
                                                                   ctgggttttGGATTAGAAATAGCAATACAAATGAGCACTTTCGGGTCATTTCGTGGATTAAAACCAGCGGGGTTGAAAAAAAAAATAAATAAACCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGtat                                                                                   <  1:146194/140‑1 (MQ=255)
                                                                   cggCGGTATGGATTAGAAATAGCCATACCCATGAGCGCTTTCTGGTCAATTCGTGGATTAAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGtatta                                                                                 <  2:670700/148‑1 (MQ=255)
                                                                     gggttATGGATTAAAAATAGCCATACCAATGAGCCCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTATAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGta                                                                                    <  1:843609/141‑1 (MQ=255)
                                                                     gCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTattatt                                                                               <  2:906455/150‑1 (MQ=255)
                                                                     gCGGTATGAATTAGAAATAGCCATACAAATGAGCGCTTTCTGGTAAATTCGTGGATTAAAAACAGCAGCGGTGAAAAAAAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGtattat                                                                                <  1:1108343/149‑1 (MQ=255)
                                                                         tATGGATTAGAAATAGCCATGCAAATGAGCTCTTCCTTGTAAATTGGTGTATTGAAAACAGCAGGGGTGAAAAAAAAAGGTCATTCATCATAACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTaca                                                                           <  2:1036734/150‑1 (MQ=255)
                                                                          aTGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACAcc                                                                         >  2:765498/1‑151 (MQ=255)
                                                                            ggATTAGAAATAGCCATCCCAATGAGCGTTTTCTGGTCAATTCGTGGATTGAAAACACCAGCGGTGTAAAAAAAAATTCATTCATAATAACCCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGtattat                                                                                <  2:405651/142‑1 (MQ=255)
                                                                                tAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTcc                                                                    >  1:989152/1‑150 (MQ=255)
                                                                                tAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTcc                                                                    >  2:1141376/1‑150 (MQ=255)
                                                                                 aaaaaTAGCAATACCAATGAGCGTTTTCTGGTCATTTCGTGGATTGAAAACAGCGGGGGTTAAAAAAAAAATTCATCCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTcc                                                                    <  1:582279/147‑1 (MQ=255)
                                                                                          cATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGgc                                                           >  1:674515/1‑149 (MQ=255)
                                                                                             aCCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGgcagc                                                        >  2:235870/1‑149 (MQ=255)
                                                                                              ccAATGAGCGTTTTCTGTTAATTTGGTGTTTTAAAACCAACAGGGGTAAAAAAAAAAGTTCTTCCACCTTTTCTCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCt                                                              <  1:992827/142‑1 (MQ=255)
                                                                                                 aTGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGgcagc                                                        >  2:1198708/1‑145 (MQ=255)
                                                                                                  ttaGCGCTTTCGGTTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATCCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGcagcag                                                      >  1:81918/3‑146 (MQ=255)
                                                                                                  tGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGt                                                 >  2:79758/1‑151 (MQ=255)
                                                                                                   gagctttttCTGGTCATTTCTTGGATTGAAAACAGCAGCGGTGAAAAAGAAAAGTCATCCTTCTTTCCACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTg                                                <  1:717964/145‑1 (MQ=255)
                                                                                                   gagctttttCTGGTAAATTCGTGGATTGAAAACAGCAGCGGTGAAAAAAAAAAGTCTTTCATCATTACGCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCg                                                  <  2:848795/143‑1 (MQ=255)
                                                                                                   gAGCGCTTTCTGGTCAATTCGTGGATTAAAAACAGCAGCGGTATAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGt                                                 <  1:476905/150‑1 (MQ=255)
                                                                                                     gcgcTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATc                                                   >  2:301879/1‑146 (MQ=255)
                                                                                                      cgcTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGt                                             >  1:586741/1‑151 (MQ=255)
                                                                                                              gttCAATTCGTGAATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTctct                                         <  1:351672/145‑1 (MQ=255)
                                                                                                               gtaaaTTGGTGAATTGAAAACAGCGGGGGTGAAAAAAAAAATCCATCCATCTTTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGtt                                      <  1:970628/146‑1 (MQ=255)
                                                                                                                  atttCGTGGATTAAAACCAGCGGCGGTGAAAAAAAAAAGTCTTTCACCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTgg                                  <  1:752529/148‑1 (MQ=255)
                                                                                                                    ttCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATg                               >  2:822424/1‑151 (MQ=255)
                                                                                                                       gTGGATTAAAAACAGCGGCGGTAAAAAAAAAAAGTCATCCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAAt                            <  1:661810/151‑1 (MQ=255)
                                                                                                                        gaggtggAAAGCAGCAGGGGTGTAAAGGAAAAATGCTTAGATATCACGCCGCCACTGGTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTAGTTACACCGTTCCGCCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGaa                             <  1:1107015/143‑1 (MQ=255)
                                                                                                                            ttaaaaaCAGCAGCGGTTAAAAAAAAAAGTCATCAATCATCACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATg                           <  2:323322/144‑1 (MQ=255)
                                                                                                                             tGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTa                        >  1:1032712/1‑149 (MQ=255)
                                                                                                                               aaaaCAGCAGGGGTGAAAAAAAAAAGTCATCCAGCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGAc                    <  1:966146/151‑1 (MQ=255)
                                                                                                                               aaaaCAGAAGGGGTGAAAAAAAAAAGTCATTATAAATCACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGa                     <  1:724975/150‑1 (MQ=255)
                                                                                                                                  aCAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGAt                 >  2:161438/1‑151 (MQ=255)
                                                                                                                                   cagcagCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATc                >  1:504685/1‑151 (MQ=255)
                                                                                                                                       agcGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCtt            >  2:7123/1‑151 (MQ=255)
                                                                                                                                        gcGGTGAAAAAGAAAATTCATCCATCATTCCACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGa                  <  1:466887/144‑1 (MQ=255)
                                                                                                                                         ctttttcccagaaaagTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCg          >  2:525775/10‑151 (MQ=255)
                                                                                                                                          ggTATAAAAAAAAAGTAATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCg          <  1:232242/150‑1 (MQ=255)
                                                                                                                                          ggTATAAAAAAAAACAAATCCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCg          <  1:77807/150‑1 (MQ=255)
                                                                                                                                           gtgtAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGt        <  2:657042/151‑1 (MQ=255)
                                                                                                                                               aaaaaaaaagCCATCCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGAt    <  1:618350/151‑1 (MQ=255)
                                                                                                                                                 aagaaaagTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATaa  <  2:49222/151‑1 (MQ=255)
                                                                                                                                                  agaaaagTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATaa  >  2:1181211/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACCAATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGCAGCGGTGTAAAAGAAAAGTCATTCATCATTACACCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAA  >  NC_000913/559034‑559329

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