Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,358 T→A 8.3% L17Q (CTG→CAG)  ybgE → putative inner membrane protein in cydABX‑ybgE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,3580TA 8.3% 252.0 / 10.3 92L17Q (CTG→CAG) ybgEputative inner membrane protein in cydABX‑ybgE operon
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base A (3/5);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

TTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTG  >  NC_000913/774215‑774506
                                                                                                                                                 |                                                                                                                                                    
ttCTGGGCACGCTTCTTGCCTGGTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGgc                                                                                                                                                <  2:309643/150‑1 (MQ=255)
ttCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGgc                                                                                                                                                <  2:322849/150‑1 (MQ=255)
  cTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGgcgc                                                                                                                                              <  1:794985/150‑1 (MQ=255)
   tGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCt                                                                                                                                             >  2:1304041/1‑150 (MQ=255)
     ggAACGCTTCTTGCCTGTTCGTGTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAAGCAGGCAAAGCCTGTCAAGAAGACAGCTGAGGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCttt                                                                                                                                           <  1:115743/150‑1 (MQ=255)
      gAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCGCAGGGGGCGCttgt                                                                                                                                          <  1:400229/150‑3 (MQ=255)
      gAACGCTTCTTCCCTGGTCGTTTGGGGTAATCACCGCGCTGGCGCGTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTACGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCGCAGGGGGCGCttgt                                                                                                                                          <  1:203161/150‑3 (MQ=255)
       accGCTTCTTGCCGGGTCGTTGGGGCAAAACACCGCGCTGGCGCTTGAACACGTCGAAACACGCAAAGCCGGTCAAGACGACATCTGATGAGTAAGATTATCGCGACTATGTATGCGGTAATGGACTAGC‑G‑CCCCCTGCGGGCGCTTTcc                                                                                                                                         <  2:98412/148‑1 (MQ=255)
       aaCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACCCGTCGAAGCAGGCAAAGCCGGGCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCGCAGGGGGCGCttgtc                                                                                                                                         <  1:1444674/150‑4 (MQ=255)
        aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCGCAGGGGGCGCttgtcc                                                                                                                                        <  2:350313/150‑5 (MQ=255)
        aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCt                                                                                                                                        <  1:1025423/150‑1 (MQ=255)
        aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGGCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCt                                                                                                                                        <  1:1103118/150‑1 (MQ=255)
          gCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCGCAGGGGGCGCttgtccat                                                                                                                                      <  1:566785/150‑7 (MQ=255)
          gCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTc                                                                                                                                      <  1:322605/150‑1 (MQ=255)
            ttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGt                                                                                                                                    >  2:1367565/1‑150 (MQ=255)
            ttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGt                                                                                                                                    >  2:687180/1‑150 (MQ=255)
             tcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTg                                                                                                                                   <  2:1328834/150‑1 (MQ=255)
             tcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTAGGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTg                                                                                                                                   >  1:69443/1‑150 (MQ=255)
              cttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGa                                                                                                                                  >  2:1517386/1‑150 (MQ=255)
                 gCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATgg                                                                                                                               >  1:1017528/1‑150 (MQ=255)
                 gCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATgg                                                                                                                               <  1:1099506/150‑1 (MQ=255)
                   cTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGgcg                                                                                                                             <  2:1548682/150‑1 (MQ=255)
                    tGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGgcgc                                                                                                                            <  1:712957/150‑1 (MQ=255)
                        cGTTTGGGGTAATCTCCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAATCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAATC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCt                                                                                                                        <  1:1197919/150‑1 (MQ=255)
                         gTTTGGGGTAATCACCGCGCTGGCGCTTGAACTCGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTCTCGCGACCTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTg                                                                                                                       <  1:197615/150‑1 (MQ=255)
                         gTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTg                                                                                                                       >  1:1185297/1‑150 (MQ=255)
                           ttGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGtt                                                                                                                     <  1:606409/150‑1 (MQ=255)
                            tGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTa                                                                                                                    >  2:392468/1‑150 (MQ=255)
                             ggggTAATCACCGCGCTGGCGCTGGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAg                                                                                                                   <  2:1323778/150‑1 (MQ=255)
                              gggTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGc                                                                                                                  >  1:741994/1‑150 (MQ=255)
                                gTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAg                                                                                                                <  2:591469/150‑1 (MQ=255)
                                  aaTCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGa                                                                                                              <  2:1017528/150‑1 (MQ=255)
                                    tCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGatg                                                                                                            >  1:1229099/1‑150 (MQ=255)
                                              ggCGCTTGAACACGTCGAATCATGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTgg                                                                                                  <  2:297786/150‑1 (MQ=255)
                                              ggCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTgg                                                                                                  >  1:1031617/1‑150 (MQ=255)
                                               gcgcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTggg                                                                                                 <  1:607694/150‑1 (MQ=255)
                                                 ggttGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGAAATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAc                                                                                               <  2:1519805/148‑1 (MQ=255)
                                                 gcTTGAACACTTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAc                                                                                               <  2:1519808/150‑1 (MQ=255)
                                                 gcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGTAGACATCTGATGAGTATGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAc                                                                                               <  2:582464/150‑1 (MQ=255)
                                                 gcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAc                                                                                               <  1:1218845/150‑1 (MQ=255)
                                                    tGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCTC‑CCCCCTGCGGGCGCTTTACTTCGTGATGGCGCTTCTGATAGCAGGATGTATGTTTTGGGAccc                                                                                             <  1:1541476/150‑1 (MQ=255)
                                                    tGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑GCCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAccc                                                                                             <  1:108881/150‑1 (MQ=255)
                                                    tGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCa                                                                                            >  2:613927/1‑150 (MQ=255)
                                                    tGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCa                                                                                            <  2:1538787/150‑1 (MQ=255)
                                                     gAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGGATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAt                                                                                           >  2:1607660/1‑150 (MQ=255)
                                                        cacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtctc                                                                                        <  1:17821/150‑1 (MQ=255)
                                                        cacGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGA‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAtctc                                                                                        >  1:1008261/1‑150 (MQ=255)
                                                            tCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGttt                                                                                    <  1:1125567/150‑1 (MQ=255)
                                                              gAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTg                                                                                  >  1:359604/1‑150 (MQ=255)
                                                              gAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTg                                                                                  >  1:1103637/1‑150 (MQ=255)
                                                              gAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTg                                                                                  >  2:1496178/1‑150 (MQ=255)
                                                              caaTCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTg                                                                                  <  1:1449253/149‑1 (MQ=255)
                                                                  cAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGCCTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcgc                                                                              >  1:1286422/1‑150 (MQ=255)
                                                                  cAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcgc                                                                              >  2:999464/1‑150 (MQ=255)
                                                                  cAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGAAATGGACAAGC‑G‑CCCCCTGCGGGCGCGTTCCTTCGTGAACGCGCATCTGTTAGAAGGATGTATGTTTTTTGACCCATCTCTTTTTGCcgc                                                                              >  2:785860/1‑150 (MQ=255)
                                                                   aGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcgcg                                                                             >  1:1361432/1‑150 (MQ=255)
                                                                    ggCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGa                                                                            >  2:1633651/1‑150 (MQ=255)
                                                                              gTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGa                                                                  >  1:1499922/1‑150 (MQ=255)
                                                                                cAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCATTCGTGATGGCGCTTCTGATAGCAGGATGAATGTTTTGCGACCCAACTCGTATTGCCGCGAAGACCAGTGAAc                                                                >  1:93707/1‑150 (MQ=255)
                                                                                cAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCACCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTTTGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAAc                                                                <  2:337896/150‑1 (MQ=255)
                                                                                  agaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAAGTg                                                              >  1:202414/1‑150 (MQ=255)
                                                                                  agaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTg                                                              >  1:904031/1‑150 (MQ=255)
                                                                                   gaagaCATCAGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCAGCGGGCGCATTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTATGACTCGAAGACCACAGAACtca                                                             >  1:162439/1‑148 (MQ=255)
                                                                                   caagaCTTCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTgg                                                             <  1:1616135/149‑1 (MQ=255)
                                                                                    aagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGa                                                            <  1:666591/150‑1 (MQ=255)
                                                                                     agaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGaa                                                           <  1:379456/150‑1 (MQ=255)
                                                                                        cATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATc                                                        >  2:907434/1‑150 (MQ=255)
                                                                                        cATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑C‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATc                                                        >  1:1712468/1‑150 (MQ=255)
                                                                                         aTCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCt                                                       >  1:1237468/1‑150 (MQ=255)
                                                                                          tCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTg                                                      <  1:1517386/150‑1 (MQ=255)
                                                                                           cTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGATATCTgg                                                     >  1:1160931/1‑150 (MQ=255)
                                                                                           cTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTgg                                                     <  2:741994/150‑1 (MQ=255)
                                                                                            tgatgaGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGc                                                    >  1:777295/1‑150 (MQ=255)
                                                                                                 aGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGt                                               >  2:1578658/1‑150 (MQ=255)
                                                                                                  gTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                              <  1:1346762/150‑1 (MQ=255)
                                                                                                  gTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCGTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGtt                                              <  2:321119/150‑1 (MQ=255)
                                                                                                              cgACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCTTGGGTTATTGGTGAtgt                                  >  1:314128/1‑150 (MQ=255)
                                                                                                                aCTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTgg                                >  2:386391/1‑150 (MQ=255)
                                                                                                                 cTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTggg                               >  1:159044/1‑150 (MQ=255)
                                                                                                                  tttGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGc                              <  1:427540/150‑1 (MQ=255)
                                                                                                                  tttGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTGCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGc                              <  1:1064858/150‑1 (MQ=255)
                                                                                                                           ggTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGc                     >  2:452657/1‑150 (MQ=255)
                                                                                                                               gaaggaaaGC‑G‑CCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCAGATGTGGGCCGTCTGTGCTGgt                 >  2:257039/7‑150 (MQ=255)
                                                                                                                                 ggACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGGTGgtgt               >  2:1307667/1‑150 (MQ=255)
                                                                                                                                 ggACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGgtgt               >  2:679894/1‑150 (MQ=255)
                                                                                                                                   gaaaGC‑G‑CCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGa             >  2:410392/3‑150 (MQ=255)
                                                                                                                                        c‑g‑cCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCGCGGTTTGCCGCGAAGACCAGTGCACTGGAAATCGGGCATGGTTTATTGCTGATGTGGTCCGTCTGTGCTGGTGTGATTCAc        >  1:258098/1‑150 (MQ=255)
                                                                                                                                            cccccTGCGGGCGCTTTCCTTCGGGCTGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACgg      <  1:1097518/150‑1 (MQ=255)
                                                                                                                                            cccccGGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACgg      <  2:69443/150‑1 (MQ=255)
                                                                                                                                              cccTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCCCGGCg    >  1:734169/1‑150 (MQ=255)
                                                                                                                                                cTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGCCCCATCTCGTTTTGACGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTg  >  2:952659/1‑150 (MQ=255)
                                                                                                                                                cTGCGGGCGCTTTCCTTCGTGATGGCGCTGCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTg  <  1:1091323/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGC‑G‑CCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGTG  >  NC_000913/774215‑774506

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: ‑