Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,809,534 T→C intergenic (+41/‑83) ygaY → / → ygaZ pseudogene, major facilitator transporter superfamily;putative transport; Not classified; putative transport protein/putative L‑valine exporter, norvaline resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,809,5340TC97.2% 249.6 / ‑2.1 71intergenic (+41/‑83)ygaY/ygaZpseudogene, major facilitator transporter superfamily;putative transport; Not classified; putative transport protein/putative L‑valine exporter, norvaline resistance protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (34/35);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGCAACATGGCGGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGTGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGCAAAGA  >  NC_000913/2809387‑2809678
                                                                                                                                                   |                                                                                                                                                
ggCAACATGGCGGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGgcgt                                                                                                                                                <  2:986103/150‑1 (MQ=255)
ggCAACATGGCGGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGtt                                                                                                                                               <  1:430040/151‑1 (MQ=255)
    aCATGGCGGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTaaga                                                                                                                                           <  2:254780/151‑1 (MQ=255)
         gCGGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGccc                                                                                                                                      <  1:761411/151‑1 (MQ=255)
           ggTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCAt                                                                                                                                    >  1:974035/1‑151 (MQ=255)
            gTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATc                                                                                                                                   <  2:863763/151‑1 (MQ=255)
               gggCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAg                                                                                                                                 <  1:462414/150‑1 (MQ=255)
                   tGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTg                                                                                                                             >  1:1018323/1‑150 (MQ=255)
                    ggCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGAt                                                                                                                           >  2:331051/1‑151 (MQ=255)
                      cGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATaa                                                                                                                         <  2:951572/151‑1 (MQ=255)
                          tGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGtt                                                                                                                     >  2:711748/1‑151 (MQ=255)
                               ggctggcGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAg                                                                                                                  <  1:1162664/149‑1 (MQ=255)
                                gctggcGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGata                                                                                                               <  1:813672/151‑1 (MQ=255)
                                    gcgcGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATattca                                                                                                           <  1:1029610/151‑1 (MQ=255)
                                    gcgcGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATattca                                                                                                           <  1:967106/151‑1 (MQ=255)
                                       cgacgaTTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATattcattc                                                                                                        >  1:159601/1‑151 (MQ=255)
                                       cgacgaTTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATattcattc                                                                                                        >  2:1041876/1‑151 (MQ=255)
                                         acgaTTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATattcattc                                                                                                        >  1:810101/1‑149 (MQ=255)
                                         acgaTTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAg                                                                                                      <  2:1064628/151‑1 (MQ=255)
                                          cgaTTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGt                                                                                                     <  1:481309/151‑1 (MQ=255)
                                                cccTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTAttt                                                                                               <  2:393549/151‑1 (MQ=255)
                                                    ggTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTtata                                                                                            >  1:273091/1‑150 (MQ=255)
                                                       tAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTataata                                                                                         <  1:541328/150‑1 (MQ=255)
                                                       tAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAatat                                                                                        >  1:703621/1‑151 (MQ=255)
                                                       tAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAatat                                                                                        >  2:541721/1‑151 (MQ=255)
                                                        aaCTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAatat                                                                                        <  2:618318/150‑1 (MQ=255)
                                                        aaCTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATAtt                                                                                       <  2:247004/151‑1 (MQ=255)
                                                            tACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAAc                                                                                    >  1:680782/1‑150 (MQ=255)
                                                              ctggtctggtGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACaa                                                                                  <  1:807176/150‑1 (MQ=255)
                                                               tggtctggtGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATc                                                                                >  1:777299/1‑151 (MQ=255)
                                                                  tctggtGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGtt                                                                             >  1:769632/1‑151 (MQ=255)
                                                                    tggtggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGTGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTa                                                                            >  1:897611/1‑150 (MQ=255)
                                                                     ggtggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGtt                                                                             <  2:373474/148‑1 (MQ=255)
                                                                     ggtggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTaa                                                                           <  2:781483/150‑1 (MQ=255)
                                                                     ggtggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTa                                                                            <  1:646776/149‑1 (MQ=255)
                                                                     ggtggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTa                                                                            >  2:66397/1‑149 (MQ=255)
                                                                     ggtggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAg                                                                          <  2:856348/151‑1 (MQ=255)
                                                                        ggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGt                                                                       <  2:361987/151‑1 (MQ=255)
                                                                        ggCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGt                                                                       <  2:728089/151‑1 (MQ=255)
                                                                          cGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTac                                                                     >  2:719218/1‑151 (MQ=255)
                                                                           gAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTac                                                                     >  1:226640/1‑150 (MQ=255)
                                                                            aGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTacac                                                                   <  2:673035/151‑1 (MQ=255)
                                                                               ttttCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACActct                                                                >  1:778318/1‑151 (MQ=255)
                                                                                tttCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACActct                                                                <  1:730757/150‑1 (MQ=255)
                                                                                   cATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATgg                                                            <  2:225645/151‑1 (MQ=255)
                                                                                      cGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGaa                                                          >  2:852898/1‑150 (MQ=255)
                                                                                       gTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGaa                                                          >  1:348331/1‑149 (MQ=255)
                                                                                         cAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGcc                                                      >  1:175082/1‑151 (MQ=255)
                                                                                               gCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTc                                                >  2:1141917/1‑151 (MQ=255)
                                                                                                ccGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGGTATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTcc                                               >  2:1141532/1‑151 (MQ=255)
                                                                                                  gCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCcac                                             >  2:95840/1‑151 (MQ=255)
                                                                                                    aaaTTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAg                                           <  2:810281/151‑1 (MQ=255)
                                                                                                      aTTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGcc                                         >  2:950793/1‑151 (MQ=255)
                                                                                                              aTCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGTGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTg                                 >  1:490879/1‑151 (MQ=255)
                                                                                                              aTCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTg                                 >  1:1191155/1‑151 (MQ=255)
                                                                                                                  gCAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTc                             <  1:330971/151‑1 (MQ=255)
                                                                                                                   cAACCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCg                            >  2:691284/1‑151 (MQ=255)
                                                                                                                     aCCGTTTTCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGc                          <  2:613026/151‑1 (MQ=255)
                                                                                                                            tCGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCtt                    <  1:1041630/150‑1 (MQ=255)
                                                                                                                             cGTGGCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTc                   >  2:728636/1‑150 (MQ=255)
                                                                                                                                ggCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATgg               >  2:244781/1‑151 (MQ=255)
                                                                                                                                 gCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATgg               >  2:225156/1‑150 (MQ=255)
                                                                                                                                 gCTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATgg               >  1:718174/1‑150 (MQ=255)
                                                                                                                                  cTTTACACTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGa              >  1:506373/1‑150 (MQ=255)
                                                                                                                                      acacTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGa          <  2:226688/150‑1 (MQ=255)
                                                                                                                                      acacTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGAt         <  2:1050144/151‑1 (MQ=255)
                                                                                                                                        acTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGc       >  2:284637/1‑151 (MQ=255)
                                                                                                                                        acTTATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGc       <  1:98173/151‑1 (MQ=255)
                                                                                                                                          ttATAAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGCa      >  2:998498/1‑150 (MQ=255)
                                                                                                                                           tataAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGCaa     >  2:262270/1‑150 (MQ=255)
                                                                                                                                           tataAGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGCaa     <  1:456237/150‑1 (MQ=255)
                                                                                                                                              aaGGGCGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGCAAAGa  <  2:351408/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGCAACATGGCGGTTGGGCTGGCGTTTGTCTGGCTGGCGCGACGATTGCCCTGGTTAACTTACTGGTCTGGTGGCGAGGTTTTCATCGTCAGGAAGCCGCAAATTAAGCAATCAGCAACCGTTTTCGTGGCTTTACACTTATAAGGGTGTTAAGAAGCCCATCAGTCTGATAAGGTTAAGATATTCATTCAGTCTATTTATAATATTAACAATCGTTAAGCGTACACTCTATGGAAAGCCCTACTCCACAGCCTGCTCCTGGTTCGGCGACCTTCATGGAAGGATGCAAAGA  >  NC_000913/2809387‑2809678

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