Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,917,692 G→T 5.1% P192T (CCG→ACG)  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,6920GT 5.1% 279.1 / 5.1 78P192T (CCG→ACG) atpGF1 sector of membrane‑bound ATP synthase, gamma subunit
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base T (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

GGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAG  >  NC_000913/3917560‑3917835
                                                                                                                                    |                                                                                                                                               
ggTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAgctg                                                                                                                                  >  4:155676/1‑148 (MQ=255)
      cAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGAGCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGAt                                                                                                                            <  3:155676/148‑1 (MQ=255)
        aCCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATgg                                                                                                                          >  4:915202/1‑148 (MQ=255)
         ccACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGCTGCCGTTACCGgcat                                                                                                                                        >  3:928869/1‑132 (MQ=255)
         ccACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGCTGCCGTTACCGgcat                                                                                                                                        <  4:928869/133‑2 (MQ=255)
                         cTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGCTGGTCGGAACCTg                                                                                                                <  3:81718/141‑1 (MQ=255)
                         cTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTg                                                                                                                >  4:81718/1‑141 (MQ=255)
                           gagaTTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGgt                                                                                                       >  1:145488/1‑148 (MQ=255)
                           gagaTTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGgt                                                                                                       >  4:816316/1‑148 (MQ=255)
                               ttCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGt                                                                                                                         >  3:477098/1‑126 (MQ=255)
                               ttCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGt                                                                                                                         <  4:477098/126‑1 (MQ=255)
                                tCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                           >  1:843058/1‑107 (MQ=255)
                                tCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                           <  2:843058/107‑1 (MQ=255)
                                      tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACATGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                           <  2:328295/101‑1 (MQ=255)
                                      tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGa                                                                                                                     <  1:955416/123‑1 (MQ=255)
                                      tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGa                                                                                                                     >  2:955416/1‑123 (MQ=255)
                                      tAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGAGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACgg                                                                                                                                           >  1:328295/1‑101 (MQ=255)
                                       aaCGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAAttt                                                                                           <  2:145488/148‑1 (MQ=255)
                                         cgacgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACa                                                                                                           >  2:1097846/1‑130 (MQ=255)
                                         cgacgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACa                                                                                                           <  1:1097846/130‑1 (MQ=255)
                                         cgacgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATttgt                                                                                         >  3:304645/1‑148 (MQ=255)
                                           acgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAcc                                                                                                                  >  3:179153/1‑121 (MQ=255)
                                           acgCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAcc                                                                                                                  <  4:179153/121‑1 (MQ=255)
                                                   gggTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTtaataa                                                                                               <  1:572447/132‑1 (MQ=255)
                                                   gggTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTtaataa                                                                                               >  2:572447/1‑132 (MQ=255)
                                                        tCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTaaa                                                                          <  4:304645/148‑1 (MQ=255)
                                                         ccAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAc                                                                                                                   <  3:147793/106‑1 (MQ=255)
                                                         ccAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAAc                                                                                                                   >  4:147793/1‑106 (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                                                                                                                            >  2:1127569/1‑87 (MQ=255)
                                                                   cTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGAt                                                                                                                            <  1:1127569/87‑1 (MQ=255)
                                                                    ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                                 <  2:815852/129‑1 (MQ=255)
                                                                    ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                                 >  3:750086/1‑129 (MQ=255)
                                                                    ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                                 <  4:750086/129‑1 (MQ=255)
                                                                    ttCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACa                                                                                 >  1:815852/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                                                                                                          <  2:209402/85‑1 (MQ=255)
                                                                                       ggTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACAt                                                                                                          >  1:209402/1‑85 (MQ=255)
                                                                                         tAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGc                                                           >  1:285214/1‑130 (MQ=255)
                                                                                         tAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGc                                                           <  2:285214/130‑1 (MQ=255)
                                                                                              ccAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATttgtt                                                                                        >  1:793576/1‑96 (MQ=255)
                                                                                              ccAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATttgtt                                                                                        <  2:793576/96‑1 (MQ=255)
                                                                                                aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGa                                                                                                                             >  4:78462/1‑57 (MQ=255)
                                                                                                aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGa                                                                                                                             <  3:78462/57‑1 (MQ=255)
                                                                                                aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATgg                                                                                                        <  1:471342/78‑1 (MQ=255)
                                                                                                aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATgg                                                                                                        >  2:471342/1‑78 (MQ=255)
                                                                                                aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgt                                                     <  2:284550/129‑1 (MQ=255)
                                                                                                aGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgt                                                     >  1:284550/1‑129 (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                    <  2:156279/146‑1 (MQ=255)
                                                                                                                atcatcatcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGa                    >  1:156279/1‑146 (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                                                     >  2:654613/1‑105 (MQ=255)
                                                                                                                        catcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgt                                                     <  1:654613/105‑1 (MQ=255)
                                                                                                                            tGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCgg                         <  4:139590/129‑1 (MQ=255)
                                                                                                                            tGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCgg                         >  3:139590/1‑129 (MQ=255)
                                                                                                                              aTGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGa    >  3:129313/1‑148 (MQ=255)
                                                                                                                                gCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAg  <  4:129313/148‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                               
GGTTTTCAACCACGCCCTGATAAACCTGAGATTCGACATAACGACGCAGCAGGGTATCCAGCAACGCCTTCGGATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAG  >  NC_000913/3917560‑3917835

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: