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.8% 278.8 / ‑3.9 81coding (573/1125 nt)smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (35/45);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCG  >  NC_000913/3432856‑3433106
                                                                                                                                     |                                                                                                                        
cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccc‑c                                                                                                                     >  2:193383/1‑137 (MQ=255)
cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccc‑c                                                                                                                     <  1:193383/137‑1 (MQ=255)
         aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAACTGTAAGCAAGGGGTGGAACATCGAGGGGAACTTCCGAGACGAGAGCG‑CCC‑CCACTGTTCAAg                                                                                                          >  1:729356/1‑139 (MQ=255)
           cACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCa                                                                                                        >  1:679772/1‑139 (MQ=255)
                   aaaGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCa                                                                                                <  2:1285255/139‑1 (MQ=255)
                   aaaGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCGCC‑CTGTTCAAGCAGACTGGCa                                                                                                >  2:317956/1‑139 (MQ=255)
                     aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGc                                                                                              >  2:846385/1‑139 (MQ=255)
                       ccGCTT‑CACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAg                                                                                           <  2:396856/139‑1 (MQ=255)
                               accaccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCa                                                                                    >  1:657058/1‑139 (MQ=255)
                                     aGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                              >  1:544894/1‑139 (MQ=255)
                                     aGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                              <  2:248095/139‑1 (MQ=255)
                                       tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                            >  2:1180056/1‑139 (MQ=255)
                                       tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                            <  2:647485/139‑1 (MQ=255)
                                       tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                            <  1:273596/139‑1 (MQ=255)
                                          acCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         <  2:679772/139‑1 (MQ=255)
                                           ccTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGa                                                                        >  1:1007958/1‑139 (MQ=255)
                                            cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                              <  2:600390/132‑1 (MQ=255)
                                            cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                              >  1:600390/1‑132 (MQ=255)
                                            cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                       <  2:935685/139‑1 (MQ=255)
                                              ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                     >  2:1245819/1‑139 (MQ=255)
                                              ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                     >  2:424024/1‑139 (MQ=255)
                                              ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATGGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                     >  2:432938/1‑139 (MQ=255)
                                               tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                    <  2:544894/139‑1 (MQ=255)
                                               tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                    <  2:939136/139‑1 (MQ=255)
                                               tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                    <  2:729356/139‑1 (MQ=255)
                                                aCTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                   <  1:1245819/139‑1 (MQ=255)
                                                aCTTAGACCACCGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                   <  1:846385/139‑1 (MQ=255)
                                                 cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATg                                                                  <  2:657058/139‑1 (MQ=255)
                                                      aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         >  1:445933/1‑127 (MQ=255)
                                                      aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         <  2:445933/127‑1 (MQ=255)
                                                      aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                             >  1:56659/1‑139 (MQ=255)
                                                        cACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTaa                                                           <  1:1145058/139‑1 (MQ=255)
                                                          cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccc‑cc‑cctg                                                                                                               >  2:783537/1‑81 (MQ=255)
                                                          cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccc‑cc‑cctg                                                                                                               <  1:783537/84‑4 (MQ=255)
                                                           tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAg                                                                                                          <  1:1210630/89‑1 (MQ=255)
                                                           tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAg                                                                                                          >  2:1210630/1‑89 (MQ=255)
                                                           tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                        >  1:426355/1‑139 (MQ=255)
                                                             aTAATGCGATTTCTTCGTTGGAAATTGTAAGCAAGGGGTTGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      <  2:1081247/139‑1 (MQ=255)
                                                             aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      <  2:56659/139‑1 (MQ=255)
                                                             aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      <  1:1138556/139‑1 (MQ=255)
                                                             aTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      <  2:1007958/139‑1 (MQ=255)
                                                                     aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      >  1:11028/1‑131 (MQ=255)
                                                                     aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      <  2:11028/131‑1 (MQ=255)
                                                                     aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                              >  1:1072224/1‑139 (MQ=255)
                                                                       ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCa                                                                                                        <  2:352409/79‑1 (MQ=255)
                                                                       ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCa                                                                                                        >  1:352409/1‑79 (MQ=255)
                                                                       ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                    >  1:1231976/1‑115 (MQ=255)
                                                                       ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                    <  2:1231976/115‑1 (MQ=255)
                                                                         cttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGAc                                                                                                     >  2:19854/1‑80 (MQ=255)
                                                                         cttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGAc                                                                                                     <  1:19854/80‑1 (MQ=255)
                                                                         cttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACa                                          >  1:1190073/1‑139 (MQ=255)
                                                                            cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                 <  1:755615/129‑1 (MQ=255)
                                                                            cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                 >  2:755615/1‑129 (MQ=255)
                                                                                   aaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                <  2:1190073/139‑1 (MQ=255)
                                                                                     ttGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                              <  1:432938/139‑1 (MQ=255)
                                                                                       gTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                        >  1:330105/1‑111 (MQ=255)
                                                                                       gTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                        <  2:330105/111‑1 (MQ=255)
                                                                                       gTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGcc                            <  2:426355/139‑1 (MQ=255)
                                                                                        tAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                        >  1:925470/1‑110 (MQ=255)
                                                                                        tAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                        <  2:925470/110‑1 (MQ=255)
                                                                                          aGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAtt                         <  2:416317/139‑1 (MQ=255)
                                                                                               ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         <  2:215333/86‑1 (MQ=255)
                                                                                               ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         >  1:215333/1‑86 (MQ=255)
                                                                                               ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCa                           <  2:1124618/132‑1 (MQ=255)
                                                                                               ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCa                           >  1:1124618/1‑132 (MQ=255)
                                                                                                 ggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGt                  >  1:1096926/1‑139 (MQ=255)
                                                                                                     gAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      <  1:484959/99‑1 (MQ=255)
                                                                                                     gAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      >  2:484959/1‑99 (MQ=255)
                                                                                                      aaCATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         >  1:541910/1‑79 (MQ=255)
                                                                                                      aaCATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         <  2:541910/79‑1 (MQ=255)
                                                                                                          tCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgctgct         >  1:844201/1‑139 (MQ=255)
                                                                                                           cGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         <  2:316121/74‑1 (MQ=255)
                                                                                                           cGAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                         >  1:316121/1‑74 (MQ=255)
                                                                                                            gAGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGTATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt       <  2:1096926/139‑1 (MQ=255)
                                                                                                             aGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGt                                                                                          >  1:326234/1‑55 (MQ=255)
                                                                                                             aGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGt                                                                                          <  2:326234/55‑1 (MQ=255)
                                                                                                             aGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                       <  1:97447/74‑1 (MQ=255)
                                                                                                             aGGGGAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                       >  2:97447/1‑74 (MQ=255)
                                                                                                               gggAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGAGTGGCAGCCAGTCGGGCATGACGGCGGGGATGTATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATg    <  2:773689/139‑1 (MQ=255)
                                                                                                                 gAAATTCCGAGACGAGAGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATgcg  <  1:424024/139‑1 (MQ=255)
                                                                                                                         gagacgagaGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                        <  1:328062/109‑1 (MQ=255)
                                                                                                                         gagacgagaGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                        >  2:328062/1‑109 (MQ=255)
                                                                                                                           gacgagaGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGt                                                                >  2:512166/1‑67 (MQ=255)
                                                                                                                           gacgagaGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGt                                                                <  1:512166/67‑1 (MQ=255)
                                                                                                                           gacgagaGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      <  1:1284739/77‑1 (MQ=255)
                                                                                                                           gacgagaGCGCCCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                      >  2:1284739/1‑77 (MQ=255)
                                                                                                                                     |                                                                                                                        
CACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCC‑CC‑CTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCG  >  NC_000913/3432856‑3433106

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