Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,432,994 (C)6→7 coding (567/1125 nt) smf ← DNA recombination‑mediator A family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,432,9881.C98.5% 234.2 / ‑3.6 67coding (573/1125 nt)smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (32/34);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTAC  >  NC_000913/3432850‑3433110
                                                                                                                                           |                                                                                                                           
tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                              >  1:224603/1‑139 (MQ=255)
      cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGcccccc                                                                                                                      <  1:257516/139‑1 (MQ=255)
      cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGcccccc                                                                                                                      >  2:52677/1‑139 (MQ=255)
               aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAg                                                                                                              >  1:1320385/1‑139 (MQ=255)
               aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccccctgtt                                                                                                                 <  2:1138429/135‑6 (MQ=255)
               aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccccctgtt                                                                                                                 >  1:1138429/1‑130 (MQ=255)
                           aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGc                                                                                                  >  1:542532/1‑139 (MQ=255)
                           aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGc                                                                                                  <  2:542532/139‑1 (MQ=255)
                                    caccaccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGc                                                                                         <  1:589641/139‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                <  2:224603/139‑1 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                >  1:726688/1‑139 (MQ=255)
                                             tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                <  2:1294120/139‑1 (MQ=255)
                                               agtaaaggtgtAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGa                                                                                                          >  2:731988/11‑111 (MQ=255)
                                               agtaaaggtgtAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGa                                                                                                          <  1:731988/101‑1 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                                 >  2:900833/1‑139 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGATGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             >  2:835014/1‑123 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                                  <  2:1179987/118‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                                  >  1:1179987/1‑118 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             >  2:624312/1‑123 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             <  1:624312/123‑1 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             <  1:835014/123‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcgc                                                                                                                            >  2:1182781/1‑75 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcgc                                                                                                                            <  1:1182781/75‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             >  1:25698/1‑122 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                             <  2:25698/122‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                            <  1:52677/139‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                            >  2:1345155/1‑139 (MQ=255)
                                                                   aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:1045800/139‑1 (MQ=255)
                                                                   aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:726688/139‑1 (MQ=255)
                                                                   aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  1:900833/139‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCa                                                                                                            <  2:610810/84‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCa                                                                                                            >  1:610810/1‑84 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       <  1:1415244/121‑1 (MQ=255)
                                                                        gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       >  2:1415244/1‑121 (MQ=255)
                                                                             ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATa                                                >  2:807500/1‑139 (MQ=255)
                                                                                  cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcgc                                                                                                                            <  1:1257673/58‑1 (MQ=255)
                                                                                  cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcgc                                                                                                                            >  2:1257673/1‑58 (MQ=255)
                                                                                      ggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGc                                       >  1:557475/1‑139 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       >  2:1021282/1‑103 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       <  1:1021282/103‑1 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     >  2:975466/1‑121 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     <  1:975466/121‑1 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    >  1:905926/1‑138 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    <  2:905926/138‑1 (MQ=255)
                                                                                              tAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:1157687/112‑1 (MQ=255)
                                                                                              tAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          >  1:1157687/1‑112 (MQ=255)
                                                                                                aGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       <  2:625221/97‑1 (MQ=255)
                                                                                                aGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                       >  1:625221/1‑97 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGAc                                                                                                         <  1:772881/58‑1 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGAc                                                                                                         >  2:772881/1‑58 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                <  1:548762/83‑1 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                >  2:548762/1‑83 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCt                        >  2:634458/1‑139 (MQ=255)
                                                                                                     ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCt                        >  1:545508/1‑139 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                            <  2:666340/134‑1 (MQ=255)
                                                                                                      gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                            >  1:666340/1‑134 (MQ=255)
                                                                                                          ggAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAg                   <  1:634458/139‑1 (MQ=255)
                                                                                                                tCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          >  2:690409/1‑94 (MQ=255)
                                                                                                                tCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  1:690409/94‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATg                                                                      <  2:502022/81‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATg                                                                      >  1:502022/1‑81 (MQ=255)
                                                                                                                  gAGGGGAAATTCCGGGACGAGTGCGCACCCCCTGTTCAAACCGACTGGCAGGCAGTAGGGCATCAGGGCGGGGATCCATGGTACTAGGTCCATTCCCCAATACAGCAATGCGGACGCCAGTTAGGTGTTAGGCTGCttt           >  2:880685/1‑139 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                                            <  2:1103740/57‑1 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                                            >  1:1103740/1‑57 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     <  2:893804/96‑1 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                     >  1:893804/1‑96 (MQ=255)
                                                                                                                     gggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATg        <  1:807500/139‑1 (MQ=255)
                                                                                                                      ggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                  <  1:152122/112‑1 (MQ=255)
                                                                                                                      ggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                  >  2:152122/1‑112 (MQ=255)
                                                                                                                           ttCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTac  <  2:557475/139‑1 (MQ=255)
                                                                                                                              cgagacgagaGCGCCCCCCCTGTTCAAGCAGACTGGc                                                                                                     <  2:652981/37‑1 (MQ=255)
                                                                                                                              cgagacgagaGCGCCCCCCCTGTTCAAGCAGACTGGc                                                                                                     >  1:652981/1‑37 (MQ=255)
                                                                                                                                agacgagaGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGACGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    <  1:171646/100‑1 (MQ=255)
                                                                                                                                agacgagaGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGACGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                    >  2:171646/1‑100 (MQ=255)
                                                                                                                                      gagCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          >  1:1351370/1‑72 (MQ=255)
                                                                                                                                      gagCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                          <  2:1351370/72‑1 (MQ=255)
                                                                                                                                           |                                                                                                                           
TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTAC  >  NC_000913/3432850‑3433110

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