Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,257,952 T→C D703G (GAT→GGT)  tdcE ← putative formate acetyltransferase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,257,9520TC99.5% 627.4 / ‑4.9 202D703G (GAT→GGT) tdcEputative formate acetyltransferase 3
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (107/94);  total (107/95)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 210 total aligned reads displayed.

AGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATT  >  NC_000913/3257862‑3258041
                                                                                          |                                                                                         
    tAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCtt                                                                                     >  1:1872777/1‑93 (MQ=255)
    tAGGATATATTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCtt                                                                                     >  2:4569830/1‑93 (MQ=255)
     aGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTc                                                                                    <  2:3789922/93‑1 (MQ=255)
     aGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTc                                                                                    >  1:200842/1‑93 (MQ=255)
       gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCgt                                                                                  >  1:1640424/1‑93 (MQ=255)
       gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCgt                                                                                  <  1:2971150/93‑1 (MQ=255)
         tatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCgtgg                                                                                >  2:4117195/1‑93 (MQ=255)
         tatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCgtgg                                                                                >  2:1058697/1‑93 (MQ=255)
          atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCgtggt                                                                               >  2:1012004/1‑93 (MQ=255)
           tttttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCgtggtg                                                                              <  2:3521346/93‑1 (MQ=255)
           tttttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCgtggtg                                                                              >  1:915957/1‑93 (MQ=255)
               tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGaaa                                                                          <  2:693018/93‑1 (MQ=255)
               tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGaaa                                                                          <  1:949846/93‑1 (MQ=255)
               tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGaaa                                                                          <  1:3827331/93‑1 (MQ=255)
                ccGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAAt                                                                         <  2:3043143/93‑1 (MQ=255)
                 cGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATa                                                                        <  1:2683110/93‑1 (MQ=255)
                  gggTTCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATAc                                                                       <  1:3195363/93‑1 (MQ=255)
                  gggTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATAc                                                                       >  1:4133065/1‑93 (MQ=255)
                   ggTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATAcc                                                                      >  2:1282808/1‑93 (MQ=255)
                     tGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCa                                                                    <  2:3487381/93‑1 (MQ=255)
                     tGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCa                                                                    >  1:3345351/1‑93 (MQ=255)
                       cTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATc                                                                  <  2:219142/93‑1 (MQ=255)
                        tCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATcc                                                                 >  1:3249257/1‑93 (MQ=255)
                         cGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCca                                                                <  1:1751791/93‑1 (MQ=255)
                           aTGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCcagc                                                              >  2:2897545/1‑93 (MQ=255)
                            tGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCcagca                                                             <  2:647508/93‑1 (MQ=255)
                            tGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCcagca                                                             >  1:1819006/1‑93 (MQ=255)
                              gCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGTCCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAgg                                                           <  2:3124665/93‑1 (MQ=255)
                              gCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAgg                                                           <  2:3206191/93‑1 (MQ=255)
                                aTCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGcc                                                         >  2:2370924/1‑93 (MQ=255)
                                 tCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCg                                                        >  2:4249616/1‑93 (MQ=255)
                                 tCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCg                                                        <  2:2628522/93‑1 (MQ=255)
                                     gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACaa                                                    >  1:718026/1‑93 (MQ=255)
                                     gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGCTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACaa                                                    >  2:2283752/1‑93 (MQ=255)
                                      cagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAg                                                   <  1:1519078/93‑1 (MQ=255)
                                        gcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggt                                                 <  2:3160593/93‑1 (MQ=255)
                                         caTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggtt                                                >  1:868422/1‑93 (MQ=255)
                                          aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttg                                               <  1:4330688/93‑1 (MQ=255)
                                          aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttg                                               <  1:2908242/93‑1 (MQ=255)
                                          aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttg                                               >  2:1469914/1‑93 (MQ=255)
                                          aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttg                                               >  2:1372080/1‑93 (MQ=255)
                                          aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttg                                               >  2:4670157/1‑93 (MQ=255)
                                          aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttg                                               <  1:1070937/93‑1 (MQ=255)
                                           tttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttgg                                              >  1:3466450/1‑93 (MQ=255)
                                           tttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttgg                                              >  1:1589344/1‑93 (MQ=255)
                                           tttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttgg                                              >  1:4029885/1‑93 (MQ=255)
                                            ttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttggt                                             <  2:1162426/93‑1 (MQ=255)
                                             tCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttggtt                                            <  2:3897496/93‑1 (MQ=255)
                                             tCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAggttggtt                                            >  1:4491282/1‑93 (MQ=255)
                                              cccGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGttt                                           <  2:3006942/93‑1 (MQ=255)
                                               ccGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt                                          <  2:2683917/93‑1 (MQ=255)
                                                cGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa                                         >  1:3674088/1‑93 (MQ=255)
                                                cGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa                                         <  1:2230196/93‑1 (MQ=255)
                                                  aTTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACg                                       <  1:4672982/93‑1 (MQ=255)
                                                  aTTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACg                                       >  2:3470685/1‑93 (MQ=255)
                                                  aTTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCCACAAGGTTGGTTTTACg                                       >  1:389808/1‑93 (MQ=255)
                                                   ttCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGt                                      >  1:529654/1‑93 (MQ=255)
                                                   ttCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGt                                      >  1:1428583/1‑93 (MQ=255)
                                                    tCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTa                                     >  1:4554810/1‑93 (MQ=255)
                                                    tCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTa                                     <  2:3593906/93‑1 (MQ=255)
                                                    tCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTa                                     <  1:3309597/93‑1 (MQ=255)
                                                     cATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc                                    >  1:425541/1‑93 (MQ=255)
                                                     cATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc                                    >  1:2492123/1‑93 (MQ=255)
                                                     cATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGCTGGTTTTACGTAc                                    >  1:54960/1‑93 (MQ=255)
                                                     cATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAAAACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc                                    >  1:3186114/1‑93 (MQ=255)
                                                      aTTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACt                                   >  1:2780340/1‑93 (MQ=255)
                                                      aTTACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACt                                   >  1:101177/1‑93 (MQ=255)
                                                        tACGTTGACGTTGAGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTgg                                 <  1:4032542/93‑1 (MQ=255)
                                                          cgttgacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAt                               <  1:1225021/93‑1 (MQ=255)
                                                           gttgacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtc                              >  2:4021366/1‑93 (MQ=255)
                                                           gttgacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtc                              >  1:2976394/1‑93 (MQ=255)
                                                            ttgacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtct                             >  2:1573783/1‑93 (MQ=255)
                                                            ttgacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCATGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtct                             <  1:56515/93‑1 (MQ=255)
                                                              gacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtcttc                           <  2:4190961/93‑1 (MQ=255)
                                                              gacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtcttc                           <  2:1935308/93‑1 (MQ=255)
                                                              gacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtcttc                           <  1:3728796/93‑1 (MQ=255)
                                                              gacgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGAtcttc                           >  2:1451914/1‑93 (MQ=255)
                                                                cgttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATcttctt                         <  2:3698008/93‑1 (MQ=255)
                                                                 gttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCttt                        <  2:1798495/93‑1 (MQ=255)
                                                                 gttgaGGTGTTGACCGCCTTCGACACCCCCTTCGTGGTGGAAATACCCATCCACCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCttt                        >  2:3309137/1‑93 (MQ=255)
                                                                  ttgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTg                       >  2:3341073/1‑93 (MQ=255)
                                                                   tgaGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGc                      >  2:554576/1‑93 (MQ=255)
                                                                     aGGTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGccc                    <  1:3878782/93‑1 (MQ=255)
                                                                      ggTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCa                   <  1:208141/93‑1 (MQ=255)
                                                                      ggTGTTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCa                   >  2:3658175/1‑93 (MQ=255)
                                                                       gtgtTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAg                  <  2:2553852/93‑1 (MQ=255)
                                                                        tgtTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgc                 <  2:3958045/93‑1 (MQ=255)
                                                                         gtTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                <  2:337871/93‑1 (MQ=255)
                                                                         gtTGACCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                >  2:2229008/1‑93 (MQ=255)
                                                                           tGACCGCCTGCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGcgcc              <  1:4376846/93‑1 (MQ=255)
                                                                             aCCGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGcgccgc            >  1:2250667/1‑93 (MQ=255)
                                                                               cGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAg          <  1:984370/93‑1 (MQ=255)
                                                                               cGCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAg          >  1:2369312/1‑93 (MQ=255)
                                                                                gCCTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAgg         >  2:1701444/1‑93 (MQ=255)
                                                                                  cTTCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGaa       <  2:4310480/93‑1 (MQ=255)
                                                                                   ttCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAAc      >  2:3104633/1‑93 (MQ=255)
                                                                                   ttCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAAc      >  2:2083622/1‑93 (MQ=255)
                                                                                   ttCGACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTATGCCCAGCGCCGCAGGAAc      <  2:4543801/93‑1 (MQ=255)
                                                                                      gACACCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGAt   >  1:2466648/1‑93 (MQ=255)
                                                                                       acacCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGAtt  <  2:4657497/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATT  >  NC_000913/3257862‑3258041

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