Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,432,994 (C)6→8 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.C100.0% 330.4 / NA 92G191G (GGC→GGGsmfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (44/48);  total (44/48)
*NC_0009133,432,9882.C96.7% 322.4 / 1.7 92coding (573/1125 nt)smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (43/46);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATG  >  NC_000913/3432854‑3433104
                                                                                                                                       ||                                                                                                                    
gTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cccc                                                                                                                  <  2:538505/139‑1 (MQ=255)
  cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cccccc                                                                                                                <  1:148365/139‑1 (MQ=255)
      aGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cccccccctg                                                                                                            >  2:1304254/1‑135 (MQ=255)
           aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cccccccctgttcaa                                                                                                       >  1:996927/1‑130 (MQ=255)
            ccACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAg                                                                                                        <  1:1249657/139‑1 (MQ=255)
                    cAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTgg                                                                                                >  2:591701/1‑139 (MQ=255)
                     caaGCCGCTTCCACCACCAGTAAACCTTTACTTAGACCACTGATAATGCGGTTTCTTCGTGGGAAATTGTAAGCAAGGGGGGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGc                                                                                               <  2:1424954/138‑1 (MQ=255)
                          cGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCa                                                                                          >  2:1291005/1‑139 (MQ=255)
                           gCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                         <  1:591701/139‑1 (MQ=255)
                            cTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGt                                                                                        >  1:1471409/1‑139 (MQ=255)
                             ttCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTc                                                                                       <  2:1450/139‑1 (MQ=255)
                                         tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggc                                                                           >  2:538616/1‑139 (MQ=255)
                                          acacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcgg                                                                         <  1:1592118/139‑1 (MQ=255)
                                          acacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTCGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                          <  1:985867/139‑1 (MQ=255)
                                             ccTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                       <  1:497577/139‑1 (MQ=255)
                                             ccTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                       >  2:497577/1‑139 (MQ=255)
                                              cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGACGGGCATGACGGCGGGGa                                                                      >  2:1082017/1‑139 (MQ=255)
                                              cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cc                                                                                                                    >  1:497252/1‑91 (MQ=255)
                                              cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cc                                                                                                                    <  2:497252/91‑1 (MQ=255)
                                                ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATg                                                                    >  2:1302819/1‑139 (MQ=255)
                                                   cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 >  2:872506/1‑139 (MQ=255)
                                                   cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 <  1:1388083/139‑1 (MQ=255)
                                                   cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 >  2:1126408/1‑139 (MQ=255)
                                                   cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 <  1:872506/139‑1 (MQ=255)
                                                    ttAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATg                                                                <  2:127406/139‑1 (MQ=255)
                                                        aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                 <  2:282846/118‑1 (MQ=255)
                                                        aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                 >  1:282846/1‑118 (MQ=255)
                                                        aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAt                                                            >  2:555238/1‑139 (MQ=255)
                                                        aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTGCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAt                                                            >  1:927456/1‑139 (MQ=255)
                                                           aCTGTTAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTgg                                                                                                >  1:937536/1‑100 (MQ=255)
                                                           aCTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTgg                                                                                                <  2:937536/100‑1 (MQ=255)
                                                            cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAg                                                                                             >  2:180300/1‑102 (MQ=255)
                                                            cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAg                                                                                             <  1:180300/102‑1 (MQ=255)
                                                            cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                 <  1:1087697/114‑1 (MQ=255)
                                                            cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                 >  2:1087697/1‑114 (MQ=255)
                                                             tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGTAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                          <  2:1055564/120‑1 (MQ=255)
                                                             tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                  <  2:1596115/128‑1 (MQ=255)
                                                             tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                  >  1:1596115/1‑128 (MQ=255)
                                                             tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGt                                                       >  2:1423895/1‑139 (MQ=255)
                                                             tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGt                                                       >  2:1345790/1‑139 (MQ=255)
                                                             tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGt                                                       >  1:363968/1‑139 (MQ=255)
                                                             tGATAATCCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGTAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                          >  1:1055564/1‑120 (MQ=255)
                                                              gATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                      <  2:1114815/139‑1 (MQ=255)
                                                                tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  1:538616/139‑1 (MQ=255)
                                                                tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  1:1126408/139‑1 (MQ=255)
                                                                tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  1:1490821/139‑1 (MQ=255)
                                                                tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  1:356704/139‑1 (MQ=255)
                                                                tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  2:1471409/139‑1 (MQ=255)
                                                                  aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                        <  1:89046/69‑1 (MQ=255)
                                                                  aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                        >  2:89046/1‑69 (MQ=255)
                                                                  aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGNCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCAtt                                                  >  2:1228050/1‑139 (MQ=255)
                                                                  aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCAtt                                                  >  1:1395966/1‑139 (MQ=255)
                                                                   tGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                                                 <  1:1304254/139‑1 (MQ=255)
                                                                   tGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                                                 <  2:656821/139‑1 (MQ=255)
                                                                   tGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTcc                                                >  1:1032407/1‑139 (MQ=255)
                                                                    gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTCGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                               <  2:1032407/139‑1 (MQ=255)
                                                                     cGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                               <  1:1345790/139‑1 (MQ=255)
                                                                        tttCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            <  2:363968/139‑1 (MQ=255)
                                                                         ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                       <  2:525495/111‑1 (MQ=255)
                                                                         ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                                       >  1:525495/1‑111 (MQ=255)
                                                                              cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                                      <  2:316861/91‑1 (MQ=255)
                                                                              cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                                      >  1:316861/1‑91 (MQ=255)
                                                                              cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                  >  2:571108/1‑111 (MQ=255)
                                                                              cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                  <  1:571108/111‑1 (MQ=255)
                                                                                    aaaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCt                                >  1:683234/1‑139 (MQ=255)
                                                                                      aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cccccccc                                                                                                              >  2:940095/1‑55 (MQ=255)
                                                                                      aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cccccccc                                                                                                              <  1:940095/57‑3 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                         <  2:707586/69‑1 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                         >  1:707586/1‑69 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 >  2:1444764/1‑93 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 <  1:1444764/93‑1 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  2:535637/106‑1 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    >  1:535637/1‑106 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  2:897515/106‑1 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    >  1:897515/1‑106 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                              <  1:132280/128‑1 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                              >  2:132280/1‑128 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                      >  1:1134396/1‑136 (MQ=255)
                                                                                                 ggggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                      <  2:1134396/136‑1 (MQ=255)
                                                                                                    gTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                        >  1:102156/1‑35 (MQ=255)
                                                                                                      ggAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTaa              >  2:1195937/1‑139 (MQ=255)
                                                                                                        aaCATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCTTTTACCTGTAAgg            <  1:1082017/139‑1 (MQ=255)
                                                                                                           aTCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAg                                                                                                        <  1:1298714/44‑1 (MQ=255)
                                                                                                           aTCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAg                                                                                                        >  2:1298714/1‑44 (MQ=255)
                                                                                                            tCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            >  1:1038109/1‑103 (MQ=255)
                                                                                                            tCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            <  2:1038109/103‑1 (MQ=255)
                                                                                                             cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 >  1:1548822/1‑81 (MQ=255)
                                                                                                             cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                               <  2:274731/99‑1 (MQ=255)
                                                                                                             cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                               >  1:274731/1‑99 (MQ=255)
                                                                                                             cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgctgct       <  2:683234/139‑1 (MQ=255)
                                                                                                               aGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                          <  1:1190744/70‑1 (MQ=255)
                                                                                                               aGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                          >  2:1190744/1‑70 (MQ=255)
                                                                                                               aGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt     <  1:1302819/139‑1 (MQ=255)
                                                                                                                 gggAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 <  1:1080677/77‑1 (MQ=255)
                                                                                                                 gggAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 >  2:1080677/1‑77 (MQ=255)
                                                                                                                 gggAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                      <  2:559479/88‑1 (MQ=255)
                                                                                                                 gggAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                      >  1:559479/1‑88 (MQ=255)
                                                                                                                  ggAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATg  >  2:1388833/1‑139 (MQ=255)
                                                                                                                              acgagaGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                          >  2:809658/1‑55 (MQ=255)
                                                                                                                              acgagaGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                          <  1:809658/55‑1 (MQ=255)
                                                                                                                                       ||                                                                                                                    
GTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATG  >  NC_000913/3432854‑3433104

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