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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,432,9881.C99.3% 1212.9 / ‑0.2 292coding (573/1125 nt)smfprotein Smf
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (161/129);  total (163/129)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 321 total aligned reads displayed.

CCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATT  >  NC_000913/3432888‑3433083
                                                                                                     |                                                                                                
ccaccaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                   >  2:6376336/1‑101 (MQ=255)
  accaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑cc                                                                                               >  1:6376341/1‑101 (MQ=255)
     aGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑ccccc                                                                                            >  2:6376343/1‑101 (MQ=255)
     aGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑ccccc                                                                                            >  1:6376347/1‑101 (MQ=255)
        acacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑ccccccct                                                                                         >  2:6376350/1‑99 (MQ=255)
          acCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑ccccccctgt                                                                                       >  2:6376354/1‑97 (MQ=255)
          acCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑ccccccctgt                                                                                       >  1:6376356/1‑97 (MQ=255)
           ccTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑ccccccctgtt                                                                                      <  2:6376212/101‑6 (MQ=255)
           ccTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑ccccccctgtt                                                                                      >  2:6376358/1‑96 (MQ=255)
             tttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCa                                                                                     >  1:6376365/1‑101 (MQ=255)
              ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCaa                                                                                    >  2:6376367/1‑101 (MQ=255)
              ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCaa                                                                                    <  1:6375765/101‑1 (MQ=255)
              ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCaa                                                                                    <  2:6375844/101‑1 (MQ=255)
               tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTc                                                                                      >  1:6376369/1‑98 (MQ=255)
                 cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCa                                                                                 >  1:6376372/1‑101 (MQ=255)
                 cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCa                                                                                 <  1:6375925/101‑1 (MQ=255)
                  ttAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAg                                                                                <  2:6375977/101‑1 (MQ=255)
                   tAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGa                                                                               >  2:6376376/1‑101 (MQ=255)
                    aGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGAc                                                                              >  2:6376379/1‑101 (MQ=255)
                    aGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGAc                                                                              <  1:6376066/101‑1 (MQ=255)
                       ccACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTgg                                                                           >  1:6376382/1‑101 (MQ=255)
                        cACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGc                                                                          >  2:6376387/1‑101 (MQ=255)
                        cACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGc                                                                          >  1:6376386/1‑101 (MQ=255)
                        cACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGc                                                                          <  2:6376096/101‑1 (MQ=255)
                         aCTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCa                                                                         <  1:6375872/101‑1 (MQ=255)
                           tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGAc                                                                              >  1:6376392/1‑94 (MQ=255)
                              tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                    >  2:6376398/1‑101 (MQ=255)
                              tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                    >  2:6376403/1‑101 (MQ=255)
                              tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                    >  1:6376397/1‑101 (MQ=255)
                              tAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                    <  2:6375608/101‑1 (MQ=255)
                                aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTc                                                                  >  2:6376406/1‑101 (MQ=255)
                                aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTc                                                                  >  1:6376407/1‑101 (MQ=255)
                                 tGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                 <  1:6375917/101‑1 (MQ=255)
                                 tGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCg                                                                 <  2:6376044/101‑1 (MQ=255)
                                  gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAg                                                                        <  1:6375959/93‑1 (MQ=255)
                                  gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCgg                                                                <  1:6376089/101‑1 (MQ=255)
                                  gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCgg                                                                >  1:6376410/1‑101 (MQ=255)
                                    gATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGc                                                              >  2:6376413/1‑101 (MQ=255)
                                    gATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGc                                                              <  1:6376329/101‑1 (MQ=255)
                                       ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                            >  1:6376418/1‑100 (MQ=255)
                                       ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATg                                                           >  2:6376417/1‑101 (MQ=255)
                                        tcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGa                                                          <  1:6375830/101‑1 (MQ=255)
                                             gTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                     >  1:6376426/1‑101 (MQ=255)
                                             gTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                     >  2:6376427/1‑101 (MQ=255)
                                              tGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcgg                                                    >  1:6376430/1‑101 (MQ=255)
                                                ggAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTc                                                                  >  2:6376434/1‑85 (MQ=255)
                                                 gAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGa                                                 <  2:6376021/101‑1 (MQ=255)
                                                  aaaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGt                                                                   >  2:6376437/1‑82 (MQ=255)
                                                  aaaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                <  2:6375944/101‑1 (MQ=255)
                                                     ttGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                             <  1:6376142/101‑1 (MQ=255)
                                                     ttGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                             <  2:6375866/101‑1 (MQ=255)
                                                      tGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTc                                                                  >  1:6376446/1‑79 (MQ=255)
                                                      tGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                  >  2:6376445/1‑95 (MQ=255)
                                                      tGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                            >  1:6376441/1‑101 (MQ=255)
                                                      tGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                            <  2:6375855/101‑1 (MQ=255)
                                                       gTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGGAGCCAGTCGGGCATGACGGCGGGGAt                                                <  2:6375906/96‑1 (MQ=255)
                                                         aaGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGt                                         >  1:6376450/1‑101 (MQ=255)
                                                             aaGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTa                                     >  1:6376455/1‑101 (MQ=255)
                                                             aaGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTa                                     >  1:6376458/1‑101 (MQ=255)
                                                             aaGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTa                                     <  1:6375992/101‑1 (MQ=255)
                                                              aGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTaa                                    >  2:6376459/1‑101 (MQ=255)
                                                                gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGt                                  >  1:6376462/1‑101 (MQ=255)
                                                                 ggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                 <  2:6375873/101‑1 (MQ=255)
                                                                 ggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                 <  2:6375964/101‑1 (MQ=255)
                                                                 ggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                 <  1:6376033/101‑1 (MQ=255)
                                                                    ggAACATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCAt                              >  1:6376465/1‑101 (MQ=255)
                                                                      aaCATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                            >  1:6376470/1‑101 (MQ=255)
                                                                      aaCATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                            <  2:6375991/101‑1 (MQ=255)
                                                                       aCATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTcc                                >  2:6376475/1‑96 (MQ=255)
                                                                       aCATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTcc                           <  2:6375936/101‑1 (MQ=255)
                                                                        cATCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                          >  1:6376476/1‑101 (MQ=255)
                                                                         aTCGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTcccc                         >  1:6376479/1‑101 (MQ=255)
                                                                           cGAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCgggg                                                  >  2:6376484/1‑74 (MQ=255)
                                                                            gAGGGGAAATTCCGAGACGAGAGCGCCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                       >  1:6376488/1‑101 (MQ=255)
                                                                            gAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATg                                                           <  2:6375774/64‑1 (MQ=255)
                                                                            gAGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAAt                      >  1:6376485/1‑101 (MQ=255)
                                                                             aGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATa                     <  2:6375981/101‑1 (MQ=255)
                                                                             aGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATa                     <  2:6375918/101‑1 (MQ=255)
                                                                             aGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATa                     >  2:6376491/1‑101 (MQ=255)
                                                                             aGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATa                     <  1:6375954/101‑1 (MQ=255)
                                                                             aGGGGAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATa                     <  1:6376062/101‑1 (MQ=255)
                                                                              ggggAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATAc                    <  1:6375741/101‑1 (MQ=255)
                                                                                ggAAATTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAg                  <  1:6376112/101‑1 (MQ=255)
                                                                                   aaTTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGt                                  <  1:6376025/82‑1 (MQ=255)
                                                                                   aaTTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCaa               <  2:6376156/101‑1 (MQ=255)
                                                                                   aaTTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCACCCAGTCGGGCATGACGGCGGGGAt                                                >  2:6376497/1‑68 (MQ=255)
                                                                                    aTTCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAAt              <  1:6375705/101‑1 (MQ=255)
                                                                                     ttCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATg             <  1:6376195/101‑1 (MQ=255)
                                                                                      tCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                              >  1:6376500/1‑67 (MQ=255)
                                                                                      tCCGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGc            <  2:6375957/101‑1 (MQ=255)
                                                                                       ccGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCt           <  1:6375941/101‑1 (MQ=255)
                                                                                       ccGAGACGAGAGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCt           >  2:6376501/1‑101 (MQ=255)
                                                                                         gagacgagaGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCa                >  2:6376506/1‑94 (MQ=255)
                                                                                          agacgagaGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGAc        <  2:6375891/101‑1 (MQ=255)
                                                                                           gacgagaGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTTGGGCATGAcggcg                                                     >  2:6376512/1‑55 (MQ=255)
                                                                                           gacgagaGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg       >  2:6376509/1‑101 (MQ=255)
                                                                                              gagaGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATgg                                          >  2:6376518/1‑63 (MQ=255)
                                                                                              gagaGCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGc            >  1:6376514/1‑93 (MQ=255)
                                                                                                gagCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAtt  >  1:6376523/1‑101 (MQ=255)
                                                                                                gagCGC‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAtt  <  1:6376183/101‑1 (MQ=255)
                                                                                                     |                                                                                                
CCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATT  >  NC_000913/3432888‑3433083

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