Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 810,116 A→G Y258C (TAC→TGC)  bioB → biotin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913810,1160AG100.0% 283.1 / NA 75Y258C (TAC→TGC) bioBbiotin synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/45);  total (30/45)

GGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGCA  >  NC_000913/809970‑810262
                                                                                                                                                  |                                                                                                                                                  
ggCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTg                                                                                                                                                <  1:278519/151‑1 (MQ=255)
   aaaGCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCTGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTgcgc                                                                                                                                             <  2:976437/151‑1 (MQ=255)
   aaaCCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTtgcgtgcg                                                                                                                                              <  2:609748/150‑1 (MQ=255)
    aaCCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGcc                                                                                                                                            >  1:689217/1‑151 (MQ=255)
         gCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCttt                                                                                                                                         <  1:404264/149‑1 (MQ=255)
         gCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCttt                                                                                                                                         <  1:460982/149‑1 (MQ=255)
              cgccgcCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCgg                                                                                                                                  >  1:126875/1‑151 (MQ=255)
               gccgccGGAAAGCGTTCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGACGATGCCTTTGATTTTATTCACACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGcc                                                                                                                                    >  2:827382/1‑148 (MQ=255)
               gccgccGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCgg                                                                                                                                  >  1:1139524/1‑150 (MQ=255)
               gccgccGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGa                                                                                                                                 <  1:11124/151‑1 (MQ=255)
                    cgGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGa                                                                                                                            <  2:573881/151‑1 (MQ=255)
                     ggAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAg                                                                                                                           <  1:292218/151‑1 (MQ=255)
                        aaGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGc                                                                                                                          >  2:1100303/1‑149 (MQ=255)
                         aGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGa                                                                                                                       <  1:861621/151‑1 (MQ=255)
                          gCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGAt                                                                                                                      <  2:490214/151‑1 (MQ=255)
                           cGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATg                                                                                                                     <  1:833771/151‑1 (MQ=255)
                                cAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacg                                                                                                                 >  2:471538/1‑150 (MQ=255)
                                 aaTCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGCTCGCGCGGATCATGATGCCAACCTCTTGCGGGCGCCTTTCTGCCGGACGCGAGCAGATgaacga                                                                                                                >  2:312797/1‑150 (MQ=255)
                                       cATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGAc                                                                                                          <  1:437887/150‑1 (MQ=255)
                                       cATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACt                                                                                                         <  1:1019077/151‑1 (MQ=255)
                                         tGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTc                                                                                                        >  1:207886/1‑150 (MQ=255)
                                             ggtgaaggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAgg                                                                                                     <  2:175212/149‑1 (MQ=255)
                                              gtgaaggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGa                                                                                                  <  2:1101724/151‑1 (MQ=255)
                                                 aaggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAtgt                                                                                               <  2:712474/151‑1 (MQ=255)
                                                  aggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAtgtg                                                                                              <  2:157587/151‑1 (MQ=255)
                                                  aggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAtgtg                                                                                              >  1:1071889/1‑151 (MQ=255)
                                                              cgccgcTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAgg                                                                                  <  2:572618/151‑1 (MQ=255)
                                                                ccgcTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGgcg                                                                                <  2:212097/151‑1 (MQ=255)
                                                                  gcTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGgcgc                                                                               >  1:623249/1‑150 (MQ=255)
                                                                  gcTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGgcgc                                                                               >  2:80883/1‑150 (MQ=255)
                                                                    ttGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCaa                                                                             >  2:982861/1‑150 (MQ=255)
                                                                    ttGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCaa                                                                             >  2:457353/1‑150 (MQ=255)
                                                                     tGCCGATAACGATGCTGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAAc                                                                           >  2:771340/1‑151 (MQ=255)
                                                                     tGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCaaa                                                                            <  2:644677/150‑1 (MQ=255)
                                                                       ccGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTc                                                                         >  2:799304/1‑151 (MQ=255)
                                                                       ccGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTc                                                                         <  2:253172/151‑1 (MQ=255)
                                                                          aTAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGAt                                                                      <  1:136352/151‑1 (MQ=255)
                                                                            aaCGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGAttt                                                                    <  2:559343/151‑1 (MQ=255)
                                                                            aaCGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGAtt                                                                     <  2:975424/150‑1 (MQ=255)
                                                                               gatgatGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCt                                                                 <  2:238085/151‑1 (MQ=255)
                                                                                 tgatgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTAc                                                               >  1:1202658/1‑151 (MQ=255)
                                                                                  gatgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTAc                                                               <  2:182975/150‑1 (MQ=255)
                                                                                  gatgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACg                                                              <  1:1090715/151‑1 (MQ=255)
                                                                                  gatgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACg                                                              <  2:1047309/151‑1 (MQ=255)
                                                                                   atgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACgg                                                             <  1:126643/151‑1 (MQ=255)
                                                                                       cGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTg                                                          <  2:644836/150‑1 (MQ=255)
                                                                                       cGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTg                                                          <  1:1103992/150‑1 (MQ=255)
                                                                                           gCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAAc                                                     >  1:768691/1‑151 (MQ=255)
                                                                                             cTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAActg                                                   >  1:819821/1‑151 (MQ=255)
                                                                                             cTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAAc                                                     <  2:548106/149‑1 (MQ=255)
                                                                                                  aTTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGAc                                              >  2:1202826/1‑151 (MQ=255)
                                                                                                     ttATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCa                                            <  2:1115290/150‑1 (MQ=255)
                                                                                                     ttATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCAc                                           >  1:1003630/1‑151 (MQ=255)
                                                                                                      tATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCAc                                           <  2:408739/150‑1 (MQ=255)
                                                                                                         tCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGcc                                        <  1:53390/150‑1 (MQ=255)
                                                                                                          cGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGc                                         >  1:766994/1‑148 (MQ=255)
                                                                                                           gCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGaa                                     <  1:1219314/151‑1 (MQ=255)
                                                                                                           gCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGa                                      >  1:908848/1‑150 (MQ=255)
                                                                                                            cACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGaa                                     <  2:612418/150‑1 (MQ=255)
                                                                                                            cACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAAt                                    >  1:1156745/1‑151 (MQ=255)
                                                                                                            cACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAAt                                    <  1:657574/151‑1 (MQ=255)
                                                                                                             aCCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATc                                   >  1:39428/1‑151 (MQ=255)
                                                                                                              ccATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATcc                                  <  1:771348/151‑1 (MQ=255)
                                                                                                              ccATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATcc                                  <  1:1024681/151‑1 (MQ=255)
                                                                                                                   gCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGgaa                              <  2:501831/150‑1 (MQ=255)
                                                                                                                          cgcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATa                       >  1:709143/1‑150 (MQ=255)
                                                                                                                           gcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATaaa                     <  2:248884/151‑1 (MQ=255)
                                                                                                                             ggATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGa                   <  2:992574/151‑1 (MQ=255)
                                                                                                                                tCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCt                >  1:910816/1‑151 (MQ=255)
                                                                                                                                 cATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCt                <  2:488345/150‑1 (MQ=255)
                                                                                                                                    gatgCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCa             >  1:475229/1‑150 (MQ=255)
                                                                                                                                        ccAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGt        <  1:286787/151‑1 (MQ=255)
                                                                                                                                          aaCCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGtt       >  2:433854/1‑150 (MQ=255)
                                                                                                                                          aaCCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTc      >  2:509058/1‑151 (MQ=255)
                                                                                                                                               ctTGCGTGCGCCTTTCTGCCGGACGCGGGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGc   <  2:185206/149‑1 (MQ=255)
                                                                                                                                               ctTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGCa  >  1:64537/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAACTGTTCCGCA  >  NC_000913/809970‑810262

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