Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,036,103 C→T 6.3% H112Y (CAT→TAT)  hyaD → putative hydrogenase 1 maturation protease HyaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,036,1030CT 6.3% 296.8 / 6.8 95H112Y (CAT→TAT) hyaDputative hydrogenase 1 maturation protease HyaD
Reads supporting (aligned to +/- strand):  ref base C (39/50);  new base T (3/3);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAA  >  NC_000913/1035954‑1036252
                                                                                                                                                     |                                                                                                                                                     
tGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGAc                                                                                                                                                       >  2:1338756/1‑150 (MQ=255)
tGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGAc                                                                                                                                                       <  2:1549830/150‑1 (MQ=255)
     gtcGATGCCATTGACTACGGGCTGGATCCTTGAACGCTGCGAAGCTATCGAGTAGAAGTTATTCCGGCTTATCTCAGCTCGATTAATATGAGCCTGCATCAGTAAAGTTTCTCCGAAGTGTTGGCGCTGGCGGTTTTCCGCTGACATCTg                                                                                                                                                  <  2:973238/149‑1 (MQ=255)
     ctcGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTg                                                                                                                                                  <  2:97452/150‑1 (MQ=255)
     ctcGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTg                                                                                                                                                  <  1:1235005/150‑1 (MQ=255)
       cGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGGACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGcc                                                                                                                                                <  2:1182708/150‑1 (MQ=255)
       cGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGcc                                                                                                                                                <  1:181772/150‑1 (MQ=255)
         ctGCCATTGTCTACGGGCTGGAACCTGGAACCCTGCGAACCTATGCCTGATAACGCACTGCCGCTTATCTCAGCGCGATGAAACTGAGCCTGCATCAGAACAGCTTCTCCGTAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAg                                                                                                                                              <  2:947968/149‑1 (MQ=255)
         aTGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTACCAAGTGTTGGCGCTGGCGGATATCAGAGGACATCTGCCAg                                                                                                                                              >  2:1134448/1‑150 (MQ=255)
             cATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcaca                                                                                                                                          <  2:289100/150‑1 (MQ=255)
             cATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGcaca                                                                                                                                          >  1:460771/1‑150 (MQ=255)
                tGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcgccaac                                                                                                                                       <  2:1312739/150‑7 (MQ=39)
                 gACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAtt                                                                                                                                      >  1:1578610/1‑150 (MQ=255)
                 aaCTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAtt                                                                                                                                      <  1:1069304/149‑1 (MQ=255)
                   cTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGc                                                                                                                                    <  1:1000503/150‑1 (MQ=255)
                   cTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGc                                                                                                                                    <  1:545028/150‑1 (MQ=255)
                   cTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGc                                                                                                                                    <  1:1552823/150‑1 (MQ=255)
                   cTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGc                                                                                                                                    <  1:406864/150‑1 (MQ=255)
                         gCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                                                                            >  2:1295098/1‑135 (MQ=255)
                         gCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCAGATGTCCGCGGATATCCGCCAGcg                                                                                                                                            <  1:1295098/136‑2 (MQ=255)
                             gaacctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGt                                                                                                                          >  2:509862/1‑150 (MQ=255)
                              aacctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTc                                                                                                                         >  1:189566/1‑150 (MQ=255)
                                 ctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                                                                      >  2:1363925/1‑150 (MQ=255)
                                 ctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                                                                      >  1:1386168/1‑150 (MQ=255)
                                 ctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGc                                                                                                                      <  2:952912/150‑1 (MQ=255)
                                 ctggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGGCGGTCTGc                                                                                                                      <  2:486120/150‑1 (MQ=255)
                                  tggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCa                                                                                                                     >  2:1049546/1‑150 (MQ=255)
                                   ggaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCaa                                                                                                                    >  2:612675/1‑150 (MQ=255)
                                    gaacGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAAc                                                                                                                   <  2:1470719/150‑1 (MQ=255)
                                        gCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGc                                                                                                               >  1:1420095/1‑150 (MQ=255)
                                              gagataagCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCt                                                                                                         >  1:1605640/8‑150 (MQ=255)
                                              aaCCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCATCGTCGGTCTGCAACCCGCAATGCt                                                                                                         >  1:1432487/1‑150 (MQ=255)
                                               aCCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCt                                                                                                                        >  2:755629/1‑134 (MQ=255)
                                               aCCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCt                                                                                                                        <  1:755629/134‑1 (MQ=255)
                                                ccTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcg                                                                                                       <  1:173733/150‑1 (MQ=255)
                                                ccTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcg                                                                                                       <  1:862950/150‑1 (MQ=255)
                                                  tATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTcgac                                                                                                     >  2:58895/1‑150 (MQ=255)
                                                      ccGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACt                                                                                                 <  1:1557069/150‑1 (MQ=255)
                                                       cGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTa                                                                                                >  2:360344/1‑150 (MQ=255)
                                                           gaACGCATTCCGGCTGATCTGAGCGCGAAGAAAATGAGCCTGCATCAGAAGAGTTTCTCCGAAGTGTTGGCGCTGTCGGATATTCTCTGACATCTGTCAGCACATATTTTCTTCGTAGGGCTTCAATCTGAAATTCTCGAAGTCTAtggc                                                                                            <  2:706569/150‑1 (MQ=255)
                                                             aCGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTAcggcgg                                                                                          <  1:961351/150‑1 (MQ=255)
                                                              cGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGt                                                                                         >  1:1267203/1‑150 (MQ=255)
                                                                cATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAg                                                                                       >  1:870430/1‑150 (MQ=255)
                                                                  ttCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGCCATCTGCCAGCACATATTGCCCTCGTCGGTCTGCACCCCGCAAAGGTCGACGACTACGGCGGTAGcc                                                                                     >  2:1050200/1‑150 (MQ=255)
                                                                  ttCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGcc                                                                                     <  1:1049546/150‑1 (MQ=255)
                                                                  ttCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGcc                                                                                     >  2:1575338/1‑150 (MQ=255)
                                                                     cGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGa                                                                                  >  1:189409/1‑150 (MQ=255)
                                                                      ggCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAg                                                                                 <  2:622437/150‑1 (MQ=255)
                                                                      ggCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAg                                                                                 <  2:1578610/150‑1 (MQ=255)
                                                                           aTCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAAc                                                                            <  1:125094/150‑1 (MQ=255)
                                                                            tctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACt                                                                           <  1:835449/150‑1 (MQ=255)
                                                                             ctcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTg                                                                          <  2:1134501/150‑1 (MQ=255)
                                                                              tcAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTgg                                                                         <  1:538307/150‑1 (MQ=255)
                                                                               cAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGc                                                                        >  1:1338536/1‑150 (MQ=255)
                                                                                   gcgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACgg                                                                    >  1:1088551/1‑150 (MQ=255)
                                                                                    cgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACggg                                                                   >  1:806763/1‑150 (MQ=255)
                                                                                    cgAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACggg                                                                   <  2:480461/150‑1 (MQ=255)
                                                                                     gaagaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGa                                                                  <  2:1432487/150‑1 (MQ=255)
                                                                                      aagaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAg                                                                 <  1:191476/150‑1 (MQ=255)
                                                                                        gaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCa                                                               <  2:1226437/150‑1 (MQ=255)
                                                                                        gaaAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCa                                                               <  1:248398/150‑1 (MQ=255)
                                                                                         aaaaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCaa                                                              >  2:1538254/1‑150 (MQ=255)
                                                                                         aaaaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCaa                                                              >  1:450113/1‑150 (MQ=255)
                                                                                          aaaTGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAAc                                                             >  2:1444226/1‑150 (MQ=255)
                                                                                               aGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGccc                                                        <  1:1575338/150‑1 (MQ=255)
                                                                                                    gCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGc                                                   <  1:963892/150‑1 (MQ=255)
                                                                                                    gCATCAGAACAGTTGCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGc                                                   <  1:1197932/150‑1 (MQ=255)
                                                                                                     cATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCg                                                  >  2:952115/1‑150 (MQ=255)
                                                                                                        cAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGaa                                               >  2:467713/1‑150 (MQ=255)
                                                                                                         aGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAAc                                              >  1:954928/1‑150 (MQ=255)
                                                                                                         aGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAAc                                              >  1:617195/1‑150 (MQ=255)
                                                                                                          gAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACa                                             >  2:1477777/1‑150 (MQ=255)
                                                                                                          gAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACa                                             >  2:773063/1‑150 (MQ=255)
                                                                                                           aaCAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAg                                            <  2:384277/150‑1 (MQ=255)
                                                                                                             cAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAggc                                          <  2:1420095/150‑1 (MQ=255)
                                                                                                             cAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAggc                                          <  1:1047137/150‑1 (MQ=255)
                                                                                                             cAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAggc                                          >  1:658080/1‑150 (MQ=255)
                                                                                                                    tcCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTg                                   >  1:972765/1‑150 (MQ=255)
                                                                                                                     ccGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTgg                                  >  2:328294/1‑150 (MQ=255)
                                                                                                                     ccGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATAATGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTgg                                  >  2:43430/1‑150 (MQ=255)
                                                                                                                              tGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtg                         >  1:1039927/1‑150 (MQ=255)
                                                                                                                              tGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtg                         <  2:513629/150‑1 (MQ=255)
                                                                                                                               ggCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgc                        <  1:994136/150‑1 (MQ=255)
                                                                                                                                gtgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgct                       >  2:28987/3‑150 (MQ=255)
                                                                                                                                 tgcTGGCAGATGTCCGCGGATATCCGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTtgctg                      >  1:1213680/2‑150 (MQ=255)
                                                                                                                                     ggCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgcc                  <  2:460771/150‑3 (MQ=255)
                                                                                                                                      gCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTgg                 <  1:360344/150‑1 (MQ=255)
                                                                                                                                      gCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTgg                 <  1:1366233/150‑1 (MQ=255)
                                                                                                                                       cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCTAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTggg                <  2:1386168/150‑1 (MQ=255)
                                                                                                                                       cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTggg                >  1:419781/1‑150 (MQ=255)
                                                                                                                                       cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTggg                >  2:652803/1‑150 (MQ=255)
                                                                                                                                       cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgccag                <  1:328294/150‑5 (MQ=255)
                                                                                                                                       cGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCAAGCTGCgccag                <  1:224060/150‑5 (MQ=255)
                                                                                                                                             tCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATtg          <  1:403206/150‑1 (MQ=255)
                                                                                                                                                gcgGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGc       <  1:58895/150‑1 (MQ=255)
                                                                                                                                                  ggACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCg     <  2:450113/150‑1 (MQ=255)
                                                                                                                                                  ggACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCg     <  2:1088551/150‑1 (MQ=255)
                                                                                                                                                     cATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCaa  >  1:689437/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCGAAGAAAATGAGCCTGCATCAGAACAGTTTCTCCGAAGTGTTGGCGCTGGCGGATATCCGCGGACATCTGCCAGCACATATTGCCCTCGTCGGTCTGCAACCCGCAATGCTCGACGACTACGGCGGTAGCCTGAGCGAACTGGCACGGGAGCAACTGCCCGCTGCGGAACAGGCGGCGCTGGCGCAGCTTGCTGCGTGGGGAATTGTGCCGCAA  >  NC_000913/1035954‑1036252

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