Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,535,7850AG45.1% 97.4 / 22.1 71F25S (TTC→TCC) yddEPhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base A (7/32);  new base G (32/0);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 9.94e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTG  >  NC_000913/1535635‑1535933
                                                                                                                                                       |                                                                                                                                                    
gcgTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGa                                                                                                                                                      <  1:326510/151‑1 (MQ=255)
          aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACAcc                                                                                                                                            <  1:907349/151‑1 (MQ=255)
                    ttCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCGCAGAAtt                                                                                                                                   >  1:279120/1‑150 (MQ=255)
                      cAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCGCAggaat                                                                                                                                   >  1:749396/1‑144 (MQ=255)
                      cAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCAGaaattc                                                                                                                                 >  1:24026/1‑146 (MQ=255)
                       aaCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAaacaaaa                                                                                                                                              >  1:332575/1‑134 (MQ=255)
                       aaCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCGccggaa                                                                                                                                    >  1:617846/1‑141 (MQ=255)
                        aCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCcgc                                                                                                                              <  1:561832/151‑1 (MQ=255)
                           gTTGGCGTAAAGTGGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGGAAAAAAAACCGGcgaaatggcgc                                                                                                                               >  2:187860/1‑138 (MQ=255)
                           gTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACACGCAGGATAGTcgcg                                                                                                                             >  1:300081/1‑149 (MQ=255)
                           gTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGaa                                                                                                                           <  1:33914/151‑1 (MQ=255)
                             tGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCCCAGCATTGCCGCCaaa                                                                                                                          >  1:107850/1‑150 (MQ=255)
                              ggCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACg                                                                                                                        <  1:444450/151‑1 (MQ=255)
                                   aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGGCGGAATTGcgcggaaaggg                                                                                                                       >  1:25003/1‑138 (MQ=255)
                                     aGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGa                                                                                                                 <  1:376341/151‑1 (MQ=255)
                                     aGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGa                                                                                                                 <  2:1023655/151‑1 (MQ=255)
                                        aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAAATGCCGCGGAACGGGTGTGAg                                                                                                                >  2:357103/1‑149 (MQ=255)
                                        aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGt                                                                                                              <  1:703883/151‑1 (MQ=255)
                                        aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGt                                                                                                              <  1:499255/151‑1 (MQ=255)
                                           ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAAATGTCGCGGAACGGGtgggggg                                                                                                               >  2:85232/1‑142 (MQ=255)
                                            gATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTaaaa                                                                                                          <  2:1019616/151‑1 (MQ=255)
                                            gATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTaaaa                                                                                                          <  2:464611/151‑1 (MQ=255)
                                             aTGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCGGCCGAATTTCCGCGGAACGGGtggggggggaaa                                                                                                          >  2:574367/1‑140 (MQ=255)
                                             aTGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAg                                                                                                         <  2:329539/151‑1 (MQ=255)
                                              tGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGAAATGCCGCGGAAAGGGtgggggggaga                                                                                                           >  1:668124/1‑139 (MQ=255)
                                               gcgcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGc                                                                                                        <  1:666035/150‑1 (MQ=255)
                                                   acaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAATTGGCGCGGCACGGGTGTGGGGTgaaa                                                                                                          >  1:169995/1‑144 (MQ=255)
                                                    caTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGAc                                                                                                  <  2:771664/151‑1 (MQ=255)
                                                    caTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACGCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGaa                                                                                                  <  1:172489/151‑2 (MQ=255)
                                                        gCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGAAATTCCGCGGAAAGGGTGTGGGGGGAAAAGATCGACgggg                                                                                              >  1:1005832/1‑151 (MQ=255)
                                                          tGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGGAAAAAAACCCGCAGAATTGCCGCGCAAAGGGTGTGGGGGAAAAAAAAAGACGTgg                                                                                              >  1:60385/1‑149 (MQ=255)
                                                          tGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACAACCGCAGAATTTGCGCGAAACGGGTGTGAGGGAAAAAAAACGACgggg                                                                                              >  2:1015887/1‑149 (MQ=255)
                                                          tGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACACGCCGAAATGGCGGGAAAAGGGTGGGAGGGAAAAgaaaagcagggg                                                                                              >  1:745391/1‑138 (MQ=255)
                                                            tcgtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCGCCGAAATGCCGCGGAACGGGTGTGGGGGAAAAACAACCAAGTGGGac                                                                                           >  2:694195/1‑150 (MQ=255)
                                                                cgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTg                                                                                      <  2:1010837/151‑1 (MQ=255)
                                                                 gCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGc                                                                                     <  1:363305/151‑1 (MQ=255)
                                                                        agcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCAGGATTGCCGCGGAAAGGGTGTGGGGGAAAAAAATCGACGTGGTGCAACTg                                                                                      >  1:800083/1‑143 (MQ=255)
                                                                             aaaaGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCAGCAGAAATGCCGCGCAAAGGGTGTGGGGGAAAAGAATAGACGGGGGAGACCTACggggttatt                                                                            >  1:823712/1‑141 (MQ=255)
                                                                             aaaaGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAAc                                                                         <  2:327563/151‑1 (MQ=255)
                                                                                  cAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCGGCAGGATTGCCGCGGAACGGGTGTGGGGGAAAAAAAACGAAGGGGGACACCCGCGGTTTCATTAAAGttt                                                                     >  2:457313/1‑149 (MQ=255)
                                                                                        ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCAGAATTGTCCCGAAACGGGTGTGGGGGGAAAACAACCAAGGGGGAAACCCGCGGTTTTAGTTACGtt                                                                      >  1:385771/1‑143 (MQ=255)
                                                                                        ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCtgt                                                               <  1:654126/150‑1 (MQ=255)
                                                                                          gAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACAACCGCAGGATTGCCCCGAAAACGGTGTGGGGGAAAAAGATCCACGTGGGAAACCCGCGGGTTCATTAAAGTTCTCCTGTGAc                                                            >  2:918/1‑151 (MQ=255)
                                                                                           aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCAGAAATGCCGCGGAACGGGTGTGAGGGAAAAACAACGACGGGGGACAACCGCGGGTTTATTAAACTTCTcc                                                                  >  2:838926/1‑144 (MQ=255)
                                                                                           aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAATTGCCGCGGAAAAGGTGGGGGGGAAAAAAAACGACGAGGGAGACCCGCGGGTTCATTAACGTTCTCCTGTGAc                                                            >  1:142276/1‑150 (MQ=255)
                                                                                              ggCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCGCAGAAATGCCGCGAAACGGGTGTGGGGGAAAAAAATCGACGGGGGACACCCGCGGGTTTATTTAAGTTCTCCTGTGGCTgg                                                         >  2:640708/1‑150 (MQ=255)
                                                                                                ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACAACGGCAGAATTGCCCCGAAACGGGTGGGCGGGGAAAACCACGACGGGGGACGCCCGCGGGTTCTTTAACGGTCTCCTGTGACCgggg                                                       >  1:324107/1‑148 (MQ=255)
                                                                                                ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGgaga                                                      <  2:748149/151‑1 (MQ=255)
                                                                                                ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGgaga                                                      <  2:914818/151‑1 (MQ=255)
                                                                                                    aCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCGCAGAAATGCCGCGAAAAGGGTGTGAGGGGAAAAAATCGACGGGGGACACCCGCGGGTTTATTAAAGTTCTCCCGTTAATggggga                                                     >  2:451559/1‑144 (MQ=255)
                                                                                                    aCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCAGCAGGATTGCCGCGGAAAGGGTGTGAGGGGAAAACATCGACGTGGGACACCCGCGGGTTTATTAAAGTTTTCCTGTTAATGGGGaaat                                                   >  1:969048/1‑148 (MQ=255)
                                                                                                         cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATa                                             <  1:966970/151‑1 (MQ=255)
                                                                                                         cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAAAACCGGCAGAAATGCCGCGAAACGGGTGTGGGGGAAAAGCATCCACGTGGGACAACTGCGGGTTCATTAAAGTTCTCCCGTGGCTGGGGAAATATCAt                                              >  2:49656/1‑150 (MQ=255)
                                                                                                            gCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGcc                                          <  2:502927/151‑1 (MQ=255)
                                                                                                             cGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAAAACCCGCAGGAATGGCGCGCAAAGGGTGTGGGGGGAAAAAAACGACGGGGTACACCCGCGGGTTTATTAACCTTCTCCCTGGACTGGAGAAAAATCATAAcc                                          >  1:413834/1‑150 (MQ=255)
                                                                                                                 aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCAGAATTGCCGCGGAACGGGTGTGGGGGAAAAACAACCACGGGGGACACCCGCGGGTTTATTTAAGTTTTCCTCTGACTGGGGAACTATCATAAcc                                          >  2:638108/1‑146 (MQ=255)
                                                                                                                 aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCa                                      <  1:85225/150‑1 (MQ=255)
                                                                                                                  aGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAg                                    <  2:533490/151‑1 (MQ=255)
                                                                                                                   gCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGGATTGCCGCGAAAAGGGTGTGGGGGAAAAAAAACGACGGGGGACACCCGCGGGTTTATTTAAGGTCTCCTTTGGCTTGGGAACTATCATAAcc                                          >  1:487947/1‑144 (MQ=255)
                                                                                                                     tGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGAATTGCCGCGGAACGGGTGTGGGGGAAAAACATCGACGGGGGAAACCTGCGGGTTTATTAAAGTTCTCCCGTGGCTGGGGAAATATCATAGCCTGCAAGTgg                                 >  2:199024/1‑151 (MQ=255)
                                                                                                                        aTTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACAACGGCCGAATTGTCGCGGAACAGGTGTGGGGGAAAAGCAACCACGCGGGACACCCCCGGGTTTATTAAATATATCCTGTGACTGGGGAAATATCCTAGCCTGCAAGTg                                  >  1:106367/1‑147 (MQ=255)
                                                                                                                             cgcTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGagag                         <  2:871629/151‑1 (MQ=255)
                                                                                                                              gcTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGc                        <  2:5452/151‑1 (MQ=255)
                                                                                                                                 tCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGa                      <  1:881229/150‑1 (MQ=255)
                                                                                                                                  cGCTAAGATTATCGGC‑GGGGAAAAAAACACCCGCCGAAATGCCGCGAAAAGGGTGTGGGGGAAAAAAAACCAACTGGGAAACCCGCGGTTTCATTAAACTTCTCCTGTGACTGGGGAACAATCATAACCTGCaaaggggcgggg                           >  2:407264/1‑134 (MQ=255)
                                                                                                                                            aTCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCg           <  2:60391/150‑1 (MQ=255)
                                                                                                                                            aTCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCg           <  1:482524/150‑1 (MQ=255)
                                                                                                                                               ggc‑ggGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGGACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGcc        >  2:135980/1‑150 (MQ=255)
                                                                                                                                                 c‑ggGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg     <  2:1005998/151‑1 (MQ=255)
                                                                                                                                                   ggggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAggg    >  1:813919/1‑151 (MQ=255)
                                                                                                                                                    gggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt   <  1:1019144/151‑1 (MQ=255)
                                                                                                                                                     ggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg  <  2:617491/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTG  >  NC_000913/1535635‑1535933

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