Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,506,318 (C)5→6 coding (1134/1227 nt) yhfW ← putative mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,506,3131.C88.8% 778.4 / 90.5 259coding (1139/1227 nt)yhfWputative mutase
Reads supporting (aligned to +/- strand):  ref base . (11/18);  new base C (122/108);  total (133/126)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-02

Only 100 of 282 total aligned reads displayed.

AAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGG  >  NC_000913/3506221‑3506405
                                                                                             |                                                                                            
aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                               >  2:3109582/1‑93 (MQ=255)
  aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑cc                                                                                            <  1:5273433/93‑1 (MQ=255)
  aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑cc                                                                                            <  1:3348926/93‑1 (MQ=255)
   cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑ccc                                                                                           <  2:6130636/93‑1 (MQ=255)
    tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑cccc                                                                                          >  1:7294986/1‑93 (MQ=255)
    tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑cccc                                                                                          >  2:3540505/1‑93 (MQ=255)
      aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑cccccc                                                                                        <  2:705145/93‑2 (MQ=255)
        aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑ccccccac                                                                                      >  1:5602345/1‑90 (MQ=255)
         gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑ccccccaca                                                                                     >  1:3125587/1‑89 (MQ=255)
           gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATc                                                                                    <  2:5023521/93‑1 (MQ=255)
           gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCa                                                                                   <  2:361096/93‑1 (MQ=255)
           gtgtCTACTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATc                                                                                    <  1:1442912/93‑1 (MQ=255)
             gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAg                                                                                  >  2:6687251/1‑93 (MQ=255)
              tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGa                                                                                 >  1:7476357/1‑93 (MQ=255)
               ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGaa                                                                                >  2:7299130/1‑93 (MQ=255)
                 ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAg                                                                              >  2:4442473/1‑93 (MQ=255)
                 ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAg                                                                              <  2:970318/93‑1 (MQ=255)
                  cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGc                                                                             >  2:3004793/1‑93 (MQ=255)
                    gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgt                                                                           <  1:3630258/93‑1 (MQ=255)
                    gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGCAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCgtg                                                                          <  1:6212664/93‑1 (MQ=255)
                     cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtg                                                                          <  1:1123595/93‑1 (MQ=255)
                     cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtg                                                                          <  2:3843457/93‑1 (MQ=255)
                     cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtg                                                                          >  2:4129768/1‑93 (MQ=255)
                      aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtgg                                                                         <  1:6064742/93‑1 (MQ=255)
                      aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtgg                                                                         <  1:4573148/93‑1 (MQ=255)
                         ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTgc                                                                      <  2:4172071/93‑1 (MQ=255)
                              gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                 <  2:4609128/93‑1 (MQ=255)
                              gggaggAAAGAAAACAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCacac                                                                >  2:7490571/1‑93 (MQ=255)
                               ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCacac                                                                <  2:2396346/93‑1 (MQ=255)
                                gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACAcc                                                               >  1:3804741/1‑93 (MQ=255)
                                 aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                             <  1:6596335/93‑1 (MQ=255)
                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                             <  2:6271832/93‑1 (MQ=255)
                                    aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                           <  1:4450905/93‑1 (MQ=255)
                                    aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGAGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                          <  2:4976141/93‑1 (MQ=255)
                                     aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                          <  1:4535085/93‑1 (MQ=255)
                                      agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                         >  1:2961944/1‑93 (MQ=255)
                                      agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                         >  1:7151890/1‑93 (MQ=255)
                                       gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                        <  2:1745908/93‑1 (MQ=255)
                                         aaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                      >  2:1945428/1‑93 (MQ=255)
                                         aaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                      >  2:3388423/1‑93 (MQ=255)
                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                     >  2:1217097/1‑93 (MQ=255)
                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                    <  2:2354911/93‑1 (MQ=255)
                                           agagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                    >  1:3684460/1‑93 (MQ=255)
                                           agagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                    >  2:6556145/1‑93 (MQ=255)
                                            gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                   >  1:6480263/1‑93 (MQ=255)
                                            gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                  >  1:1860676/1‑93 (MQ=255)
                                             agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                  <  1:2255697/93‑1 (MQ=255)
                                                gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                               >  2:6336217/1‑93 (MQ=255)
                                                gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                               >  2:4067195/1‑93 (MQ=255)
                                                gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                               >  2:1543200/1‑93 (MQ=255)
                                                gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                               >  1:7362907/1‑93 (MQ=255)
                                                gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                               >  1:5203371/1‑93 (MQ=255)
                                                 aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                              >  1:94769/1‑93 (MQ=255)
                                                 aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                              >  1:4022655/1‑93 (MQ=255)
                                                   cATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcatc                                            <  1:5706737/93‑1 (MQ=255)
                                                     ttttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCGCATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                          >  1:4315444/1‑93 (MQ=255)
                                                     ttttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                          <  2:5793404/93‑1 (MQ=255)
                                                     ttttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                          >  2:2316721/1‑93 (MQ=255)
                                                     ttttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                          <  2:6230717/93‑1 (MQ=255)
                                                      tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                        <  1:6432415/93‑1 (MQ=255)
                                                       ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                        <  1:1879398/93‑1 (MQ=255)
                                                       ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                        <  1:2074856/93‑1 (MQ=255)
                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                      >  2:937359/1‑93 (MQ=255)
                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                      >  1:5535783/1‑93 (MQ=255)
                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                      <  2:5237115/93‑1 (MQ=255)
                                                           ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                    <  2:561767/93‑1 (MQ=255)
                                                            gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                   <  2:2619727/93‑1 (MQ=255)
                                                            gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                   <  2:5542438/93‑1 (MQ=255)
                                                            gcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGa                                   >  1:5055822/1‑93 (MQ=255)
                                                              ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                 >  1:6254577/1‑93 (MQ=255)
                                                               gtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                >  2:1671229/1‑93 (MQ=255)
                                                                  gcgcgcGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                             <  2:6901850/93‑1 (MQ=255)
                                                                   cgcgcgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                            <  2:3950709/93‑1 (MQ=255)
                                                                     cgcgGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                         <  1:2845008/93‑1 (MQ=255)
                                                                      gcgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                         <  1:4866685/93‑1 (MQ=255)
                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                        <  2:2777609/93‑1 (MQ=255)
                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                        <  2:4711518/93‑1 (MQ=255)
                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                        <  2:2462171/93‑1 (MQ=255)
                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                        >  2:2253519/1‑93 (MQ=255)
                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                        <  1:6847499/93‑1 (MQ=255)
                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                        >  1:67040/1‑93 (MQ=255)
                                                                        ggAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                       >  2:5340092/1‑93 (MQ=255)
                                                                         gAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                      >  2:5911001/1‑93 (MQ=255)
                                                                           aaaaTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc                    >  2:2063859/1‑93 (MQ=255)
                                                                           aaaaTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc                    >  1:533711/1‑93 (MQ=255)
                                                                            aaaTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                   <  2:3249175/93‑1 (MQ=255)
                                                                             aaTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac                  >  1:2379592/1‑93 (MQ=255)
                                                                              aTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacc                 >  1:4766819/1‑93 (MQ=255)
                                                                              aTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacc                 >  2:4824006/1‑93 (MQ=255)
                                                                               ttCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                >  1:1499406/1‑93 (MQ=255)
                                                                                tCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcaccac               <  2:3813202/93‑1 (MQ=255)
                                                                                 cacacaCGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt              >  1:2618088/1‑93 (MQ=255)
                                                                                 cacacaCGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt              <  1:7317679/93‑1 (MQ=255)
                                                                                   cacacGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc            >  1:3526082/1‑93 (MQ=255)
                                                                                    acacGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg           >  2:7096988/1‑93 (MQ=255)
                                                                                    acacGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc          <  1:1045071/93‑1 (MQ=255)
                                                                                      acGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg         >  1:5907103/1‑93 (MQ=255)
                                                                                        ggTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg       <  1:1742838/93‑1 (MQ=255)
                                                                                        ggTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg       <  1:589198/93‑1 (MQ=255)
                                                                                            g‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGTGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg  <  2:2660064/93‑1 (MQ=255)
                                                                                             |                                                                                            
AAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGG  >  NC_000913/3506221‑3506405

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