Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1090AC41.1% 161.9 / 30.6 92G219G (GGT→GGG) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (44/9);  new base C (0/37);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 4.51e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATT  >  NC_000913/4537960‑4538252
                                                                                                                                                          |                                                                                                                                                   
aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCAc                                                                                                                                                    >  1:642998/1‑151 (MQ=255)
 tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCAcc                                                                                                                                                   >  1:1380427/1‑151 (MQ=255)
    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑tt‑t                                                                                                                                              <  1:1396790/151‑1 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑tt‑t                                                                                                                                              >  1:942688/1‑150 (MQ=255)
     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑Ta                                                                                                                                             <  1:1127560/151‑1 (MQ=255)
         aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑c                                                                                                                                         >  1:1133685/1‑150 (MQ=255)
         aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑c                                                                                                                                         >  1:407290/1‑150 (MQ=255)
           aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑caca                                                                                                                                      <  2:800340/151‑1 (MQ=255)
            ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACaa                                                                                                                                     <  2:1203552/151‑1 (MQ=255)
             cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑c                                                                                                                                   >  2:292496/1‑151 (MQ=255)
               ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTTGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑ccg                                                                                                                                 <  1:1137766/151‑1 (MQ=255)
               ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑ccg                                                                                                                                 <  2:214822/151‑1 (MQ=255)
               ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑ccg                                                                                                                                 >  2:399533/1‑151 (MQ=255)
                 aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑ccgcc                                                                                                                               <  1:69304/151‑1 (MQ=255)
                           tCTGGTGATGAGACGGAAGCAAGTTTATTCCATTTAAAATTATTCCCGGTGAAATCAAGTTAAAATACGCCCTCCTTCCGCAAT‑CCTGTTTTGGC‑TCCGCC‑AT‑TAAAAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGcc                                                                                                                      <  1:801052/150‑1 (MQ=255)
                             tGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGtt                                                                                                                   >  2:293831/1‑151 (MQ=255)
                                     agaCGGGAGAAATTTTATCCTATTTTAAATTATCACCGGTAAAATAAAGTTAAAATCGGGCATCCGTTCCCACT‑CCTGTTTTGGC‑TCCGCC‑AT‑AAAAACGCCGGG‑TTTTACCCCT‑TT‑TATT‑CACAC‑CCGCCGCACCAGCCGTTCCGTAc                                                                                                             <  1:365315/149‑1 (MQ=255)
                                           gAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                      >  2:1275461/1‑150 (MQ=255)
                                               atttttttcccttttaAAATTATCCCCGGGAAAATCAATTACAAATCCGGCACCCTTTCCCAAC‑CGTGGTTT‑GC‑CTCCCC‑AA‑TAAAAAACCCGG‑TTTTTTCCCCTTT‑TATT‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                 <  1:1389226/139‑1 (MQ=255)
                                                 gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                               >  2:422012/1‑151 (MQ=255)
                                                 gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                               >  2:720058/1‑151 (MQ=255)
                                                  cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                               >  2:821585/1‑150 (MQ=255)
                                                  cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                               >  1:954502/1‑150 (MQ=255)
                                                     aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGAACCACCCGTTCCGTACAAGGGCGTTTCGCCACCg                                                                                           >  1:273093/1‑151 (MQ=255)
                                                      ttCCATTTTAATTTATTCCCGGTAAAATAAAGTTAAAATACGCCTCCCGTTCGCAAT‑CCTGTTTTGGC‑TTCGCA‑TT‑AAAAAAGCCAGT‑TTTTATCCCC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                           <  2:89785/150‑1 (MQ=255)
                                                          aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                      >  2:361424/1‑151 (MQ=255)
                                                          aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCAATTCGCCAGCGTAACt                                                                                      >  2:989436/1‑151 (MQ=255)
                                                           ttttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc                                                                                     >  2:955957/1‑151 (MQ=255)
                                                            tttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                    >  2:421198/1‑151 (MQ=255)
                                                               aaaTTATTCCGGGGAAAATAAATTTAAAATGCGGCATCCGTTCGAAAC‑CCTGTTTTGGT‑TCCGCC‑AT‑AAATAAGCCGGT‑TTTTATCCCT‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                 <  2:784903/151‑1 (MQ=255)
                                                               aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                 >  1:983876/1‑151 (MQ=255)
                                                               aaaTTATCCCGGGTAAAAACAATTAAAAATCCGGC‑TCCCTTCCCAATCCCGGGTTTGCC‑TCCGCA‑AT‑AAATAACCCAGT‑TTTTA‑CCCC‑TT‑TTTT‑CCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                  <  1:123164/150‑1 (MQ=255)
                                                                    attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                            >  1:622046/1‑151 (MQ=255)
                                                                         cGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                        >  2:474203/1‑150 (MQ=255)
                                                                         cGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                        >  1:892840/1‑150 (MQ=255)
                                                                          ggaaaaaTCAAGTTAAACTCCGCCTCCCTTCCGCATT‑CTTGGTTTGGC‑TTCCCA‑AT‑AAAAACGCAAGT‑TTTTATCCCC‑TT‑TTTC‑CCCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                        <  1:1268160/145‑1 (MQ=255)
                                                                          gataaaaTCAGGTTCAAATCCGCCACCCGTTCAAACT‑CTTGG‑TTGTT‑TTCCCC‑AA‑TAAAAAACCAAG‑TTTTTCCCCC‑TT‑TTTTCCCCAT‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                               <  1:204733/138‑1 (MQ=255)
                                                                           gTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                     >  1:1129558/1‑151 (MQ=255)
                                                                            tGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                    >  2:832809/1‑151 (MQ=255)
                                                                             gAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑ATTT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                    >  1:98817/1‑150 (MQ=255)
                                                                              aatacaattaaaaaTACGGCCCCCTTCCGC‑AC‑CCTTGTTTTGG‑TTCCCC‑AT‑TAAAAAACCCGGTTTTTACCCCT‑TT‑TATC‑CCCAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                    <  1:1409798/138‑1 (MQ=255)
                                                                              aaaTCAATTTAAAATCCGGCACCCGTTCGCATT‑CCTGTTTTGGC‑TCCCCC‑AT‑TAATAAGCCAGT‑TTTTATCCCT‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                     <  2:1246886/148‑1 (MQ=255)
                                                                              aaaTCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                  >  1:1117627/1‑151 (MQ=255)
                                                                                 tCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                >  2:792681/1‑150 (MQ=255)
                                                                                   aactaaaaaTACGACACCCTTCCCAAAT‑CCTGGTTTGCT‑TTCCCA‑AT‑AAATAAGCATGT‑TTTTACCCCC‑TT‑TATC‑AACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                              <  2:492474/144‑1 (MQ=255)
                                                                                      ttCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                          >  2:972232/1‑151 (MQ=255)
                                                                                       tCAAATACGGCATCCGTTCGCAAT‑CCTGGCTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                          >  2:1382857/1‑150 (MQ=255)
                                                                                        aaaaTCCGCCACCCGTCCG‑AAT‑CCCGGTTTGGG‑TTCCCC‑AT‑TAAAAACGCAGGTTTTTACCCCT‑TT‑TATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                        <  2:1289634/150‑1 (MQ=255)
                                                                                         aaaTACGGCATCCTTTCGCAAT‑CCTGTTTTGGC‑TTCCCC‑AT‑AAAAAAGCCGGG‑TTTTATCCCT‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                        <  2:1017699/150‑1 (MQ=255)
                                                                                         aaaTACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                       >  2:1267314/1‑151 (MQ=255)
                                                                                         aaaTACGCCACCCGTCCGAAAC‑CCTGTTTTGGT‑TCCCCC‑AT‑AAATAGCCCAGT‑TTTTATCCCC‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                           <  2:1265121/147‑1 (MQ=255)
                                                                                          aaaaccgccacctttCCCAAT‑CCTTGGTTTGC‑TTCCCC‑CA‑TAAAAACCCAGGTTTTTACCCCC‑TT‑TATC‑CCCAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                          <  1:1122995/135‑1 (MQ=255)
                                                                                             aCGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                    >  1:598498/1‑150 (MQ=255)
                                                                                               ggcccccGTCCGCAAC‑CCTGTTTTGGT‑TCCCCC‑AA‑AAATACCCCAGG‑TTTTCTCCCT‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                 <  1:1275222/146‑1 (MQ=255)
                                                                                               gccATCCTTCCCCAAT‑CCTGTTTTGGC‑TCCCCA‑AT‑TAAAAAGCCAGT‑TTTTACCCCT‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                     <  2:259407/145‑1 (MQ=255)
                                                                                                gCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                >  1:263036/1‑151 (MQ=255)
                                                                                                 cATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                               >  1:1028508/1‑151 (MQ=255)
                                                                                                    ccttccgaaaT‑CCTGTTTTGGC‑TTCGCC‑TT‑TAATAAGCAGGG‑TTTTATCCCC‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTa                                              <  2:110651/141‑1 (MQ=255)
                                                                                                    ccttccGCAAC‑CCT‑GTTTGGC‑TTCCCCAAT‑AAATACCCCAGT‑TTTTACCCCT‑TT‑TATC‑CCCAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                             <  1:227248/145‑1 (MQ=255)
                                                                                                    ccGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                            >  2:925124/1‑151 (MQ=255)
                                                                                                    ccGTTCGCAAT‑CCTGGTTTGGC‑TCCGCT‑AT‑AAAAAAGCCAGT‑TTTTATCCCC‑TT‑TATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                 <  1:1169428/146‑1 (MQ=255)
                                                                                                     cGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                           >  2:330772/1‑151 (MQ=255)
                                                                                                       ttCGCAAC‑CCTGTTTTGCC‑TCCGCC‑AT‑TAAAAAGCCAGT‑TTTCACCCCC‑TT‑TATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                          <  1:525878/150‑1 (MQ=255)
                                                                                                         cccact‑CCTGTTTTGCC‑TCCGCA‑AT‑AAAAACCCCAGT‑TTTTACCCCT‑TT‑TATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                          <  1:361370/143‑1 (MQ=255)
                                                                                                            aaT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                    >  2:898921/1‑151 (MQ=255)
                                                                                                             atccgggttttGGC‑TCCCCC‑TT‑TAATCA‑CCAGG‑TTTTAACCCC‑TTGTATT‑AACAG‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                     <  2:851472/142‑1 (MQ=255)
                                                                                                                ccTGTTTTGTC‑TTCTCT‑AT‑TAAAACCCTAGG‑TTTTTTCATC‑TT‑AATC‑CACGC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                   <  2:1076849/149‑1 (MQ=255)
                                                                                                                 ccggTTTGGCTTTCCCC‑TT‑TAATAAGCCAGT‑TTTTCCCCCT‑TT‑TATA‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                    <  2:511352/146‑1 (MQ=255)
                                                                                                                 cTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                >  1:105826/1‑151 (MQ=255)
                                                                                                                    gTTTGGC‑CTCCCC‑ATGTAAAATGCCGGT‑TTTTATCCTT‑TT‑TATC‑CAAAC‑CCGCCGCACCAGCCGTGCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                  <  1:1379440/147‑1 (MQ=255)
                                                                                                                     tttGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                            >  1:1236190/1‑151 (MQ=255)
                                                                                                                     tttGGC‑TCCGCC‑AT‑TAAAAACCCAGT‑TTTTTTC‑CC‑TT‑TTTT‑CCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                 <  1:783781/146‑1 (MQ=255)
                                                                                                                       tttgttttccctaT‑TAAAAAGCCTGT‑TTTTACCCCT‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                             <  2:622153/138‑1 (MQ=255)
                                                                                                                       tGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                          >  1:1133955/1‑151 (MQ=255)
                                                                                                                         gC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                        >  1:137536/1‑151 (MQ=255)
                                                                                                                         gC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                        >  1:624518/1‑151 (MQ=255)
                                                                                                                            ttCACC‑AT‑TAATAAGCGTTT‑TTTTATCCCT‑TT‑TATA‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                       <  2:842352/150‑1 (MQ=255)
                                                                                                                              cGCA‑AT‑AAAAAAGCCAGT‑TTTTATCCCC‑TT‑TATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                    <  2:1161329/151‑1 (MQ=255)
                                                                                                                                cc‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                  >  1:1179288/1‑151 (MQ=255)
                                                                                                                                cc‑AT‑TAAAAAGCCAGG‑TTTTTTCCCC‑TT‑TTTT‑CCCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                  <  2:273140/151‑1 (MQ=255)
                                                                                                                                     taaaaaaGCCAGT‑TTTTATCCCT‑TT‑TATC‑AAAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                  <  1:712081/143‑1 (MQ=255)
                                                                                                                                       aataaGCCAGT‑TTTTACCCCT‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                      <  2:302604/142‑1 (MQ=255)
                                                                                                                                       aataaGCCAGG‑TTTTACCCCC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa             <  2:964301/151‑1 (MQ=255)
                                                                                                                                       aaaaacgccggtttttACCCCT‑TT‑TATC‑CAAAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                         <  1:619926/128‑1 (MQ=255)
                                                                                                                                       aaaaaccaaGG‑TTTTACCCCT‑TT‑TATC‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa               <  2:430205/141‑1 (MQ=255)
                                                                                                                                        aaaaacccggtttttCTCCCC‑TT‑TTTC‑CCCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                       <  1:1181754/130‑1 (MQ=255)
                                                                                                                                              ctttcttttACCCCC‑TT‑TATC‑CAACT‑CCTCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACGCCATTGCTGTGTTGAGGGCTCAAAAGACAACAGAAACTTATTCAAGAAATAATCTTCTGCtttt                      <  1:202141/131‑1 (MQ=255)
                                                                                                                                              ccggtttttACCCCC‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt         <  2:1281904/144‑1 (MQ=255)
                                                                                                                                                 gtttttATC‑CC‑TT‑TATTCCACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc     <  2:98830/148‑1 (MQ=255)
                                                                                                                                                 g‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt   >  1:462694/1‑151 (MQ=255)
                                                                                                                                                  tttttATCCCT‑TT‑TATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc     <  2:1133892/148‑1 (MQ=255)
                                                                                                                                                   ttttATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt  >  1:887014/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                                   
ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAAT‑CCTGGTTTGGC‑TTCGCC‑AT‑TAATAAGCCAGG‑TTTTATCACC‑TT‑TATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATT  >  NC_000913/4537960‑4538252

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