Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,917,686 A→C 5.1% S194A (TCA→GCA)  atpG ← F1 sector of membrane‑bound ATP synthase, gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,917,6860AC 5.1% 292.2 / 5.5 79S194A (TCA→GCA) atpGF1 sector of membrane‑bound ATP synthase, gamma subunit
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGA  >  NC_000913/3917543‑3917833
                                                                                                                                               |                                                                                                                                                   
ggcCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGcc                                                                                                                                                 <  4:783849/148‑1 (MQ=255)
                 ggTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAgctg                                                                                                                                >  4:155676/1‑148 (MQ=255)
                       cAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGAGCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGAt                                                                                                                          <  3:155676/148‑1 (MQ=255)
                         aCCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATgg                                                                                                                        >  4:915202/1‑148 (MQ=255)
                          ccACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGCTGCCGTTACCGgcat                                                                                                                                      <  4:928869/133‑2 (MQ=255)
                          ccACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGCTGCCGTTACCGgcat                                                                                                                                      >  3:928869/1‑132 (MQ=255)
                                    aTAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGc                                                                                                                                                  >  3:597230/1‑111 (MQ=255)
                                    aTAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGc                                                                                                                                                  <  4:597230/111‑1 (MQ=255)
                                          cTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGCTGGTCGGAACCTg                                                                                                              <  3:81718/141‑1 (MQ=255)
                                          cTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTg                                                                                                              >  4:81718/1‑141 (MQ=255)
                                            gagaTTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGgt                                                                                                     >  4:816316/1‑148 (MQ=255)
                                            gagaTTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGgt                                                                                                     >  1:145488/1‑148 (MQ=255)
                                                ttCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGt                                                                                                                       >  3:477098/1‑126 (MQ=255)
                                                ttCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGt                                                                                                                       <  4:477098/126‑1 (MQ=255)
                                                 tCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                         <  2:843058/107‑1 (MQ=255)
                                                 tCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                         >  1:843058/1‑107 (MQ=255)
                                                       tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACATGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                         <  2:328295/101‑1 (MQ=255)
                                                       tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGa                                                                                                                   >  2:955416/1‑123 (MQ=255)
                                                       tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGa                                                                                                                   <  1:955416/123‑1 (MQ=255)
                                                       tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGAGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                         >  1:328295/1‑101 (MQ=255)
                                                        aaCGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAAttt                                                                                         <  2:145488/148‑1 (MQ=255)
                                                          cgacgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACa                                                                                                         <  1:1097846/130‑1 (MQ=255)
                                                          cgacgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACa                                                                                                         >  2:1097846/1‑130 (MQ=255)
                                                          cgacgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATttgt                                                                                       >  3:304645/1‑148 (MQ=255)
                                                            acgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAcc                                                                                                                >  3:179153/1‑121 (MQ=255)
                                                            acgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAcc                                                                                                                <  4:179153/121‑1 (MQ=255)
                                                                    gggTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTtaataa                                                                                             >  2:572447/1‑132 (MQ=255)
                                                                    gggTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTtaataa                                                                                             <  1:572447/132‑1 (MQ=255)
                                                                         tCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTaaa                                                                        <  4:304645/148‑1 (MQ=255)
                                                                          ccAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAc                                                                                                                 >  4:147793/1‑106 (MQ=255)
                                                                          ccAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAc                                                                                                                 <  3:147793/106‑1 (MQ=255)
                                                                          ccAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAg                                                                       <  3:816316/148‑1 (MQ=255)
                                                                             gCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGAt                                                                                                                          >  2:323981/1‑94 (MQ=255)
                                                                             gCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGAt                                                                                                                          <  1:323981/94‑1 (MQ=255)
                                                                             gCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCtt                                                                    <  3:915202/148‑1 (MQ=255)
                                                                              cAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTg                                                                   >  3:31076/1‑148 (MQ=255)
                                                                                   ccTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAg                                                              <  4:31076/148‑1 (MQ=255)
                                                                                    cTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGAt                                                                                                                          <  1:1127569/87‑1 (MQ=255)
                                                                                    cTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGAt                                                                                                                          >  2:1127569/1‑87 (MQ=255)
                                                                                     ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                               >  1:815852/1‑129 (MQ=255)
                                                                                     ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                               <  2:815852/129‑1 (MQ=255)
                                                                                     ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                               <  4:750086/129‑1 (MQ=255)
                                                                                     ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                               >  3:750086/1‑129 (MQ=255)
                                                                                           tCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCtt                                                      >  1:688494/1‑148 (MQ=255)
                                                                                                  cGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCtt                                                      >  3:482404/1‑141 (MQ=255)
                                                                                                  cGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCtt                                                      <  4:482404/141‑1 (MQ=255)
                                                                                                    tACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGCTGCCGTTACCGgca                                                                                                                                       >  4:81372/1‑58 (MQ=38)
                                                                                                    tACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGCTGCCGTTACCGgca                                                                                                                                       <  3:81372/58‑1 (MQ=38)
                                                                                                        ggTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACAt                                                                                                        >  1:209402/1‑85 (MQ=255)
                                                                                                        ggTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACAt                                                                                                        <  2:209402/85‑1 (MQ=255)
                                                                                                          tAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGc                                                         >  1:285214/1‑130 (MQ=255)
                                                                                                          tAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGc                                                         <  2:285214/130‑1 (MQ=255)
                                                                                                               ccAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATttgtt                                                                                      <  2:793576/96‑1 (MQ=255)
                                                                                                               ccAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATttgtt                                                                                      >  1:793576/1‑96 (MQ=255)
                                                                                                                 aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGa                                                                                                                           <  3:78462/57‑1 (MQ=255)
                                                                                                                 aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGa                                                                                                                           >  4:78462/1‑57 (MQ=255)
                                                                                                                 aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATgg                                                                                                      >  2:471342/1‑78 (MQ=255)
                                                                                                                 aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATgg                                                                                                      <  1:471342/78‑1 (MQ=255)
                                                                                                                 aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgt                                                   >  1:284550/1‑129 (MQ=255)
                                                                                                                 aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgt                                                   <  2:284550/129‑1 (MQ=255)
                                                                                                                  ggATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgtcg                                                 >  3:555770/1‑130 (MQ=255)
                                                                                                                  ggATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgtcg                                                 <  4:555770/130‑1 (MQ=255)
                                                                                                                  ggATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAgg                                             <  2:526297/134‑1 (MQ=255)
                                                                                                                  ggATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAgg                                             >  1:526297/1‑134 (MQ=255)
                                                                                                                  ggATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACt                               <  2:688494/148‑1 (MQ=255)
                                                                                                                              cAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCg                   >  1:769103/1‑148 (MQ=255)
                                                                                                                                 atcatcatcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGa                  >  1:156279/1‑146 (MQ=255)
                                                                                                                                 atcatcatcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGa                  <  2:156279/146‑1 (MQ=255)
                                                                                                                                  tcatcatcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACttt                             <  1:1102469/134‑1 (MQ=255)
                                                                                                                                  tcatcatcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACttt                             >  2:1102469/1‑134 (MQ=255)
                                                                                                                                    atcatcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGt             <  2:769103/148‑1 (MQ=255)
                                                                                                                                      catcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGa                  >  2:402661/1‑141 (MQ=255)
                                                                                                                                      catcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGa                  <  1:402661/141‑1 (MQ=255)
                                                                                                                                       atcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTc                                                                 >  2:209080/1‑93 (MQ=255)
                                                                                                                                       atcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTc                                                                 <  1:209080/93‑1 (MQ=255)
                                                                                                                                         catcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgt                                                   <  1:654613/105‑1 (MQ=255)
                                                                                                                                         catcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgt                                                   >  2:654613/1‑105 (MQ=255)
                                                                                                                                             tGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCgg                       <  4:139590/129‑1 (MQ=255)
                                                                                                                                             tGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCgg                       >  3:139590/1‑129 (MQ=255)
                                                                                                                                               aTGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGa  >  3:129313/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCCTGCTCGCTGGCCAGGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGA  >  NC_000913/3917543‑3917833

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