Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,9000AG40.3% 133.8 / 17.5 78noncoding (72/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (8/38);  new base G (31/0);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 2.31e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925751‑926048
                                                                                                                                                        |                                                                                                                                                           
gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGAt                                                                                                                                                            <  1:1228836/151‑1 (MQ=255)
  gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTc                                                                                                                                                          <  1:70474/151‑1 (MQ=255)
         gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑GA‑GGGGGGGGTTCGCAcccc                                                                                                                                                   >  1:1174684/1‑151 (MQ=255)
           cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAAccccca                                                                                                                                                 <  2:871894/151‑2 (MQ=255)
           cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                 <  1:252153/151‑1 (MQ=255)
           cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                 <  2:167642/151‑1 (MQ=255)
           cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                 <  2:458743/151‑1 (MQ=255)
            gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGT‑TA‑GGGGGGG‑TTTGAAcccccc                                                                                                                                                 >  1:848230/1‑149 (MQ=38)
              aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATa                                                                                                                                              <  1:535478/151‑1 (MQ=255)
                     gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTATATATGGC‑GGG‑GA‑AGGGGGGAATTTAAccccccgcttgcttcc                                                                                                                                       >  1:6372/1‑140 (MQ=37)
                        tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑GA‑GGGGGGGAGTTGTAcccccccaaccgttgcg                                                                                                                                      >  1:1238771/1‑137 (MQ=37)
                          gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTCGAACCCCCCCTACG‑TTGCCGttt                                                                                                                                  >  2:781764/1‑148 (MQ=255)
                            aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑aca                                                                                                                              <  2:1173465/151‑1 (MQ=255)
                            aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑aca                                                                                                                              <  2:282449/151‑1 (MQ=255)
                                cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGGACCCCCCACACG‑TTGCCGTAT‑acac‑a                                                                                                                           >  2:50954/1‑149 (MQ=25)
                                  caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑Cttt                                                                                                                      <  2:615061/151‑1 (MQ=255)
                                    ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTcc                                                                                                                    <  2:269127/151‑1 (MQ=255)
                                       aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAgg                                                                                                                 <  2:4665/151‑1 (MQ=255)
                                       aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAgg                                                                                                                 <  2:811856/151‑1 (MQ=255)
                                            aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GG‑GGGGGGGGGTTGAACCCCCCAAACC‑CTTCCCTAT‑ACAC‑ACCTTTCCAg                                                                                                                  >  2:379838/1‑146 (MQ=21)
                                             cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGc                                                                                                            <  2:530494/150‑1 (MQ=255)
                                               ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTCGAAACCCCCCTAAG‑TTGTCGTGT‑ACAA‑A‑CCTTTCACGCGGGCTcc                                                                                                         >  1:494476/1‑151 (MQ=21)
                                                      cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGc                                                                                                   <  1:344557/150‑1 (MQ=255)
                                                      cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑GA‑GGGGGGGGGTCGAACCCCCCGTACC‑TTGCCGGAT‑ACAC‑A‑CTTTCCATGCGGGGTCCCTCCGcc                                                                                                  >  2:433086/1‑151 (MQ=18)
                                                       ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTGGAAACCCCCCGATG‑TTGCCGTAT‑TAAC‑A‑ATTTCCACGCGGGCttc                                                                                                         >  2:849157/1‑141 (MQ=18)
                                                        caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                                <  2:1204694/151‑1 (MQ=255)
                                                         accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGAACCCCCCGTACG‑ATGCCGCAT‑ACAC‑A‑CCTTTCTTGCGGGCTCCTTCcaccac                                                                                                >  1:111762/1‑150 (MQ=18)
                                                           caTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCg                                                                                             <  2:1090080/151‑1 (MQ=255)
                                                             tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTCGAACCCCCCATACG‑TTGCCGTAT‑AAAC‑A‑CCTTCCAGGCGGGCTCCTTCAGCCAccccg                                                                                            >  1:663014/1‑146 (MQ=18)
                                                             tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑TA‑GGGGGGGGGTTTAACCCCCCACAAG‑TTTCCGGAT‑ACAC‑A‑GTTTTCAAACGTGCCCCTTcccccgcccagg                                                                                           >  1:801741/1‑140 (MQ=16)
                                                              gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                          <  1:620516/151‑1 (MQ=255)
                                                                tCACTGCCGTACAAACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                         <  1:1152603/150‑1 (MQ=255)
                                                                   cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCt                                                                                     <  2:235128/150‑1 (MQ=255)
                                                                    tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGAAACCCCCATACG‑CTGCCGTAT‑ACAC‑A‑CTTTTCAGGCGGGGGCCCTCAGCCCCACGGAGA‑cccc                                                                                    >  2:184592/1‑148 (MQ=17)
                                                                          acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GAGGGGGGGGGTTTAAACCCCCCA‑ACG‑TTGCCCGAT‑ACAC‑A‑CTTTCCAAGGGTGCTCCCTCCACCCCCCGGAAA‑CCCCA‑CCaaa                                                                             >  1:459452/1‑150 (MQ=14)
                                                                          acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCaaa                                                                             <  1:1422937/150‑1 (MQ=255)
                                                                            agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTg                                                                          <  1:685775/151‑1 (MQ=255)
                                                                            agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑GG‑GGGGGGGGGTCTAAACCCCCCTACG‑TTGCCCTAT‑AAAC‑A‑ACTTCCAGGGGGGCTGCCTCTGCCACCCGGGCA‑CCCCC‑CCCAAtac                                                                          >  2:273699/1‑149 (MQ=9)
                                                                               cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGttt                                                                       <  2:572619/151‑1 (MQ=255)
                                                                                  aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTCGGACCCCCCCTAAG‑TTTCCGGCT‑AAAC‑A‑CTTTCCCAGCGGGCTCCTTCCGCCACCCGGGCA‑CCCCA‑CCCAAAtttttt                                                                      >  1:275121/1‑149 (MQ=14)
                                                                                     aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTAc                                                                 <  1:31949/151‑1 (MQ=255)
                                                                                      aaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGAACCCCCCCCACG‑ATTCCGTAT‑AAAC‑A‑ATTTTCAGGCGGGGTCCTTCATCCCCCCCGACA‑CCCCC‑CCAAAATGTTTTTCt                                                                   >  1:950888/1‑148 (MQ=9)
                                                                                      aaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑GA‑AGGGGGGGGTTGAAACCCCGATACG‑TTTCCCTAT‑ACAC‑A‑CCTTCCCGGGGGGGTCCTTTAGCCCCCCCGACA‑CCCCA‑CCCAAATGTTTTTCTAcc                                                                >  2:360765/1‑151 (MQ=9)
                                                                                       aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑GA‑GGGGGGGGTTTGGAAACCCCCTAAG‑CTGGCGTAT‑ACAA‑A‑CTTTTCAGGGGGGCTCCTTTCAACACCCGGACA‑CCCCC‑CCAAAATTTTTTTCTatc                                                                >  1:1261593/1‑148 (MQ=2)
                                                                                         gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TGGTCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGTCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCa                                                               <  1:902409/149‑1 (MQ=255)
                                                                                         gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑GA‑GGGGGGGGTTTGAACCCCCCCCACG‑TTGCCGTAT‑ATAA‑A‑CTTTTCACGCGGGCCCCCTCAGCCACCCGGGCA‑CCCCA‑CCCAAATAttttttttcc                                                                >  2:97725/1‑142 (MQ=2)
                                                                                          cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GG‑GGGGGGGGTTCGAACCCCCGCTACG‑TTGCCGGAT‑AAAC‑A‑CTTTCCAGGCGCGCCCCCTCACCCCCCCGGGCA‑CCTCACCCAAATttttttttacccaa                                                              >  1:1141226/1‑142 (MQ=11)
                                                                                           gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAAcc                                                           <  1:369280/151‑1 (MQ=255)
                                                                                            aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑TA‑GGGGGGGAATCGGACCCCCCCTACG‑TTGCCGGAT‑TAAC‑A‑CCTTTCCGGCGGGCCCCCTCCGCCCCTCCGACA‑CCCCC‑CCAAATTTTTTTTCTTCCAAAccc                                                          >  1:1328054/1‑150 (MQ=9)
                                                                                             aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GG‑GGGGGGGGGTTGGACCCCCCCTACG‑TTGCCGGTT‑AAACAA‑CTTTCC‑GGCGGGGTCCCTTTGCCACCCCGACA‑CCCCACCCAAATATTTTTTTTACCAAAc                                                            >  2:1383954/1‑149 (MQ=2)
                                                                                              aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTTGCACCCCCCATACG‑CTGCCCTAT‑ACAC‑C‑CCTTTCCGGCGTGCTCCTTCAGCCACTCGGCCA‑CCCCA‑CCCAATTTTTTTTGTACCcaaac                                                           >  1:473322/1‑146 (MQ=11)
                                                                                               aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGA‑GA‑AGGGGGGGGTCGAAACCCCCCAATG‑TTGCCCTTT‑AACC‑A‑ATTTCTCGGGGGGCTTCTCCTGCCCCACGCACACCCCCA‑CCAAATttttttctcccacacccc                                                          >  2:955226/1‑137 (MQ=2)
                                                                                                 aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GG‑GGGGGGGGGTTGAACCCCCCCCACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGGGCTCCCTCAGCCACCCCGACA‑CCCCA‑CCCAAATGTTTTTTTACCAAAACTCCT‑gg                                                    >  2:571660/1‑151 (MQ=14)
                                                                                                  aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGATTGAAACCCCCAAACC‑TTGCCCGATAAAAC‑A‑TTTTCCA‑GCGCGGTCCTTCCGCCCCCCGGGCACCCCCA‑CCAAATT‑TTTTTCTAACAAACCCCAT‑ggg                                                   >  1:1276919/1‑151 (MQ=9)
                                                                                                   aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGt                                                  <  1:1220307/151‑1 (MQ=255)
                                                                                                        gCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGTGGG‑GGGGGGGGTTGGAACCCCC‑CTAAG‑TTTCCGTAT‑ACAC‑A‑CTTTTCCAGCGGGGTCCCTCCGCCCCCCCGGCA‑CCCCC‑CCCAAATGTTTTTGTTCCAAACCCCAT‑ggggggg                                               >  1:461317/1‑148 (MQ=2)
                                                                                                        gCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGGTCGGACCCCCCATACG‑TTGCCGTAT‑ACAC‑A‑CATTCCAGGGGGGCTCCTTCATCCACCCGGACA‑CCCCC‑CCAAAATGTTTTTCTTCCAAAACCCAT‑gggggg                                                >  1:91099/1‑148 (MQ=11)
                                                                                                           cTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GG‑GGGGGGGGTTTGAAACCCCCATACG‑TTTGCGTAT‑ACAA‑A‑CCTTTCCGGGGGGGTCCTTCCGCCACCCCGACA‑ACCCC‑CCCAAATGTTTTTTTACCAAAACCCAT‑TGGGGGGAAccg                                          >  1:1172657/1‑149 (MQ=2)
                                                                                                           cTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGG‑GAGGGGGGGGGTTGGAACCCCC‑CCACG‑TTGCCGTAT‑AAAC‑A‑CTTTTCCGGCGGGCTCCCTCCGCCCCTCGGACA‑CCCCC‑CCAAAATTTTTTTCTACCCAAACCCAT‑GGGTgggga                                             >  1:869644/1‑145 (MQ=2)
                                                                                                               cGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGCACCCCCCCTACG‑GTTCCGGAT‑ACAC‑A‑CTTTCCCGGCGGGCTCCTTTAGCCACTCGGACA‑CCCCA‑CCAAAATGTTTTTCTACCAAAACTCATGGGGGGGCAACgggcg                                       >  1:1401267/1‑149 (MQ=14)
                                                                                                                 tACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGGTTTGGACCCCCCCGACGTTTGCCG‑AT‑ACAA‑C‑CTTTTCCGGCGGGCTCCCTCAGCCACCCCGGCA‑CCCCC‑CCAAAATGTTTTTTTAACAAAACCCAC‑ggggggggaaaggg                                        >  1:147624/1‑139 (MQ=2)
                                                                                                                 tACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGgcgc                                     <  1:1074055/150‑1 (MQ=255)
                                                                                                                   cGAGCTCTTCTTTAAATATGGC‑GGG‑GC‑GGGGGGGGCTCTAACCCCCCCAACG‑TTTACGCAT‑ACACAA‑CTTTCCA‑GCGTGCTCCTTTTCCCACACGGGGA‑CCTCA‑CCAAAAAGCTTTTCTACCCAACCTCAT‑TGGGGGaaaaaggggggt                                    >  1:17781/1‑140 (MQ=2)
                                                                                                                    gAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGctact                                 <  1:466007/151‑1 (MQ=255)
                                                                                                                      gCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACtat                               <  1:50954/151‑1 (MQ=255)
                                                                                                                         cttcttTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACtata                              <  2:1113192/149‑1 (MQ=255)
                                                                                                                          ttcttTAAATATGGC‑GGT‑GA‑GGGGGGGAGTCGAAACCCCCATACG‑TTGTCGTAT‑ACAC‑A‑ATTTCCAGGCGTGCTCCTTCCGCCACCCGGGCC‑CCCCA‑CCCAAATGTTTTTCTATCAAAACTCAT‑GGGGGGGAACCGGGGGCTACTATAggg                           >  1:539065/1‑151 (MQ=17)
                                                                                                                           tcttTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAgg                            <  1:55397/149‑1 (MQ=255)
                                                                                                                                    atGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                   >  1:113316/1‑149 (MQ=255)
                                                                                                                                    atGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                  <  1:247080/150‑1 (MQ=255)
                                                                                                                                       gC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa              <  1:952058/151‑1 (MQ=255)
                                                                                                                                        c‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc             <  2:120239/151‑1 (MQ=255)
                                                                                                                                           gT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg           <  2:958431/151‑1 (MQ=255)
                                                                                                                                               a‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc         >  1:46042/1‑150 (MQ=255)
                                                                                                                                               a‑GGGGGGGATTAGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGGAAAACGGTc         >  2:209283/1‑150 (MQ=255)
                                                                                                                                                  ggggggATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa        >  2:1008577/1‑149 (MQ=255)
                                                                                                                                                    ggggATTCGAACCCCCGATACG‑TTGCCGTAG‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa    <  1:819274/151‑1 (MQ=255)
                                                                                                                                                      ggATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga  <  1:1150196/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                           
GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGC‑GGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGA  >  NC_000913/925751‑926048

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