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.C97.4% 261.1 / ‑0.6 78coding (573/1125 nt)smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (36/40);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01
*NC_0009133,432,9882.C97.4% 263.3 / ‑0.6 78coding (573/1125 nt)smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (36/40);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGA  >  NC_000913/3432850‑3433118
                                                                                                                                           ||                                                                                                                                     
tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                                         >  2:1512629/1‑139 (MQ=255)
tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                                         >  2:1645254/1‑139 (MQ=255)
      cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑ccccc                                                                                                                              >  1:955468/1‑139 (MQ=255)
      cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑ccccc                                                                                                                              >  2:1748347/1‑139 (MQ=255)
        ccAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑ccccccc                                                                                                                            >  2:1314958/1‑137 (MQ=255)
            cGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑ccccccctgtt                                                                                                                        >  1:1416708/1‑133 (MQ=255)
             gAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑ccccccctgttc                                                                                                                       >  1:898904/1‑132 (MQ=255)
                    tACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGa                                                                                                                  <  1:1645254/139‑1 (MQ=255)
                      cGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACt                                                                                                                >  2:1441552/1‑139 (MQ=255)
                       gCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTg                                                                                                               <  2:750304/139‑1 (MQ=255)
                          aaGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCC‑‑CCCCCTGTTCAAGCAGACTGGc                                                                                                             <  2:772771/139‑1 (MQ=255)
                           aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGCACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCCGACTGGCAg                                                                                                           >  1:1334423/1‑139 (MQ=255)
                           aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAg                                                                                                           >  2:819108/1‑139 (MQ=255)
                            gCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGc                                                                                                          <  1:694558/139‑1 (MQ=255)
                              cGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCa                                                                                                        >  1:1641529/1‑139 (MQ=255)
                               gCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                                       <  2:1334423/139‑1 (MQ=255)
                                    caccaccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCggg                                                                                                  <  1:1512629/139‑1 (MQ=255)
                                    caccaccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAACCAGACTGGCAGCCAGTCggg                                                                                                  <  1:1079434/139‑1 (MQ=255)
                                         ccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGa                                                                                             >  2:365877/1‑139 (MQ=255)
                                           aGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcg                                                                                           <  1:1441552/139‑1 (MQ=255)
                                            gTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                                          <  2:955468/139‑1 (MQ=255)
                                              acacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                                        <  1:668429/139‑1 (MQ=255)
                                               cacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcgg                                                                                       >  2:1407547/1‑139 (MQ=255)
                                                   tttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                                   >  2:1712645/1‑139 (MQ=255)
                                                    ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                                       <  2:1436948/118‑1 (MQ=255)
                                                    ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                                       >  1:1436948/1‑118 (MQ=255)
                                                    ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATg                                                                                  >  1:1611132/1‑139 (MQ=255)
                                                    ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATg                                                                                  >  2:969532/1‑139 (MQ=255)
                                                     tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGCGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                                 <  1:1314958/139‑1 (MQ=255)
                                                     tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                                 <  1:969532/139‑1 (MQ=255)
                                                       cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                               <  2:1416708/139‑1 (MQ=255)
                                                       cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                               <  1:1407547/139‑1 (MQ=255)
                                                         tAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCa                                                                                                                    <  2:340625/100‑1 (MQ=255)
                                                         tAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCa                                                                                                                    >  1:340625/1‑100 (MQ=255)
                                                          aGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑cccc                                                                                                                               >  2:1726605/1‑86 (MQ=255)
                                                          aGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑cccc                                                                                                                               <  1:1726605/86‑1 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg‑‑‑c                                                                                                                                     <  1:761584/80‑1 (MQ=255)
                                                            aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg‑‑‑c                                                                                                                                     >  2:761584/1‑80 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGACATTCCGAGACGAGAGCG‑‑‑c‑‑ccccccct                                                                                                                           >  1:1162214/1‑81 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                                      >  2:857756/1‑139 (MQ=255)
                                                                cTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑ccccccct                                                                                                                           <  2:1162214/84‑4 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                               >  1:960249/1‑129 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                               <  2:960249/129‑1 (MQ=255)
                                                                 tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGt                                                                     >  1:1802532/1‑139 (MQ=255)
                                                                  gATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAAGGTATTAATTc                                                                    <  1:1064989/139‑1 (MQ=255)
                                                                       tGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                                                               <  2:1641529/139‑1 (MQ=255)
                                                                         cGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                               >  2:1028941/1‑105 (MQ=255)
                                                                         cGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                               <  1:1028941/105‑1 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCCGACTGGCAGCCAg                                                                                                       >  1:945254/1‑95 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                                       <  2:945254/95‑1 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCa                                                           >  2:1277972/1‑139 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGCGGATGAATGGTATTAAGTCCATTCCCCa                                                           >  1:1140032/1‑139 (MQ=255)
                                                                           aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑CCTCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCa                                                           >  2:853530/1‑139 (MQ=255)
                                                                             ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAg                                                                                                                      <  2:1365557/78‑1 (MQ=255)
                                                                             ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAg                                                                                                                      >  1:1365557/1‑78 (MQ=255)
                                                                                 tcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                  >  1:24876/1‑126 (MQ=255)
                                                                                 tcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                  <  2:24876/126‑1 (MQ=255)
                                                                                   gTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCa                                                   >  2:1046562/1‑139 (MQ=255)
                                                                                     gggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑c                                                                                                                                  <  1:443876/56‑1 (MQ=255)
                                                                                     gggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑c‑‑c                                                                                                                                  >  2:443876/1‑56 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                                 >  2:1375586/1‑102 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                                 <  1:1375586/102‑1 (MQ=255)
                                                                                          aTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                            <  2:898904/139‑1 (MQ=255)
                                                                                                aGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTGAGTCCATTCCCCAATACAGCAATGCTGACGCCAt                                      >  1:508867/1‑139 (MQ=255)
                                                                                                        gTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGt                              <  2:1021849/139‑1 (MQ=255)
                                                                                                            aaCATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                            <  1:668622/121‑1 (MQ=255)
                                                                                                            aaCATCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                            >  2:668622/1‑121 (MQ=255)
                                                                                                               aTCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                  >  1:423887/1‑96 (MQ=255)
                                                                                                               aTCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGAGGGGATGAATGGTATTAAGTCCa                                                                  <  2:423887/96‑1 (MQ=255)
                                                                                                                tCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                                   >  1:449441/1‑78 (MQ=255)
                                                                                                                tCGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                                   <  2:449441/78‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGt                                                                                                      >  1:77351/1‑58 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGt                                                                                                      <  2:77351/58‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                                                               >  2:110316/1‑97 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                                                               <  1:110316/97‑1 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgctgct                     >  2:452218/1‑139 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgctgct                     >  2:1559243/1‑139 (MQ=255)
                                                                                                                 cGAGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgctgct                     <  2:1802532/139‑1 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                               <  1:345578/79‑1 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                               >  2:345578/1‑79 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTACGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                                    >  1:933301/1‑122 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                                    <  2:933301/122‑1 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt                   >  1:1034152/1‑139 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt                   >  2:847977/1‑139 (MQ=255)
                                                                                                                   aGGGGAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt                   <  1:857756/139‑1 (MQ=255)
                                                                                                                    ggggAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTa                  <  2:1611132/139‑1 (MQ=255)
                                                                                                                     gggAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                  >  2:701523/1‑90 (MQ=255)
                                                                                                                     gggAAATTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                                  <  1:701523/90‑1 (MQ=255)
                                                                                                                          aTTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTGAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCt            <  2:508867/139‑1 (MQ=255)
                                                                                                                          aTTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAg                                                     >  1:1521714/1‑98 (MQ=255)
                                                                                                                          aTTCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAg                                                     <  2:1521714/98‑1 (MQ=255)
                                                                                                                           ttCCGAGACGAGAGCGCC‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTa           <  1:847977/139‑1 (MQ=255)
                                                                                                                                  gagagcGCC‑‑‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGa  <  1:1466349/137‑1 (MQ=255)
                                                                                                                                           ||                                                                                                                                     
TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑‑C‑‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGA  >  NC_000913/3432850‑3433118

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