Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,892,362 G→A G50D (GGC→GAC)  queD → 6‑carboxy‑5,6,7,8‑tetrahydropterin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,892,3620GA100.0% 320.1 / NA 90G50D (GGC→GAC) queD6‑carboxy‑5,6,7,8‑tetrahydropterin synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (50/40);  total (50/40)

TATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATC  >  NC_000913/2892230‑2892496
                                                                                                                                    |                                                                                                                                      
tATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGAtt                                                                                                                                 >  1:308046/1‑140 (MQ=255)
tATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGAtt                                                                                                                                 >  1:1382076/1‑140 (MQ=255)
  tttAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTAt                                                                                                                               >  2:1200293/1‑140 (MQ=255)
  tttAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTAt                                                                                                                               >  1:341988/1‑140 (MQ=255)
  tttAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTAt                                                                                                                               >  2:667225/1‑140 (MQ=255)
  tttAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTAt                                                                                                                               <  1:1118509/140‑1 (MQ=255)
   ttAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATc                                                                                                                              >  1:1367514/1‑140 (MQ=255)
     aaaGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGa                                                                                                                            >  2:1271814/1‑140 (MQ=255)
          tttCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCg                                                                                                                       <  2:252609/140‑1 (MQ=255)
            tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCt                                                                                                                     >  2:237620/1‑140 (MQ=255)
            tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCt                                                                                                                     >  2:763773/1‑140 (MQ=255)
            tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCt                                                                                                                     >  2:1208300/1‑140 (MQ=255)
            tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCt                                                                                                                     >  2:1206456/1‑140 (MQ=255)
            tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCt                                                                                                                     >  1:225145/1‑140 (MQ=255)
            tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCt                                                                                                                     >  1:543297/1‑140 (MQ=255)
             cACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTg                                                                                                                    >  1:769032/1‑140 (MQ=255)
                   cGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTaa                                                                                                              >  2:813092/1‑140 (MQ=255)
                       gCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAgc                                                                                                          <  1:48292/140‑1 (MQ=255)
                        ccGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAgcg                                                                                                         <  1:248140/140‑1 (MQ=255)
                        ccGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAgcg                                                                                                         <  1:53938/140‑1 (MQ=255)
                                  cTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAgcg                                                                                                         >  2:1155525/1‑130 (MQ=255)
                                  cTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAgcg                                                                                                         <  1:1155525/130‑1 (MQ=255)
                                  cTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaac                                                                                               <  2:769032/140‑1 (MQ=255)
                                        acacGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaacc                                                                                              <  2:440594/135‑1 (MQ=255)
                                        acacGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaacc                                                                                              >  1:440594/1‑135 (MQ=255)
                                          acGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTAc                                                                                       >  2:1253952/1‑140 (MQ=255)
                                          acGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTAc                                                                                       >  1:676381/1‑140 (MQ=255)
                                          acGTCCCGGAAGGCCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTAc                                                                                       >  1:880869/1‑140 (MQ=255)
                                             tCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaac                                                                                               >  1:462019/1‑129 (MQ=255)
                                             tCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaac                                                                                               <  2:462019/129‑1 (MQ=255)
                                             tCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAg                                                                                    >  1:699409/1‑140 (MQ=255)
                                             tCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAg                                                                                    >  2:412656/1‑140 (MQ=255)
                                               ccGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaac                                                                                               >  2:511927/1‑127 (MQ=255)
                                               ccGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaac                                                                                               <  1:511927/127‑1 (MQ=255)
                                               ccGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAgcg                                                                                  >  1:1335887/1‑140 (MQ=255)
                                               ccGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAgcg                                                                                  >  1:818580/1‑140 (MQ=255)
                                                  gAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCt                                                                               <  2:225145/140‑1 (MQ=255)
                                                  gAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCt                                                                               <  1:131934/140‑1 (MQ=255)
                                                     gggCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGa                                                                            <  1:1206456/140‑1 (MQ=255)
                                                       gCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATc                                                                          <  2:1335887/140‑1 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt                                                                           >  2:308221/1‑138 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt                                                                           <  1:308221/138‑1 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCa                                                                         >  2:191487/1‑140 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCa                                                                         >  2:1204475/1‑140 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCa                                                                         >  1:339734/1‑140 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCa                                                                         >  2:287032/1‑140 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCa                                                                         >  1:935210/1‑140 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCa                                                                         >  1:602778/1‑140 (MQ=255)
                                                        cATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCa                                                                         >  1:459073/1‑140 (MQ=255)
                                                            aaTGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCAt                                                                     <  2:384286/140‑1 (MQ=255)
                                                                 ggTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTAtct                                                                <  1:763773/140‑1 (MQ=255)
                                                                 ggTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTAtct                                                                <  1:650780/140‑1 (MQ=255)
                                                                 ggTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTAtct                                                                <  2:944306/140‑1 (MQ=255)
                                                                  gTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGa                                                                            <  1:935697/127‑1 (MQ=255)
                                                                  gTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGa                                                                            >  2:935697/1‑127 (MQ=255)
                                                                     gCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAAt                                                            <  1:909116/140‑1 (MQ=255)
                                                                     gCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAAt                                                            <  1:412656/140‑1 (MQ=255)
                                                                      ccTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt                                                                           >  2:1334491/1‑124 (MQ=255)
                                                                      ccTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt                                                                           <  1:1334491/124‑1 (MQ=255)
                                                                      ccTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATg                                                           >  2:438345/1‑140 (MQ=255)
                                                                            cGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTc                                                     <  1:1200293/140‑1 (MQ=255)
                                                                            cGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTc                                                     <  1:667225/140‑1 (MQ=255)
                                                                               gCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAg                                                  >  1:1169814/1‑140 (MQ=255)
                                                                                cATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAgg                                                 <  2:676381/140‑1 (MQ=255)
                                                                                         aTGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGaaaa                                        >  2:334416/1‑140 (MQ=255)
                                                                                             tGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCa                                    >  2:991793/1‑140 (MQ=255)
                                                                                             tGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCa                                    >  2:982305/1‑140 (MQ=255)
                                                                                             tGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCa                                    >  1:1387683/1‑140 (MQ=255)
                                                                                              gCGTCTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCaa                                   <  1:528374/140‑1 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCa                                    >  1:1304380/1‑138 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCa                                    <  2:1304380/138‑1 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAAc                                  >  1:949767/1‑140 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAAc                                  >  1:390769/1‑140 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAAc                                  <  1:549666/140‑1 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAAc                                  >  1:229126/1‑140 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAAc                                  <  2:1169814/140‑1 (MQ=255)
                                                                                               cGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGAAATTCCAGGTCTGGAAAACCCAAc                                  >  2:659933/1‑140 (MQ=255)
                                                                                                  cTGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAg                               <  2:1382076/140‑1 (MQ=255)
                                                                                                   tGGAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGAAATTCCAGGTCTGGAAAACCCAACCAGc                              <  1:659933/140‑1 (MQ=255)
                                                                                                    ggAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCg                             <  1:438345/140‑1 (MQ=255)
                                                                                                     gAAATTACCGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGa                            <  2:459073/140‑1 (MQ=255)
                                                                                                             cGGGGAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAg                    <  2:341988/140‑1 (MQ=255)
                                                                                                               gggAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCa                  >  1:1302629/1‑140 (MQ=255)
                                                                                                               gggAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCa                  >  2:957123/1‑140 (MQ=255)
                                                                                                               gggAAGTCGATCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCa                  >  1:676642/1‑140 (MQ=255)
                                                                                                                        aTCCGCATACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGAtt         <  1:237620/140‑1 (MQ=255)
                                                                                                                              aTACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGAt   <  2:880869/140‑1 (MQ=255)
                                                                                                                              aTACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGAt   <  1:1398135/140‑1 (MQ=255)
                                                                                                                              aTACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGAt   <  1:1271814/140‑1 (MQ=255)
                                                                                                                              aTACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATc  >  2:666832/1‑140 (MQ=255)
                                                                                                                               tACGGACTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATc  <  2:949767/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATC  >  NC_000913/2892230‑2892496

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