Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,465,1670GT39.4% 108.6 / 32.9 71C29F (TGT→TTT) gspMgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (9/34);  new base T (28/0);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 2.72e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAA  >  NC_000913/3465023‑3465314
                                                                                                                                                   |                                                                                                                                                    
aTGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCt                                                                                                                                                <  1:138618/151‑1 (MQ=255)
 tGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCCCGGCGGTTTTTTTTGGGtt                                                                                                                                               >  1:1154409/1‑151 (MQ=255)
  ggAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTg                                                                                                                                              <  2:330365/151‑1 (MQ=255)
       aaaCAGAAAAAGATCAGGGCAGGGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGATCTCTGTCGCGGGATTTTTTGTTgggctgg                                                                                                                                             >  2:1178365/1‑142 (MQ=255)
         aCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCCCGGCGGTTTTTTTTGGGTTgggg                                                                                                                                           >  2:378707/1‑147 (MQ=255)
               aaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTTCCCCGGCGGTTTTTTTTGGGTTGGGGGGGTaaaaac                                                                                                                                 >  1:1008882/1‑149 (MQ=255)
                  gATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCGCCCGGCGGTTTTTTTTGGGTTGGGGGGGTaaaa                                                                                                                                   >  1:678485/1‑146 (MQ=255)
                      cGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTATCCTGCCTGGGGGTTTTTTTTGGGTTGGGGGGaaaaaaac                                                                                                                                 >  1:66864/1‑142 (MQ=255)
                       ggTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTCTTGTTGGCTTGGGGTGATAAAACCgatagat                                                                                                                          <  2:1019924/150‑1 (MQ=255)
                          cAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATAcc                                                                                                                       <  2:738543/150‑1 (MQ=255)
                            gtggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCCTCCTCCGGGGTTTTTTTTTGGTTTGGGGGGTAAAAACgaatgataac                                                                                                                       >  2:1167997/1‑146 (MQ=255)
                              ggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACa                                                                                                                   <  1:402970/150‑1 (MQ=255)
                              ggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACAt                                                                                                                  <  1:808138/151‑1 (MQ=255)
                              ggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCCCGGCGTTTTTTTTTTGCTTGGGGGGGTAAAACCCATAGATACCTACa                                                                                                                   >  1:642611/1‑150 (MQ=255)
                                tgtATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGGTTTTTTTTTGCTTGGGGGGGTAAAAACCATAGATACCTACAtt                                                                                                                 >  1:1132998/1‑150 (MQ=255)
                                   aTGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTTTCCCTCGGCGGTTTTTTTTGGGTTTGGGGGATAAAAACgatagataaccacctt                                                                                                                 >  2:1431517/1‑139 (MQ=255)
                                    tGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTTTTGGCTTGGGGGGATAAAAACGATAAATACCTACAtt                                                                                                                 >  1:202550/1‑146 (MQ=255)
                                     ggACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCagag                                                                                                           <  1:1066151/151‑1 (MQ=255)
                                             aGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTc                                                                                                  <  2:250545/151‑1 (MQ=255)
                                               gAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCGCTCCGCGTTTTTTGTTGGGTTGGGGGGGTAAAAACGATAGATACCTACAATTCCGAGGATCAATccc                                                                                                >  2:839621/1‑151 (MQ=255)
                                                 aCTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCCCCGCGTTTTTTTTTGGCTTGGGGGGATAAAAACGATAAATAACCACATTGCAGAGCATCaagccc                                                                                                 >  1:889439/1‑144 (MQ=255)
                                                    cAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCAcgc                                                                                            <  1:254107/150‑1 (MQ=255)
                                                     aGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCCCTCGGCGTTTTTTTTTTGCTTGGGGGGAAAAAAACGGAAAATACCCATAtttgaaaagatga                                                                                                      >  1:1128075/1‑129 (MQ=255)
                                                     aGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACgcgc                                                                                          <  1:26595/151‑1 (MQ=255)
                                                       ggAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCGGCGGTTTTTTTTTGCTTTGGGGGGTAAAACCCATAGAAACCTACATTTCTGAGCATCAGTCCCAcgc                                                                                            >  2:1339683/1‑147 (MQ=255)
                                                         aaaTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCaaaa                                                                                      <  1:1288611/151‑1 (MQ=255)
                                                             gATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAAtt                                                                                  <  1:1057193/151‑1 (MQ=255)
                                                              aTTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTTTTTGGTTTGGGGGGTAAAAACGATAGATACCCACATTTCAGAGCATCAGTCCCACCCGCAAAAAAtt                                                                                  >  1:1191917/1‑150 (MQ=255)
                                                              aTTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTATTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTa                                                                                 <  2:489426/151‑1 (MQ=255)
                                                                 aaaTCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTaag                                                                               <  1:656435/150‑1 (MQ=255)
                                                                     cATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAgg                                                                           <  2:474480/150‑1 (MQ=255)
                                                                     cATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCCCGGCGGTTTTTTTTGGGTTGGGGGGGTAAAACCGATAGATACCTACATTTCAGAGCATCAGTCCCCCGCGCAAAAAAATaaaaa                                                                             >  1:1019710/1‑148 (MQ=255)
                                                                      aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGAc                                                                         <  1:717877/151‑1 (MQ=255)
                                                                      aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGAc                                                                         <  1:325526/151‑1 (MQ=255)
                                                                      aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCCGGGGTTTTTTTTGGGTTTGGGGGGTAAAAACGATTGATACCTTCATTTGAGAGCAACAGTCCCCCGCGCAAAAAATTaaaaa                                                                             >  1:1292216/1‑147 (MQ=255)
                                                                         gtggGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCGGCGGTTTTTTTTTGGTTTGGGGGGTAAAACCCCAAGAAACCCACATTTCAGAGCATCAGTCCCACCCGCAAAAAAATaaaaaaggaa                                                                        >  2:1412415/1‑144 (MQ=255)
                                                                               aGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCGGCGGTTTTTTTTGGGTTGGGGGGATAAAACCGAAAGATAACTACATTTCCGAGCATCAGGCCCCCGCGCAAAAAATAAAAAAAGACATCAAAAgg                                                                >  2:1293963/1‑151 (MQ=255)
                                                                                   aaaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGc                                                            <  1:135094/151‑1 (MQ=255)
                                                                                       ccACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCGGGGGTTTTTGTTTGGTTGGGGGGGTAAAAACGATAGATTTCTACAATTCAGAAGATCAGGCCCACGCGCCAAAAAATAAAAAAGACATCAAATTGATGCAgg                                                         >  2:263053/1‑150 (MQ=255)
                                                                                        cacaTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAt                                                       <  1:757740/151‑1 (MQ=255)
                                                                                        cacaTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCCCCGCGTTTTTTTTTGGCTTGGGGGGGTAAAAACGATAGATACCTACCTTTTAGAGCATCAGTCCCCCCCGCAAAAAAATAAAAAGGACATCAAATGGATGCAgg                                                         >  2:325030/1‑149 (MQ=255)
                                                                                             cGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCGCTCGGGGGTTTTTTTTGGGTTTGGGGGATAAAAACGATAGATACCTACATTGCAGAGCATCAGTCCCCCCCGCAAAAAAATAAAAAGGACATCAAATAGATGCAGGATCaa                                                    >  2:335136/1‑149 (MQ=255)
                                                                                             cGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGc                                                  <  2:26788/151‑1 (MQ=255)
                                                                                              ggAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCg                                                 <  2:299756/151‑1 (MQ=255)
                                                                                              ggAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCCGGCGTTTTTGTTTGGTTTGGGGGGTAAAAACGGGAGAAAACCACACTTCAGAGCATCAGGCCCCCCGCCCAAAAAATAAAAAAGACATCAAAATGATGCAggggtaa                                                    >  1:128425/1‑143 (MQ=255)
                                                                                               gAAAAACAGATCGTTGCGGCACTGGCGCTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATGGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGa                                                <  2:935877/151‑1 (MQ=255)
                                                                                                 aaaaCAGATCGTTGCGGCACTGGCGGTTCTCTCTCCTGGCGGTTTTTTTTGGGTTGGGGGGATGAAACCGCTAGACACCTTTATTTCTGAGGATCTGTCCCACGCGCGAAAAATTAAAAAGGACAAAAAAATGATGCCggggcaa                                                    >  1:907773/1‑140 (MQ=255)
                                                                                                  aaaCAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCa                                             <  1:808731/151‑1 (MQ=255)
                                                                                                    aCAGATCGTTGCGGCACTGGCGGTTTTCTTGCTCGGGGCTTTTTTTTGGCTTGGGGGGGGTAAAACGCAAGATACCTACATTGCTGAGCCTTCGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGAATGAGGGTTAAGCGAggcc                                            >  2:625332/1‑147 (MQ=255)
                                                                                                       gATCGTTGCGGCACTGGCGGTTCTTTCGCTCGGCGGTTTTTTTTGGCTTGGGGGGGTAAAAACGATAGATACCCACCTTTCAGAGCATCAGGCCCCCACGCAAAAAATTAAAAAAGAAATTAAAAGGATGCAGGATCAAGCGAGCAccc                                          >  1:1102791/1‑149 (MQ=255)
                                                                                                        aTCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTTTTGGGTTGGGGGGATAAAAACGATAGATACCTACATTGCAGAGGATCAGGCCCACGCGCAAAAAAATAAGAAGGACATAAAAAGGATGCAGGATCAAACGAGCACCCATg                                       >  1:155874/1‑151 (MQ=255)
                                                                                                        aTCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATg                                       <  2:757563/151‑1 (MQ=255)
                                                                                                          cGTTGCGGCACTGGCGGTTCTCTCGCCCGGCGTTTTTTTTTGGGTTGGGGTGATTAAAAAGATAGAAACCTACATTGCAGAGCATCAGTCCCCCCCGCAAAAAAATAAGAAGGACATCAAAAGGATGCAGGATCAAGCGAGCACCCCTgg                                      >  2:118532/1‑150 (MQ=255)
                                                                                                            ttGCGGCACTGGCGGTTCTCTCGCTCGGCGGTTTTTTTTGGGTTGGGGGGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAAATAAAAAGGACATCAAAAGGATGCAGGATCAAACGAGCACCCATGGgtt                                   >  1:380313/1‑151 (MQ=255)
                                                                                                              gCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTa                                  <  1:848917/150‑1 (MQ=255)
                                                                                                              gCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTa                                  <  1:565789/150‑1 (MQ=255)
                                                                                                              gCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTAc                                 <  2:1197234/151‑1 (MQ=255)
                                                                                                               cGGCACTAGCAGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGAGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACt                                <  2:500677/151‑1 (MQ=255)
                                                                                                                  cACTGGCGGTTCTCTCGCTCGGCGTTTTTTTTTGGCTTGGGGGGGTAAAAACGATAGATACCTACATTGTAGAGCATCAGTCCCCCGCGCAAAAAATTAAAAAGGACATCAAATTGATGCAGGAACAAGCGAGCAACCATGGGTTACTCgg                             >  2:526636/1‑151 (MQ=255)
                                                                                                                      ggCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTc                           <  2:678284/149‑1 (MQ=255)
                                                                                                                         ggTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCt                      >  2:981823/1‑151 (MQ=255)
                                                                                                                         ggTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCt                      >  1:997302/1‑151 (MQ=255)
                                                                                                                            tctctcGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTgc                    <  1:824358/150‑1 (MQ=255)
                                                                                                                                ctcgctcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCtt                >  2:597570/1‑151 (MQ=255)
                                                                                                                                  cgctcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCtt                <  2:889171/149‑1 (MQ=255)
                                                                                                                                  cgctcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTa               <  1:187001/150‑1 (MQ=255)
                                                                                                                                  cgctcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTAc              <  2:1154060/151‑1 (MQ=255)
                                                                                                                                   gctcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACg             >  1:801894/1‑151 (MQ=255)
                                                                                                                                     tcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGc            >  1:408163/1‑150 (MQ=255)
                                                                                                                                     tcgGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGc            >  2:900933/1‑150 (MQ=255)
                                                                                                                                            gTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCa      <  1:1168354/151‑1 (MQ=255)
                                                                                                                                                 ttgttgGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCaa  >  1:185229/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAA  >  NC_000913/3465023‑3465314

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