Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,921,515 C→T T179T (ACG→ACA)  atpB ← F0 sector of membrane‑bound ATP synthase, subunit a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,921,5150CT100.0% 246.5 / NA 75T179T (ACG→ACA) atpBF0 sector of membrane‑bound ATP synthase, subunit a
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/39);  total (36/39)

ACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGCGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGA  >  NC_000913/3921398‑3921647
                                                                                                                     |                                                                                                                                    
aCCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAgccg                                                                                                                 >  2:1147107/1‑139 (MQ=255)
   ggCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAGCGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGccgccg                                                                                                              <  2:547294/139‑1 (MQ=255)
       tACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGc                                                                                                          >  1:1644464/1‑139 (MQ=255)
        aCATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGcc                                                                                                         <  2:405295/139‑1 (MQ=255)
           tGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCttt                                                                                                      <  1:1147107/139‑1 (MQ=255)
              tACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCAt                                                                                                   >  2:199045/1‑139 (MQ=255)
               aCCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCAtt                                                                                                  >  1:1423755/1‑139 (MQ=255)
                   aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              >  2:462121/1‑139 (MQ=255)
                   aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              <  1:156344/139‑1 (MQ=255)
                      aGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAgccgcc                                                                                                               <  2:972279/119‑1 (MQ=255)
                      aGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAgccgcc                                                                                                               >  1:972279/1‑119 (MQ=255)
                             aaCCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGaa                                                                                    <  1:1496008/139‑1 (MQ=255)
                                  aGTGAAATTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAAc                                                                                   <  1:860602/135‑1 (MQ=255)
                                  aGTGAAATTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAAc                                                                                   >  2:860602/1‑135 (MQ=255)
                                  aGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCTCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGaa                                                                               <  2:1184378/139‑1 (MQ=255)
                                  aGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGaa                                                                               >  2:122645/1‑139 (MQ=255)
                                   gTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGACTAAGTTGACAGGAAAGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGGAGCCGCCGATGCCTTTCATTTTGATACTGTAGAACAGAAt                                                                              <  1:1007268/139‑1 (MQ=255)
                                       aaCTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAgg                                                                          <  1:173069/139‑1 (MQ=255)
                                        aCTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGa                                                                         <  2:41398/139‑1 (MQ=255)
                                        aCTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGa                                                                         <  1:462121/139‑1 (MQ=255)
                                         cTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGa                                                                                                             >  2:372216/1‑102 (MQ=255)
                                         cTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGa                                                                                                             <  1:372216/102‑1 (MQ=255)
                                            gTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAActct                                                                                                                             >  2:650084/1‑83 (MQ=255)
                                            gTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAActct                                                                                                                             <  1:650084/83‑1 (MQ=255)
                                             tttGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGaa                                                                                                                     <  1:165232/90‑1 (MQ=255)
                                             tttGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGaa                                                                                                                     >  2:165232/1‑90 (MQ=255)
                                             tttGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAAt                                                                    >  2:1258855/1‑139 (MQ=255)
                                                ggACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACg                                                                 >  2:1651632/1‑139 (MQ=255)
                                                    agcatGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  1:871932/139‑1 (MQ=255)
                                                    agcagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:373485/139‑1 (MQ=255)
                                                    agcagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTCTAGAACAGAATCAGGATAAATACGCCCa                                                             <  1:1027120/139‑1 (MQ=255)
                                                      cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATg                                                                                            <  2:200859/106‑1 (MQ=255)
                                                      cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATg                                                                                            >  1:200859/1‑106 (MQ=255)
                                                      cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGcc                                                               <  1:1214681/135‑1 (MQ=255)
                                                      cagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGcc                                                               >  2:1214681/1‑135 (MQ=255)
                                                           ttACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCAt                                                      >  1:772538/1‑139 (MQ=255)
                                                           ttACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCAt                                                      >  1:404587/1‑139 (MQ=255)
                                                              ccccTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  1:298644/1‑129 (MQ=255)
                                                              ccccTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:298644/129‑1 (MQ=255)
                                                                   tCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCttt                                                                                                      >  2:494230/1‑83 (MQ=255)
                                                                   tCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCttt                                                                                                      <  1:494230/83‑1 (MQ=255)
                                                                   tCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAgg                                                                          <  2:395054/111‑1 (MQ=255)
                                                                   tCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAgg                                                                          >  1:395054/1‑111 (MQ=255)
                                                                   tCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCg                                              <  1:515177/139‑1 (MQ=255)
                                                                       ggATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTAc                                          >  1:133280/1‑139 (MQ=255)
                                                                            aaGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCa                                     >  1:838628/1‑139 (MQ=255)
                                                                                  aCAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              >  2:1124485/1‑76 (MQ=255)
                                                                                  aCAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              <  1:1124485/76‑1 (MQ=255)
                                                                                  aCAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCg                               <  2:133280/139‑1 (MQ=255)
                                                                                          gAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTa                                                                                       >  2:1547042/1‑75 (MQ=255)
                                                                                          gAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTa                                                                                       <  1:1547042/75‑1 (MQ=255)
                                                                                             cGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTg                                                                                         >  1:543075/1‑70 (MQ=255)
                                                                                             cGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTg                                                                                         <  2:543075/70‑1 (MQ=255)
                                                                                                ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAgg                                                                          >  1:1043549/1‑82 (MQ=255)
                                                                                                ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAgg                                                                          <  2:1043549/82‑1 (MQ=255)
                                                                                                ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACcaca                 >  2:1610326/1‑139 (MQ=255)
                                                                                                ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACcaca                 >  1:1564662/1‑139 (MQ=255)
                                                                                                 cAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:667218/94‑1 (MQ=255)
                                                                                                 cAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  1:667218/1‑94 (MQ=255)
                                                                                                 cAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCacac                <  1:1258855/139‑1 (MQ=255)
                                                                                                 cAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTGCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCacac                <  1:1398281/139‑1 (MQ=255)
                                                                                                    tGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAgccgc                                                                                                                <  1:508445/40‑1 (MQ=255)
                                                                                                    tGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAgccgc                                                                                                                >  2:508445/1‑40 (MQ=255)
                                                                                                          aaCGGCTGCAGTGTCAAGTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  2:947542/1‑85 (MQ=255)
                                                                                                          aaCGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  1:947542/85‑1 (MQ=255)
                                                                                                          aaCGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGtt                                       <  1:1335792/107‑1 (MQ=255)
                                                                                                          aaCGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGtt                                       >  2:1335792/1‑107 (MQ=255)
                                                                                                            cGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggc     >  1:59213/1‑139 (MQ=255)
                                                                                                            cGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACACCCAGTGCCATAGACAGCGTTACGTTTACGTCCGCAGACGGAACCACACGc              <  1:1012621/130‑1 (MQ=255)
                                                                                                            cGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACACCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGc              >  2:1012621/1‑130 (MQ=255)
                                                                                                              gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATg                                                                                            >  2:879678/1‑50 (MQ=255)
                                                                                                              gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATg                                                                                            <  1:879678/50‑1 (MQ=255)
                                                                                                               cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              <  2:1603449/47‑1 (MQ=255)
                                                                                                               cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              >  1:1603449/1‑47 (MQ=255)
                                                                                                               cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGa  <  1:1610326/139‑1 (MQ=255)
                                                                                                                     |                                                                                                                                    
ACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGCGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGA  >  NC_000913/3921398‑3921647

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