Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,955,459 T→C T517A (ACC→GCC)  torZ ← trimethylamine N‑oxide reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,955,4590TC98.4% 218.0 / NA 63T517A (ACC→GCC) torZtrimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (27/35);  minor base A (1/0);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTA  >  NC_000913/1955323‑1955595
                                                                                                                                        |                                                                                                                                        
tCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCAt                                                                                                                                       <  2:157264/140‑1 (MQ=255)
          tctcTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCa                                                                                                                             <  1:189541/140‑1 (MQ=255)
          tctcTTTTCCGCCAGGTTTGAGAAATTCCGCCAGACCGGCAAACACGTGAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCa                                                                                                                             <  2:343138/140‑1 (MQ=255)
              ttttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGTTGGTTGCTGTATTCGCCGGCCATCGTCAAGTCATTgc                                                                                                                         >  1:275721/1‑140 (MQ=255)
              ttttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTgc                                                                                                                         >  1:180200/1‑140 (MQ=255)
              ttttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTgc                                                                                                                         >  1:828735/1‑140 (MQ=255)
              ttttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGATGTAATCACCGGACATCGTCGAGTCATTgc                                                                                                                         >  1:540072/1‑140 (MQ=255)
               tttCCGCCAGGTTTGCGTAATTCCGTAAGCTCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTgcg                                                                                                                        <  1:315504/140‑1 (MQ=255)
                  ccGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTc                                                                                                                     <  1:835694/140‑1 (MQ=255)
                   cGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCa                                                                                                                    >  2:565471/1‑140 (MQ=255)
                   cGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCa                                                                                                                    <  1:659100/140‑1 (MQ=255)
                     ccAGGTTTGAGTAATTCAGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTCGCTATAATCACCGGCCATCGTCAAGTCATTGCGCTCaaa                                                                                                                  >  2:340810/1‑140 (MQ=255)
                     ccAGGTTTGAGTAATGCCGCAAGATAGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCaaa                                                                                                                  <  1:363971/140‑1 (MQ=255)
                        ggTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGa                                                                                                               >  2:584047/1‑140 (MQ=255)
                               gTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATgtggtg                                                                                                        <  1:809182/140‑1 (MQ=255)
                                tAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGa                                                                                                       <  1:449927/140‑1 (MQ=255)
                                  aTTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACg                                                                                                                <  2:640134/129‑1 (MQ=255)
                                  aTTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACg                                                                                                                >  1:640134/1‑129 (MQ=255)
                                   ttCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCg                                                                                                    >  2:734042/1‑140 (MQ=255)
                                     ccGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGt                                                                                                  <  1:739888/140‑1 (MQ=255)
                                          aGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGAt                                                                                                      <  2:38769/131‑1 (MQ=255)
                                          aGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGAt                                                                                                      >  1:38769/1‑131 (MQ=255)
                                           gATCGGCAAACATGTCAAAATCGTTACGCGCTCCAATTTGCGGAGCAACAGCCTGCTTCATTGCCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCTAGTCACTGCGCTCAAACGATGTGGTGATCGGTAATACg                                                                                            <  2:540072/140‑1 (MQ=255)
                                           gATCGGCAAACACGTCACAATCTTTACGCGCTTCAAATTCCGGAGCAATAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACg                                                                                            <  2:275721/140‑1 (MQ=255)
                                                gCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATc                                                                                       >  1:774120/1‑140 (MQ=255)
                                                 cAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCt                                                                                      <  1:29025/140‑1 (MQ=255)
                                                        gTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGCTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGtt                                                                               >  2:810116/1‑140 (MQ=255)
                                                            aaaTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGgc                                                                           <  2:36453/140‑1 (MQ=255)
                                                            aaaTCGTTACGCACTTCGAATTCCGGAGCAACAGCATGCATCATCGGCTCAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGgc                                                                           <  2:253421/140‑1 (MQ=255)
                                                                 gTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGctgctg                                                                      <  2:851305/140‑1 (MQ=255)
                                                                 gTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGctgctg                                                                      <  1:584047/140‑1 (MQ=255)
                                                                    aCGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGctgctg                                                                      >  2:51399/1‑137 (MQ=255)
                                                                    aCGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGctgctg                                                                      <  1:51399/137‑1 (MQ=255)
                                                                         cTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAg                                                              <  1:565471/140‑1 (MQ=255)
                                                                          ttCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATTTGGTGATCCGTAATACAATATCTGCATGTTCGGCTGCTGCGGTCCAGt                                                             >  2:295375/1‑140 (MQ=255)
                                                                                    ggAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGa                                                   >  2:84590/1‑140 (MQ=255)
                                                                                    ggAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGATAATACGAAATCGGCGGGTTTGGCTGCTGCGGTCCAGTAGCATTCAGa                                                   >  2:269333/1‑140 (MQ=255)
                                                                                           cAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCTTGTATGGTTGCTGCGGCCCAGTAGCATTCAGAAacgacg                                            >  2:297139/1‑140 (MQ=255)
                                                                                                tGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGatcat                                       <  2:180200/140‑1 (MQ=255)
                                                                                                 gCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGatcatc                                      >  1:722445/1‑140 (MQ=255)
                                                                                                  cTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCAtct                                     >  1:509458/1‑140 (MQ=255)
                                                                                                  cTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTACGTGTTTGGCTGGTGCGGTCCAGTAGCATTCAGAAACGACGATCAtct                                     >  2:274904/1‑140 (MQ=255)
                                                                                                         ggCACAATATGCTGGCTGCTGTAATCACCGGCCATCGGCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGAATTCAGAAACGACGATCATCTCCGGttt                              <  1:340810/140‑1 (MQ=255)
                                                                                                              gatatGCTGGTTGCTGTAATCATCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGcc                         <  1:418690/139‑1 (MQ=255)
                                                                                                              aaTATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGTTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGcc                         <  1:4785/140‑1 (MQ=255)
                                                                                                              aaTATGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGcc                         <  2:774120/140‑1 (MQ=255)
                                                                                                                tatGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGcc                         >  1:519798/1‑138 (MQ=255)
                                                                                                                tatGCTGGTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGcc                         <  2:519798/138‑1 (MQ=255)
                                                                                                                      ggTTGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCACTCAGAAACGACGATCAACTCCCGTTTCTGCCATGCTTTa                 >  2:333394/1‑140 (MQ=255)
                                                                                                                        ttGCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCCTACTTTAAt               >  1:108116/1‑140 (MQ=255)
                                                                                                                          gCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCa             >  1:108034/1‑140 (MQ=255)
                                                                                                                          gCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCa             >  1:565155/1‑140 (MQ=255)
                                                                                                                          gCTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGAAGCATTCACAAAGGACACTCATCTCAGGTTTGTGCCATGGTTTAATcc             >  2:418690/1‑139 (MQ=255)
                                                                                                                           cTGTAATCACCGGCCGTCGTCAGGTTATTTCGCTCAAGCGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAg            <  1:297139/140‑1 (MQ=255)
                                                                                                                           cTGTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAg            <  2:108034/140‑1 (MQ=255)
                                                                                                                             gTAATCTCCCGCCATCGTCAAGTCATTGCGCCCAAACGATGTGTTTATCGGTAACACGGTATTTGCCTGTTTGGCTGCTGCGGTCCAGTCGCATTCAAAAACGTCGCTCAACTCCGGTTTCTGCCATGCTTTAATCAGAc          <  2:108116/140‑1 (MQ=255)
                                                                                                                             gTAATCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGAc          <  1:810116/140‑1 (MQ=255)
                                                                                                                             acagtCACCGGCCATAGTCAAGTCATTGCTCTCAAACAATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGAc          <  1:420249/136‑1 (MQ=255)
                                                                                                                                aTCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGTTTTCTGCCATGCTTTAATCAGACGAt       <  2:722445/140‑1 (MQ=255)
                                                                                                                                aTCACCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGAt       <  1:734042/140‑1 (MQ=255)
                                                                                                                                   aCCGGCCATCGTCAATTCATTGTGCTCAAACGATGTGTTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAg    <  1:314824/140‑1 (MQ=255)
                                                                                                                                   aCCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAg    <  2:828735/140‑1 (MQ=255)
                                                                                                                                   aCCGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAAAAACGACGTTCATCTCGGGTTTCTGACATGCTTTAATCAGACGATTAg    >  2:314824/1‑140 (MQ=255)
                                                                                                                                   aCCGGCCATCGTCAAGTCATTGCGCCCAAGCGATGTGCTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAg    <  1:295375/140‑1 (MQ=255)
                                                                                                                                     cGGCCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTa  >  2:676646/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTA  >  NC_000913/1955323‑1955595

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