Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC41.0% 165.0 / 24.7 96Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (49/7);  new base C (0/39);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 7.51e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCC  >  NC_000913/3584937‑3585231
                                                                                                                                                       |                                                                                                                                                 
ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTAc                                                                                                                                                  >  2:461545/1‑151 (MQ=255)
     gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCccgc                                                                                                                                             >  1:696166/1‑151 (MQ=255)
     gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCccg                                                                                                                                              <  2:1323741/150‑1 (MQ=255)
        aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGt                                                                                                                                          >  1:1122372/1‑151 (MQ=255)
         aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGtt                                                                                                                                         >  2:507064/1‑151 (MQ=255)
           gACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTa                                                                                                                                        >  1:269188/1‑150 (MQ=255)
             cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAat                                                                                                                                     >  2:1252310/1‑151 (MQ=255)
              aGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAatc                                                                                                                                    >  1:361394/1‑151 (MQ=255)
                       attcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCg                                                                                                                            >  2:848001/1‑150 (MQ=255)
                        ttctttctCGGGAAACGAAAAAACCGCGGAAAAAATATCTAAAGGATTACAAAAAAAAAAGGGCTC‑CCTCCGTTTTCCCCTGAAAAAATCCCCCTCTCCCCCGGAAGAGGGCCGTTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGa                                                                                                                           <  1:239495/150‑1 (MQ=255)
                        ttcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCg                                                                                                                            >  2:549890/1‑149 (MQ=255)
                            ttcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                       >  1:1368623/1‑150 (MQ=255)
                             tctcTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                     >  2:983776/1‑151 (MQ=255)
                               tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                                                   >  2:98847/1‑151 (MQ=255)
                                  gggAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                >  2:638939/1‑151 (MQ=255)
                                  gggAACGGAAAACCCGCGGGAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTC‑CCTTCGGCTTGCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTAACCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                <  2:247565/151‑1 (MQ=255)
                                    gAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                              >  1:1377897/1‑151 (MQ=255)
                                       cGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                           >  2:200818/1‑151 (MQ=255)
                                        gcaaaacCCGCGGAAAAAATATCTACAGTATTACAAAAAAAAAAAGCCTC‑CCTCCGGCTTCCCCGGAAAAAATCTCCCTCTCCCCACAAAGAGGCCCGCTATCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                            <  1:283569/146‑1 (MQ=255)
                                        aacaagacccgagaaagaaTATCTACAGAATAAAAAAAAAAAAATGTCCC‑CCTTCGGCCAGCCCTAACAAATTACCCCTCACCCCAAGAAAAGGGCCGCTATCCTGG‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                           <  1:839962/136‑1 (MQ=255)
                                          cAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCg                                                                                                        >  2:1485484/1‑151 (MQ=255)
                                           aaacCCGCTGAAAAATTTTCTACAGAATTACAAAAAAAAAAAGGCTA‑CCTTCGTCTTCCCCTAAAAAAATACCCCTCTCCCCCCGAGGAGGCCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                            <  1:648349/146‑1 (MQ=255)
                                           aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                       >  2:1188246/1‑151 (MQ=255)
                                           aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                       >  2:223710/1‑151 (MQ=255)
                                               ccGCGGAAAAAATATATAAAGAATTAAAAAAAAAAAAAGGCTC‑CCTTCGCCTTCCCCTAAAAAAATACCCCTCTCCCCAGGAGGGGGGCCGCTACCCTTA‑TTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACaa                                                                                                   <  1:769429/151‑1 (MQ=255)
                                                cGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                  >  2:194934/1‑151 (MQ=255)
                                                   tGGAAAAAAACCTAAAGAATTAAAAAAAAAAAAAGTTTC‑CCTCCGGCTCCCCCTGACAAAATCCCCTTCTCCCCACGAAGAGGGCCTTAACCCTTA‑CCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                 <  1:208391/149‑1 (MQ=255)
                                                    ggAAAATTATCTAAAGAATAACAAAAAAAAAAAGGCTC‑CCTTCGCCTTGCCCTGACAAAATACCCCCCTTCCCACGAAGGGGCCCGTTAACTTTA‑GCACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                               <  1:364684/150‑1 (MQ=255)
                                                    gaaagaaTATCAAAAGAATTAC‑AAAAAAAAAAGGCCACCCTTCGGCTGCCCAAGAAAAAACCCCCTTCTCCCCAGGAAGGGGCCCGCAACCCTGA‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                               <  1:256923/148‑1 (MQ=255)
                                                      aagaaTATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                            >  1:272133/1‑151 (MQ=255)
                                                      aaaaaTATCTACAGGATTACAAAAAAAAAAAGGCTC‑CCTTCGGCTTGCCCTGAAAAAATAGCCCTCTCCCCACGAAGGGGGCCCCTAACCTTA‑TTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                            <  2:1176530/148‑1 (MQ=255)
                                                      aaaaaTATCAACAGGATTAAAAAAAAAAAAAGGCTC‑CCTTCGCCTCCCCCTGAAAAAATACCCCTCTCCCCAGGAAGGGGGCCGCTACCCTAA‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                              <  1:995318/146‑1 (MQ=255)
                                                       agaaTATCTACAGGATTAAAAAAAAAAAAAGCCTC‑CCTTCGGCTTGCCCTGAAAAAATAGCCCTCTCCCCACGAAGGGGGCCGCTAACCTTA‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                              <  1:736980/148‑1 (MQ=255)
                                                       aaaaTATCTACGGGATTAAAAAAAAAAAAAGGCTC‑CTTTCGTCTTCCCCTAAAAAAATCCCCCTCCCCCCGGGAAGAGGGCCGCTAACTTTA‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                           <  1:1007574/149‑1 (MQ=255)
                                                         aaTATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                >  1:1373455/1‑144 (MQ=255)
                                                         aaTATCTACAGAATTAAAAAAAAAAAAAGCCTC‑CCTCCGGCTTCCCCGAAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCAACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                <  2:1373536/144‑1 (MQ=255)
                                                          atatCAAAAGGATGAGAAAAAAAAAAAGTCTC‑CTTTCGGCTCGCCGAAAAAAAATCCCCCTCCCCCCAAGAAGGGGCCCGCAACCTTTG‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                         <  1:381477/150‑1 (MQ=255)
                                                            atCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                                      >  1:35269/1‑151 (MQ=255)
                                                            accaaCAGGATTACAAAAAAAAAAAGCCTC‑CCTTCGTTTTCCCCAAAAAAAATACCCCTCTCCCCACGAGGAGGCCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                                      <  1:670619/147‑1 (MQ=255)
                                                                aaaggataaaaaaaaaaaaAGGGCTC‑CTTCCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCA‑GGAGGGGGCCCCCAACCTTAGGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                    <  1:826303/138‑1 (MQ=255)
                                                                  agatttACAAAAAAAAAAAGCCTC‑CCTTCGTCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAGGGGGCCCGCTACCCTTA‑CTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                                      <  1:1013410/141‑1 (MQ=255)
                                                                   gtattaaaaaaaaaaaaaGGCTC‑CCTTCGTCTCCCCCTGACAAAAAAGCCCTCTCCCCAGGAAGGGGCCCGCTAACCTTA‑TTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                       <  1:716855/135‑1 (MQ=255)
                                                                   ggATTACGAAAAAAAAAAGGCTC‑CCTTCGGCTTGCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                                               <  2:1363215/151‑1 (MQ=255)
                                                                   ggATAACAAAAAAAAAAAGCCTT‑CCTTCGGCTTGCCCTGAAAAATTGCCCCTCTCCCCACGAAGGGGCCCGTTAACTTTA‑GTTCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                <  2:954750/150‑1 (MQ=255)
                                                                     ataaaaaaaaaaaaaaGGCTC‑CCTTGGTCTTGCCCTAAAAAAATCCCCCTCTCCCCCCGAGGGGGGCCGCTAACCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                              <  1:589164/144‑1 (MQ=255)
                                                                          gAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                        >  1:312336/1‑151 (MQ=255)
                                                                           aaaaaaaaaaGGCTC‑CCTTCGCCTTCCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                       <  2:1236168/151‑1 (MQ=255)
                                                                           aaaaaaaaaaGGCTA‑CCTGCGCCTTGCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                              <  2:1368867/144‑1 (MQ=255)
                                                                           aaaTAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                       >  2:1287071/1‑151 (MQ=255)
                                                                            aaTAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                      >  2:951167/1‑151 (MQ=255)
                                                                             aaaaaaaaGGCTC‑CCTTCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCAACCCTTA‑TTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                     <  2:879953/149‑1 (MQ=255)
                                                                             aTAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACACACCGTCCGGCCCAACCATAATGCTTtg                                                                     >  1:482395/1‑151 (MQ=255)
                                                                                  aaaGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTa                                                                 >  1:1228883/1‑150 (MQ=255)
                                                                                  aaaGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                >  1:27300/1‑151 (MQ=255)
                                                                                   agggCTC‑CCTTCGGTTTCCCCTAAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTATCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACt                                                               <  1:795084/149‑1 (MQ=255)
                                                                                    aGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                              >  2:776510/1‑151 (MQ=255)
                                                                                     ggCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                             >  2:1140916/1‑151 (MQ=255)
                                                                                     ggCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                             >  1:149175/1‑151 (MQ=255)
                                                                                       cTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                           >  1:1187582/1‑151 (MQ=255)
                                                                                       cTA‑CCTTAGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCTGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                           >  1:1228946/1‑151 (MQ=255)
                                                                                        tA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                           >  2:963824/1‑150 (MQ=255)
                                                                                        tA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                           >  1:1426350/1‑150 (MQ=255)
                                                                                           ccTTCGCCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAGGGGGGCCGCTACCTTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                             <  1:1224145/146‑1 (MQ=255)
                                                                                            cTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                       >  1:777221/1‑151 (MQ=255)
                                                                                             ttCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                      >  1:1098641/1‑151 (MQ=255)
                                                                                             ttCGGCTTCCCCTGACAAAATCCCCCTCTCCCCACGAGGAGGGCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                        <  2:628758/149‑1 (MQ=255)
                                                                                                ggCTTCCCCTGACAAAATACCCCTCTCCCCACGAAGAGGCCCGCTAACTTAA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                       <  2:331310/147‑1 (MQ=255)
                                                                                                gccTTGCCGTGAAAAAATAGCCCTCTCCCCAAGAAGAGGGCCGCTACCTTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                      <  1:1528324/146‑1 (MQ=255)
                                                                                                  cTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                 >  1:734291/1‑151 (MQ=255)
                                                                                                    tGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                               >  1:1256692/1‑151 (MQ=255)
                                                                                                      cccTGACAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                              <  2:1452174/150‑1 (MQ=255)
                                                                                                         tGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                          >  2:783495/1‑151 (MQ=255)
                                                                                                         tGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                          >  2:1139900/1‑151 (MQ=255)
                                                                                                             aaaatccccccccccccACGAAGAGGGCCGCAACCCTAA‑TCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCAc                                       <  1:191467/136‑1 (MQ=255)
                                                                                                               aaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                    >  2:976484/1‑151 (MQ=255)
                                                                                                                    cccTCTCCCCACAAGGGGGCCCGCAAACCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGa                                <  2:1279900/150‑1 (MQ=255)
                                                                                                                        ctccccACGAAGAGCCCCGCTACCCTTA‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtt                            <  1:400698/147‑1 (MQ=255)
                                                                                                                          ccccAGGAAGGGGCCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtttt                          <  1:1496192/149‑1 (MQ=255)
                                                                                                                           ccccCGAAGAGGCCCGCTACCTTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtt                            <  1:308882/147‑1 (MQ=255)
                                                                                                                            ccACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                       >  2:306500/1‑151 (MQ=255)
                                                                                                                                gAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                    >  2:1048276/1‑150 (MQ=255)
                                                                                                                                gAAGAGGCCCGCTACCCTTA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                    <  2:131056/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcag                  >  1:1240181/1‑150 (MQ=255)
                                                                                                                                   gggtgCCGCTCCCTTTA‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagca                <  1:467084/147‑1 (MQ=255)
                                                                                                                                   gagGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                 >  1:1064893/1‑150 (MQ=255)
                                                                                                                                     ggcccGCTACCTTAA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                 <  2:63275/145‑1 (MQ=255)
                                                                                                                                     ggcccGCTACCCTAA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag               <  2:962588/147‑1 (MQ=255)
                                                                                                                                        ccGCTCCCCTAA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt            <  2:1068995/150‑1 (MQ=255)
                                                                                                                                        ccGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt            <  2:1377978/150‑1 (MQ=255)
                                                                                                                                            tAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg       >  1:845244/1‑151 (MQ=255)
                                                                                                                                              aCCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt      >  1:1505462/1‑150 (MQ=255)
                                                                                                                                               ccTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTAt    >  1:807185/1‑151 (MQ=255)
                                                                                                                                               ccTTA‑GACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTa     <  1:638898/150‑1 (MQ=255)
                                                                                                                                               ccTAA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg       <  1:1109121/148‑1 (MQ=255)
                                                                                                                                                   a‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc  <  2:1130003/149‑1 (MQ=255)
                                                                                                                                                    attACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc   <  2:315170/146‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTA‑CCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTA‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCC  >  NC_000913/3584937‑3585231

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