Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 787,782 G→C 5.2% intergenic (‑187/+15) gpmA ← / ← galM phosphoglyceromutase 1/aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913787,7820GC 5.2% 274.7 / 5.1 92intergenic (‑187/+15)gpmA/galMphosphoglyceromutase 1/aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (36/50);  new base C (3/2);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCC  >  NC_000913/787634‑787930
                                                                                                                                                    |                                                                                                                                                    
ataaTTATATTGCCGCGACGAAGCATCAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCTCAGTGTTAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGa                                                                                                                                                     <  1:955398/147‑1 (MQ=255)
 tcATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGAc                                                                                                                                                    >  2:1214051/1‑150 (MQ=255)
    ttATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAg                                                                                                                                                 <  2:442682/150‑1 (MQ=255)
        atTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCg                                                                                                                                             >  2:173933/1‑150 (MQ=255)
           gCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTaa                                                                                                                                          >  1:977435/1‑150 (MQ=255)
            ccGCGACGCAGCAACAGCAATGCTAACGCATAACCAAAGCGAAAATCGTTGCTCAGTGTAAGGTTGTTTTGAATATTGAGTTGCAACTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAAc                                                                                                                                         <  2:224737/150‑1 (MQ=255)
            ccGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAAc                                                                                                                                         <  1:906586/150‑1 (MQ=255)
            ccGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAAc                                                                                                                                         >  1:1312432/1‑150 (MQ=255)
             cgcgACGAAGCAACAGGAATGCTTACGCATAACCATAGGGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATGGAGTTGCAAATATGTCGGAGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCGTCATGAGGGctttatc                                                                                                                                        <  1:713785/150‑7 (MQ=255)
             cgcgACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACa                                                                                                                                        <  1:495395/150‑1 (MQ=255)
                 acgaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCACTCATGACGAGGGCGTAACATTAc                                                                                                                                    >  1:402855/1‑150 (MQ=255)
                  cgaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACt                                                                                                                                   <  1:1312011/150‑1 (MQ=255)
                   taagcaaaagcaaTTCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGTTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTc                                                                                                                                  <  2:803193/149‑1 (MQ=255)
                   gacgCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATAATGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTc                                                                                                                                  <  1:1199439/147‑1 (MQ=255)
                   gaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTCAAGAATATGAGATAAACACATCATGACGAGGGAGTAACAGTAATc                                                                                                                                  >  1:837220/1‑150 (MQ=255)
                    aaGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCa                                                                                                                                 <  1:766935/150‑1 (MQ=255)
                     aGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAg                                                                                                                                >  1:578714/1‑150 (MQ=255)
                         aCAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAAt                                                                                                                            >  2:901588/1‑150 (MQ=255)
                          cAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATa                                                                                                                           >  2:189856/1‑150 (MQ=255)
                            gCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATaaa                                                                                                                         >  1:1112476/1‑150 (MQ=255)
                            gCAATGCTTACGCATAAACATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATaaa                                                                                                                         >  2:673652/1‑150 (MQ=255)
                             cAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAAc                                                                                                                        <  2:248164/150‑1 (MQ=255)
                                  cTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGTGTAACATTACTCAGCAATAAACTGata                                                                                                                   >  2:608282/1‑150 (MQ=255)
                                   ttACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGatat                                                                                                                  <  1:1250703/150‑1 (MQ=255)
                                      cGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTcc                                                                                                               >  2:334693/1‑150 (MQ=255)
                                       gCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCg                                                                                                              <  2:50973/150‑1 (MQ=255)
                                         aTAACCATAGCGTATATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGGGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTc                                                                                                            <  1:717258/150‑1 (MQ=255)
                                             ccATAGCGAAAATAGTGGCGCAGTGTAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCt                                                                                                       >  2:379178/1‑150 (MQ=255)
                                              cATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCt                                                                                                       <  2:1581/150‑1 (MQ=255)
                                                 aGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGa                                                                                                    >  2:866147/1‑150 (MQ=255)
                                                   cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAAt                                                                                                  <  1:260424/150‑1 (MQ=255)
                                                    gAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCTTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATa                                                                                                 >  1:296296/1‑150 (MQ=255)
                                                     aaaaTAGAGGCGCAGTGTAAGGTTGTTGTGAATAATGAGTAGCAAATATGTCGGTGTATGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTCATATTCCGTCAGGCTGGAATAc                                                                                                <  2:193848/150‑1 (MQ=255)
                                                           gTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCg                                                                                          <  1:1270519/150‑1 (MQ=255)
                                                              gcgcAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCa                                                                                       >  1:234787/1‑150 (MQ=255)
                                                               tgcAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTAGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACTTTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAg                                                                                      <  2:696952/149‑1 (MQ=255)
                                                               cgcAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAg                                                                                      <  2:466490/150‑1 (MQ=255)
                                                               cgcAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAg                                                                                      <  1:455852/150‑1 (MQ=255)
                                                                 cAGTGTAAGGTTGTTTTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGa                                                                                    <  1:611857/150‑1 (MQ=255)
                                                                    tgtAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGc                                                                                 <  1:1065020/150‑1 (MQ=255)
                                                                     gtAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCa                                                                                <  1:499142/150‑1 (MQ=255)
                                                                     gtAAGGTTGGTGTGAATATTGAGTTGCAAATATGTCGGGGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCa                                                                                <  2:134527/150‑1 (MQ=255)
                                                                       aaGGTTGTTGTGAATAATGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGCTAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAgg                                                                              <  1:643121/150‑1 (MQ=255)
                                                                           ttgttgTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGc                                                                          <  2:594495/150‑1 (MQ=255)
                                                                               tgtgAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGAAATAAACTGTTATTCCGACAGGGTGGAATACTATTCGCCTGGAGGCAGGATGCAGTc                                                                      >  2:287687/1‑150 (MQ=255)
                                                                                 tgAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCATGAAGCAGTCCg                                                                    >  2:82742/1‑150 (MQ=255)
                                                                                   aaTATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGt                                                                  >  1:1163272/1‑150 (MQ=255)
                                                                                     aatTGAGTTGCAAATAAGTCGGTTATTGCTGGTGATTTGAACAATAAGAGAAAAAGCCCTCATCACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTTGAATACTCTTCGCCAGGACGCAGGAAGCACTCCGGTTg                                                                <  1:881728/149‑1 (MQ=255)
                                                                                         gAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGc                                                            >  2:500177/1‑150 (MQ=255)
                                                                                         gAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCACCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGc                                                            <  2:883806/150‑1 (MQ=255)
                                                                                          aGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGcc                                                           >  2:150708/1‑150 (MQ=255)
                                                                                            ttGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCCTGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTACGCCAGGACGCAGGAAGCACTCCGGTAGCGGCCAt                                                         >  1:1000645/1‑150 (MQ=255)
                                                                                             tGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTATTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCAtt                                                        >  1:813458/1‑150 (MQ=255)
                                                                                                   tatGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggt                                                  <  1:99700/150‑1 (MQ=255)
                                                                                                    ttGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTATCATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtg                                                 <  2:517672/149‑1 (MQ=255)
                                                                                                     tGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCTGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtgg                                                <  1:681515/150‑1 (MQ=255)
                                                                                                     tGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtgg                                                >  2:1187263/1‑150 (MQ=255)
                                                                                                     tGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtgg                                                <  1:987142/150‑1 (MQ=255)
                                                                                                      ttCGGTGTTTGCTGGTGATTTGAACAATATGAGATAACGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATCTTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtggt                                               <  1:642281/149‑1 (MQ=255)
                                                                                                      gTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtggt                                               >  2:233270/1‑150 (MQ=255)
                                                                                                        cGGTGTTTGCTGGTGATTTGAACAATATGAGATAACGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTc                                             <  2:1002741/150‑1 (MQ=255)
                                                                                                        cGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTc                                             <  2:670056/150‑1 (MQ=255)
                                                                                                         ggTGTTTGCTGGTGATTTGAACATTATGAGATAAAGCCGTCATGACGAGGGCGTAACATTACTCAGCTATATACTGATATTCCGTGAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCATGGTGGTTCg                                            <  1:914043/150‑1 (MQ=255)
                                                                                                          gtgtTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCgg                                           <  1:111298/150‑1 (MQ=255)
                                                                                                          gtgtATGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTATTCGCCAGGACGCAGGAAGCAGTACGGTTGAGGCCATTCAGGGTGCTTCgg                                           >  2:272193/1‑150 (MQ=255)
                                                                                                           tgtATGCTGGTGATTTGAACAATATGAGATAAAGCCCGCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATCTTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCggg                                          <  2:1035061/150‑1 (MQ=255)
                                                                                                            gtTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGc                                         <  2:759083/150‑1 (MQ=255)
                                                                                                            gtTTGCAGGGTATTTGAACAATAAGAGATAAAGCCCTCATGACGAGGGCGAAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGGTGCGGCCATTCAGGGTGGTTCGGGc                                         <  1:110031/150‑1 (MQ=255)
                                                                                                            gtTAGCTGGTGATTTGAACAATATGAGATAAACCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGc                                         <  2:243609/150‑1 (MQ=255)
                                                                                                               tGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAAATGATATTCCGTCAGGCTGGAATACTCTTCGCCAGTACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGTTTCGGGCTGt                                      >  2:1283054/1‑150 (MQ=255)
                                                                                                                 cTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTcc                                    >  2:704085/1‑150 (MQ=255)
                                                                                                                  tGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCg                                   <  2:493013/150‑1 (MQ=255)
                                                                                                                   ggTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCgg                                  >  2:1220805/1‑150 (MQ=255)
                                                                                                                    gTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGt                                 <  1:173933/150‑1 (MQ=255)
                                                                                                                     tGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGAGTAACATTACTCAGCAATAAACTGATATTCCGTCAGTCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCTGGGTGGTTCGGGCTGTCCGGTa                                <  2:604407/150‑1 (MQ=255)
                                                                                                                      gATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAg                               >  2:1049013/1‑150 (MQ=255)
                                                                                                                           gAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACt                          >  2:132829/1‑150 (MQ=255)
                                                                                                                            aaCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTc                         >  2:1061776/1‑150 (MQ=255)
                                                                                                                               aaTATGAGATAAAGCCCTCATGACGAGGGTGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCt                      <  1:608282/150‑1 (MQ=255)
                                                                                                                                 aaatagaTAAACCCCTCATTACGCGTGCGTAACATTACTCAGCAATAAACTGATAATCAGTCAGGCTGGAATACTCTTCGCCTGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGGAGAAACTCGCttt                    <  2:250015/146‑1 (MQ=255)
                                                                                                                                     atgttacgCCCTCGTCATGAGGGCTTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCag                >  2:439340/8‑150 (MQ=255)
                                                                                                                                     atgttacgCCCTCGTCATGAGGGCTTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTCGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCag                >  1:93401/8‑150 (MQ=255)
                                                                                                                                     agaTAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGACCGGTAGAAACTCGCTTTCCag                <  2:1143531/150‑1 (MQ=255)
                                                                                                                                        ttacgCCCTCGTCATGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagc             >  2:454892/5‑150 (MQ=255)
                                                                                                                                        tAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagc             <  1:500177/150‑1 (MQ=255)
                                                                                                                                        tAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTGTCCAgcgc             >  1:231337/1‑150 (MQ=255)
                                                                                                                                         acaGCCCTCATGACGCGGGCGTAACATTACTCAGCAATAAACTGATCTTCCGTCAGGCTGGAATAATCTTCGCCAGGACGACGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGatcc            <  2:902731/148‑1 (MQ=255)
                                                                                                                                             cccTCATGACGCGGGCGTAACATTCCTCAGCAATAAACTGATATTCGGTCAGGCTGGAATCCTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTACCAGagccagcc        >  1:773321/1‑150 (MQ=255)
                                                                                                                                             cccTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagccagcc        <  1:31645/150‑1 (MQ=255)
                                                                                                                                               cTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCATGCTCGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGGAGAAACTCGCTTTCCAGAGCCAGCCCt      <  1:79580/150‑1 (MQ=255)
                                                                                                                                                  aTGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGc   <  1:778862/150‑1 (MQ=255)
                                                                                                                                                   tGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGcc  >  1:60213/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCC  >  NC_000913/787634‑787930

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