Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 442,390 C→T G159D (GGT→GAT)  thiJ ← 4‑methyl‑5(beta‑hydroxyethyl)‑thiazole monophosphate synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913442,3900CT99.7% 1057.3 / ‑5.8 293G159D (GGT→GAT) thiJ4‑methyl‑5(beta‑hydroxyethyl)‑thiazole monophosphate synthesis
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (192/100);  total (193/100)
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.88e-01

Only 100 of 305 total aligned reads displayed.

GCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGC  >  NC_000913/442300‑442479
                                                                                          |                                                                                         
   gccATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGc                                                                                      >  1:388495/1‑93 (MQ=255)
    ccATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcc                                                                                     >  1:7752313/1‑93 (MQ=255)
    ccATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcc                                                                                     >  2:4412871/1‑93 (MQ=255)
    ccATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcc                                                                                     >  2:2003923/1‑93 (MQ=255)
    ccATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcc                                                                                     >  2:1881068/1‑93 (MQ=255)
     cATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTGTCCGGccc                                                                                    >  1:7117262/1‑93 (MQ=255)
     cATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGccc                                                                                    <  2:243544/93‑1 (MQ=255)
     cATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGccc                                                                                    >  2:4173219/1‑93 (MQ=255)
      aTCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcccc                                                                                   >  1:6251970/1‑93 (MQ=255)
      aTCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcccc                                                                                   <  2:6852156/93‑1 (MQ=255)
      aTCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcccc                                                                                   <  2:5781753/93‑1 (MQ=255)
      aTCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGcccc                                                                                   >  1:4547959/1‑93 (MQ=255)
       tCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCct                                                                                  >  1:7239102/1‑93 (MQ=255)
        cACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCctg                                                                                 <  1:6166843/93‑1 (MQ=255)
        cACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCctg                                                                                 <  1:6570123/93‑1 (MQ=255)
        cACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCctg                                                                                 >  2:3503804/1‑93 (MQ=255)
         aCCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCctgg                                                                                <  2:6449507/93‑1 (MQ=255)
          ccAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCctggc                                                                               >  1:285231/1‑93 (MQ=255)
           cAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCctggct                                                                              >  1:5135054/1‑93 (MQ=255)
            aGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCctggctg                                                                             <  2:7274845/93‑1 (MQ=255)
              ttGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATCGCTGTATCCGGCCCCTGGCTGGt                                                                           >  2:4282926/1‑93 (MQ=255)
              ttGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGt                                                                           >  1:139340/1‑93 (MQ=255)
               tgtgATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTc                                                                          <  2:4583465/93‑1 (MQ=255)
                gtgATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCa                                                                         >  2:2808467/1‑93 (MQ=255)
                gtgATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCa                                                                         <  1:261059/93‑1 (MQ=255)
                 tgATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAg                                                                        >  1:245051/1‑93 (MQ=255)
                  gATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGc                                                                       <  1:1881251/93‑1 (MQ=255)
                   aTGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCa                                                                      >  1:2940449/1‑93 (MQ=255)
                    tGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCaa                                                                     >  2:1592734/1‑93 (MQ=255)
                     gCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAAt                                                                    <  1:7552780/93‑1 (MQ=255)
                     gCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAAt                                                                    >  2:6493047/1‑93 (MQ=255)
                        aCTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAAtttt                                                                 >  1:4198138/1‑93 (MQ=255)
                         cTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTa                                                                >  2:7140807/1‑93 (MQ=255)
                            cATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTAccc                                                             >  1:7414110/1‑93 (MQ=255)
                               ggtcTTTTTCACGCCCAACCAACAGGTCGATAATTTCCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                                                          <  1:6943601/90‑1 (MQ=255)
                                ggCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCATCAATTTTACCCGTGc                                                         <  2:5388987/93‑1 (MQ=255)
                                ggCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGc                                                         <  1:1974026/93‑1 (MQ=255)
                                ggCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGc                                                         <  1:5293529/93‑1 (MQ=255)
                                 gCTTTTTCTCGCCCAAACAACAGTTCGATAATTTTCAGACCAAAGTCTATAGCTGTATCCGGCCCCCGGCTGGTCAGCAATTTTAACCGTGCa                                                        <  2:1692518/93‑1 (MQ=255)
                                 gCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCa                                                        >  1:6386782/1‑93 (MQ=255)
                                 gCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCa                                                        >  1:39583/1‑93 (MQ=255)
                                 gCTTTTCCACGCCCAGCCAACGGGTCGATAATTGTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCa                                                        <  1:6799089/93‑1 (MQ=255)
                                    ttttCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATcc                                                     >  2:2667512/1‑93 (MQ=255)
                                     tttCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATccc                                                    <  2:6373928/93‑1 (MQ=255)
                                      ttCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCa                                                   <  2:4729163/93‑1 (MQ=255)
                                        cACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAga                                                 <  2:3272142/93‑1 (MQ=255)
                                        cACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAga                                                 <  1:3670147/93‑1 (MQ=255)
                                         aCGCCCAACCAACAGGTCGATAGTTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAgac                                                >  2:4978146/1‑93 (MQ=255)
                                         aCGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCacac                                                <  1:2037227/93‑1 (MQ=255)
                                           gCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAgacga                                              >  2:6595877/1‑93 (MQ=255)
                                            cccAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAgacgac                                             <  2:7625300/93‑1 (MQ=255)
                                              caaccaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGAcgc                                           <  2:7463486/93‑1 (MQ=255)
                                                accaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACgcgc                                         <  1:5467934/93‑1 (MQ=255)
                                                accaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACgcgc                                         <  1:1548065/93‑1 (MQ=255)
                                                 ccaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt                                        <  2:4690842/93‑1 (MQ=255)
                                                 ccaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt                                        >  2:7748276/1‑93 (MQ=255)
                                                 ccaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt                                        >  1:2315296/1‑93 (MQ=255)
                                                  caacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCtt                                       <  1:7165295/93‑1 (MQ=255)
                                                  caacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCtt                                       >  1:667782/1‑93 (MQ=255)
                                                  caacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCtt                                       >  1:5507372/1‑93 (MQ=255)
                                                   cacAGGTCGATAATTTTCAGACCAAAGTCGATACCTGTATCCGGCCCCTGGCTTGTCAGCAATTTTTCCCGTGCATCCCAGACGACGCGCTTg                                      <  2:6072136/92‑1 (MQ=255)
                                                   aacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTg                                      >  1:7686661/1‑93 (MQ=255)
                                                     cAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTc                                    >  1:437930/1‑93 (MQ=255)
                                                       ggTCGATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCa                                  <  1:5364817/93‑1 (MQ=255)
                                                           gATAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCa                              >  2:1239307/1‑93 (MQ=255)
                                                            aTAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCAt                             >  2:2982493/1‑93 (MQ=255)
                                                             tAATTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCCGCCAtt                            >  2:7302484/1‑93 (MQ=255)
                                                              aaTTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTg                           <  2:5632926/93‑1 (MQ=255)
                                                               aTTTTCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGt                          <  2:6765308/93‑1 (MQ=255)
                                                                ttttCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCTCAGACGACGCGCTTGTCCAGCCATTGtt                         >  1:4726718/1‑93 (MQ=255)
                                                                ttttCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGtt                         >  2:6027223/1‑93 (MQ=255)
                                                                ttttCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGtt                         <  1:3381162/93‑1 (MQ=255)
                                                                ttttCAGACCAAAGTCGATAGCTGTATCCGGCACCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGtt                         >  2:5416168/1‑93 (MQ=255)
                                                                 tttCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTc                        <  1:6445841/93‑1 (MQ=255)
                                                                 tttCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTc                        <  1:7387676/93‑1 (MQ=255)
                                                                   tCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTcgg                      >  1:5863498/1‑93 (MQ=255)
                                                                   tCAGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTcgg                      >  1:5932310/1‑93 (MQ=255)
                                                                     aGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTcggcg                    >  1:1246601/1‑93 (MQ=255)
                                                                     aGACCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTcggcg                    <  1:901966/93‑1 (MQ=255)
                                                                       aCCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCggg                  >  2:4103494/1‑93 (MQ=255)
                                                                       aCCAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCggg                  >  1:4333327/1‑93 (MQ=255)
                                                                        ccAAAGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGa                 >  2:2230435/1‑93 (MQ=255)
                                                                          aaaGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAAt               >  2:2879146/1‑93 (MQ=255)
                                                                            aGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAAttt             >  1:7928299/1‑93 (MQ=255)
                                                                            aGTCGATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAAttt             <  1:4798982/93‑1 (MQ=255)
                                                                                gATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTc         >  2:3659549/1‑93 (MQ=255)
                                                                                gATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTc         <  1:5993199/93‑1 (MQ=255)
                                                                                gATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTc         >  1:5746189/1‑93 (MQ=255)
                                                                                gATAGCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTc         >  1:1119926/1‑93 (MQ=255)
                                                                                    gCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTc     >  2:1472429/1‑93 (MQ=255)
                                                                                    gCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTc     >  2:4009325/1‑93 (MQ=255)
                                                                                    gCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTc     >  2:4365330/1‑93 (MQ=255)
                                                                                    gCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTc     >  1:581994/1‑93 (MQ=255)
                                                                                    gCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTATTTc     >  1:2864206/1‑93 (MQ=255)
                                                                                    gCTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGACGGGAATTTTGTCTTTc     >  2:4009329/1‑93 (MQ=255)
                                                                                     cTGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCa    <  2:5213893/93‑1 (MQ=255)
                                                                                      tGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAg   >  2:3744866/1‑93 (MQ=255)
                                                                                      tGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAg   <  1:1494781/93‑1 (MQ=255)
                                                                                      tGTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGGCTTTCAg   >  2:808559/1‑93 (MQ=255)
                                                                                       gTATCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGc  >  1:1762390/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCATTGTTCGGCGGGAATTTTGTCTTTCAGC  >  NC_000913/442300‑442479

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: