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

ATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAA  >  NC_000913/3465023‑3465314
                                                                                                                                                   |                                                                                                                                                    
aTGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCt                                                                                                                                                <  1:138618/151‑1 (MQ=255)
 tGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCC‑CGGC‑GGTTTTTTTTGGGtt                                                                                                                                               >  1:1154409/1‑151 (MQ=255)
  ggAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTg                                                                                                                                              <  2:330365/151‑1 (MQ=255)
       aaaCAGAAAAAGATCAGGGCAGGGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGATCT‑CTGTCG‑CGGG‑ATTTTTTGTTgggctgg                                                                                                                                             >  2:1178365/1‑142 (MQ=255)
         aCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCC‑CGGC‑GGTTTTTTTTGGGTTgggg                                                                                                                                           >  2:378707/1‑147 (MQ=255)
               aaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTTCCC‑CGGC‑GGTTTTTTTTGGGTTGGGGGGGTaaaaac                                                                                                                                 >  1:1008882/1‑149 (MQ=255)
                  gATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTT‑CTCGCC‑CGGC‑GGTTTTTTTTGGGTTGGGGGGGTaaaa                                                                                                                                   >  1:678485/1‑146 (MQ=255)
                      cGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTAT‑CCTGCC‑TGGG‑GGTTTTTTTTGGGTTGGGGGGaaaaaaac                                                                                                                                 >  1:66864/1‑142 (MQ=255)
                       ggTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTCTTGTTGGCTTGGGGTGATAAAACCgatagat                                                                                                                          <  2:1019924/150‑1 (MQ=255)
                          cAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATAcc                                                                                                                       <  2:738543/150‑1 (MQ=255)
                            gtggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCCTCCTC‑CGGG‑GTTTTTTTTTGGTTTGGGG‑GGTAAAAACgaatgataac                                                                                                                       >  2:1167997/1‑146 (MQ=255)
                              ggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACa                                                                                                                   <  1:402970/150‑1 (MQ=255)
                              ggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACAt                                                                                                                  <  1:808138/151‑1 (MQ=255)
                              ggtgTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCC‑CGGC‑GTTTTTTTTTTGCTTGGGGGGGTAAAACCCATAGATACCTACa                                                                                                                   >  1:642611/1‑150 (MQ=255)
                                tgtATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GGTTTTTTTTTGCTTGGGGGGGTAAAAACCATAGATACCTACAtt                                                                                                                 >  1:1132998/1‑150 (MQ=255)
                                   aTGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTT‑TTCCCT‑CGGC‑GGTTTTTTTTGGGTTTGGGGGATAAAAACgatagataaccacctt                                                                                                                 >  2:1431517/1‑139 (MQ=255)
                                    tGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTTTTGGCTTGGGGGGATAAAAACGATAAATACCTACAtt                                                                                                                 >  1:202550/1‑146 (MQ=255)
                                     ggACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCagag                                                                                                           <  1:1066151/151‑1 (MQ=255)
                                             aGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTc                                                                                                  <  2:250545/151‑1 (MQ=255)
                                               gAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTT‑CTCGCT‑CCGC‑GTTTTTTGTTGGGTTGGGGGGGTAAAAACGATAGATACCTACAATTCCGAGGATCA‑ATccc                                                                                                >  2:839621/1‑151 (MQ=255)
                                                 aCTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCC‑CCGC‑GTTTTTTTTTGGCTTGGGGGGATAAAAACGATAAATAACCACATTGCAGAGCATCaagccc                                                                                                 >  1:889439/1‑144 (MQ=255)
                                                    cAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCAcgc                                                                                            <  1:254107/150‑1 (MQ=255)
                                                     aGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTT‑CTCCCT‑CGGC‑GTTTTTTTTTTGCTTGGGGGGAAAAAAACGGAAAATACCCATAtttgaaaagatga                                                                                                      >  1:1128075/1‑129 (MQ=255)
                                                     aGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACgcgc                                                                                          <  1:26595/151‑1 (MQ=255)
                                                       ggAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCT‑CGGC‑GGTTTTTTTTTGCTTTGGGGGGTAAAACCCATAGAAACCTACATTTCTGAGCATCA‑GTCCCAcgc                                                                                            >  2:1339683/1‑147 (MQ=255)
                                                         aaaTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCaaaa                                                                                      <  1:1288611/151‑1 (MQ=255)
                                                             gATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAAtt                                                                                  <  1:1057193/151‑1 (MQ=255)
                                                              aTTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTTTTTGGTTTGGGGGGTAAAAACGATAGATACCCACATTTCAGAGCATCA‑GTCCCACCCGCAAAAAAtt                                                                                  >  1:1191917/1‑150 (MQ=255)
                                                              aTTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTATTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTa                                                                                 <  2:489426/151‑1 (MQ=255)
                                                                 aaaTCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTaag                                                                               <  1:656435/150‑1 (MQ=255)
                                                                     cATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAgg                                                                           <  2:474480/150‑1 (MQ=255)
                                                                     cATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCC‑CGGC‑GGTTTTTTTTGGGTTGGGGGGGTAAAACCGATAGATACCTACATTTCAGAGCATCA‑GTCCCCCGCGCAAAAAAATaaaaa                                                                             >  1:1019710/1‑148 (MQ=255)
                                                                      aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGAc                                                                         <  1:717877/151‑1 (MQ=255)
                                                                      aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGAc                                                                         <  1:325526/151‑1 (MQ=255)
                                                                      aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCT‑CCGG‑GGTTTTTTTTGGGTTTGGGGGGTAAAAACGATTGATACCTTCATTTGAGAGCAACA‑GTCCCCCGCGCAAAAAATTaaaaa                                                                             >  1:1292216/1‑147 (MQ=255)
                                                                         gtggGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCT‑CGGCGGTTTTTTTTTGGTTTGGGG‑GGTAAAACCCCAAGAAACCCACATTTCAGAGCATCA‑GTCCCACCCGCAAAAAAATaaaaaaggaa                                                                        >  2:1412415/1‑144 (MQ=255)
                                                                               aGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCT‑CGGC‑GGTTTTTTTTGGGTTGGGGGGATAAAACCGAAAGATAACTACATTTCCGAGCATCA‑GGCCCCCGCGCAAAAAATAAAAAAAGACATCAAAAgg                                                                >  2:1293963/1‑151 (MQ=255)
                                                                                   aaaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGc                                                            <  1:135094/151‑1 (MQ=255)
                                                                                       ccACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCT‑CGGG‑GGTTTTTGTTTGGTTGGGGGGGTAAAAACGATAGATTTCTACAATTCAGAAGATCA‑GGCCCACGCGCCAAAAAATAAAAAAGACATCAAATTGATGCAgg                                                         >  2:263053/1‑150 (MQ=255)
                                                                                        cacaTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAt                                                       <  1:757740/151‑1 (MQ=255)
                                                                                        cacaTCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCC‑CCGC‑GTTTTTTTTTGGCTTGGGGGGGTAAAAACGATAGATACCTACCTTTTAGAGCATCA‑GTCCCCCCCGCAAAAAAATAAAAAGGACATCAAATGGATGCAgg                                                         >  2:325030/1‑149 (MQ=255)
                                                                                             cGGAAAAACAGATCGTTGCGGCACTGGCGGTTTT‑CTCGCT‑CGGG‑GGTTTTTTTTGGGTTTGGGGGATAAAAACGATAGATACCTACATTGCAGAGCATCA‑GTCCCCCCCGCAAAAAAATAAAAAGGACATCAAATAGATGCAGGATCaa                                                    >  2:335136/1‑149 (MQ=255)
                                                                                             cGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGc                                                  <  2:26788/151‑1 (MQ=255)
                                                                                              ggAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCg                                                 <  2:299756/151‑1 (MQ=255)
                                                                                              ggAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCCCTCCGGC‑GTTTTTGTTTGGTTTGGGG‑GGTAAAAACGGGAGAAAACCACACTTCAGAGCATCAGGCCCCCCGCCCAAAAAA‑TAAAAAAGACATCAAAATGATGCAggggtaa                                                    >  1:128425/1‑143 (MQ=255)
                                                                                               gAAAAACAGATCGTTGCGGCACTGGCGCTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATGGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGa                                                <  2:935877/151‑1 (MQ=255)
                                                                                                 aaaaCAGATCGTTGCGGCACTGGCGGTTCT‑CTCTCC‑TGGC‑GGTTTTTTTTGGGTTGGGGGGATGAAACCGCTAGACACCTTTATTTCTGAGGATCT‑GTCCCACGCGCGAAAAATTAAAAAGGACAAAAAAATGATGCCggggcaa                                                    >  1:907773/1‑140 (MQ=255)
                                                                                                  aaaCAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCa                                             <  1:808731/151‑1 (MQ=255)
                                                                                                    aCAGATCGTTGCGGCACTGGCGGTTTT‑CTTGCT‑CGGG‑GCTTTTTTTTGGCTTGGGGGGGGTAAAACGCAAGATACCTACATTGCTGAGCCTTC‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGAATGAGGGTTAAGCGAggcc                                            >  2:625332/1‑147 (MQ=255)
                                                                                                       gATCGTTGCGGCACTGGCGGTTCT‑TTCGCT‑CGGC‑GGTTTTTTTTGGCTTGGGGGGGTAAAAACGATAGATACCCACCTTTCAGAGCATCA‑GGCCCCCACGCAAAAAATTAAAAAAGAAATTAAAAGGATGCAGGATCAAGCGAGCAccc                                          >  1:1102791/1‑149 (MQ=255)
                                                                                                        aTCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTTTTGGGTTGGGGGGATAAAAACGATAGATACCTACATTGCAGAGGATCA‑GGCCCACGCGCAAAAAAATAAGAAGGACATAAAAAGGATGCAGGATCAAACGAGCACCCATg                                       >  1:155874/1‑151 (MQ=255)
                                                                                                        aTCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATg                                       <  2:757563/151‑1 (MQ=255)
                                                                                                          cGTTGCGGCACTGGCGGTTCT‑CTCGCC‑CGGC‑GTTTTTTTTTGGGTTGGGGTGATTAAAAAGATAGAAACCTACATTGCAGAGCATCA‑GTCCCCCCCGCAAAAAAATAAGAAGGACATCAAAAGGATGCAGGATCAAGCGAGCACCCCTgg                                      >  2:118532/1‑150 (MQ=255)
                                                                                                            ttGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GGTTTTTTTTGGGTTGGGGGGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAAATAAAAAGGACATCAAAAGGATGCAGGATCAAACGAGCACCCATGGgtt                                   >  1:380313/1‑151 (MQ=255)
                                                                                                              gCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTa                                  <  1:848917/150‑1 (MQ=255)
                                                                                                              gCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTa                                  <  1:565789/150‑1 (MQ=255)
                                                                                                              gCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTAc                                 <  2:1197234/151‑1 (MQ=255)
                                                                                                               cGGCACTAGCAGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGAGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACt                                <  2:500677/151‑1 (MQ=255)
                                                                                                                  cACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTTTTGGCTTGGGGGGGTAAAAACGATAGATACCTACATTGTAGAGCATCA‑GTCCCCCGCGCAAAAAATTAAAAAGGACATCAAATTGATGCAGGAACAAGCGAGCAACCATGGGTTACTCgg                             >  2:526636/1‑151 (MQ=255)
                                                                                                                      ggCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTc                           <  2:678284/149‑1 (MQ=255)
                                                                                                                         ggTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCt                      >  2:981823/1‑151 (MQ=255)
                                                                                                                         ggTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCt                      >  1:997302/1‑151 (MQ=255)
                                                                                                                            tct‑ctcGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTgc                    <  1:824358/150‑1 (MQ=255)
                                                                                                                                ctcgct‑cgGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCtt                >  2:597570/1‑151 (MQ=255)
                                                                                                                                  cgct‑cgGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCtt                <  2:889171/149‑1 (MQ=255)
                                                                                                                                  cgct‑cgGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTa               <  1:187001/150‑1 (MQ=255)
                                                                                                                                  cgct‑cgGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTAc              <  2:1154060/151‑1 (MQ=255)
                                                                                                                                   gct‑cgGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACg             >  1:801894/1‑151 (MQ=255)
                                                                                                                                     t‑cgGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGc            >  1:408163/1‑150 (MQ=255)
                                                                                                                                     t‑cgGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGc            >  2:900933/1‑150 (MQ=255)
                                                                                                                                            gTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCa      <  1:1168354/151‑1 (MQ=255)
                                                                                                                                                 ttgttgGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCaa  >  1:185229/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCT‑CTCGCT‑CGGC‑GTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCA‑GTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAA  >  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: ‑